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PROTEIN DESIGN AUTOMATION FOR PROTEIN LIBRARIES 

This application Is a continuing application of Serial No. 09/927,790. filed on August 10, 2001 and 
claims the benefit of the filing dates of Serial Nos. 60/31 1 ,545, filed on August 1 0, 2001 , 60/324,899, 
filed on September 25, 2001, 60/351,937, filed on January 25, 2002, and 60/352,103, filed on January 
25, 2002. 

FIELD OF THE INVENTION 

The invention relates to the use of a variety of computation methods, including protein design 
automation (PDA^'^) technobgy to generate computationally prescreened secondary libraries of 
proteins, and to methods of making and methods and compositions utilizing the libraries. 

BACKGROUND OF THE INVENTION 

Directed molecular evolution may be used to create proteins and enzymes with novel functions and 
properties. Starting with a known natural protein, several rounds of mutagenesis, functional 
screening, and/or selection and propagation of successful sequences are performed. The advantage 
of this process is that it may be used to rapidly evolve any protein without knowledge of its structure. 
Several different mutagenesis strategies exist, including point mutagenesis by error-prone PGR, 
cassette mutagenesis, and DNA shuffling. These techniques have had many successes; however, 
they are all handicapped by their inability to produce more than a tiny fraction of the potential changes 
and their ability to effectively explore all possible sequences. For example, there are 20^°° possible 
amino ackJ changes for an average protein approximately 500 amino acids long. Clearly, the 
mutagenesis and functional screening of so many mutants is impossible; directed evolution provides a 
very sparse sampling of the possible sequences and hence examines only a small portion of possible 
improved proteins, typically point mutants or recombinations of existing sequences. By sampling 
randomly from the vast number of possible sequences, directed evolution is unbiased and broadly 
applicable, but inherently inefficient because it ignores all structural and biophysical knowledge of 
proteins. 
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In contrast, compufational nnethods may be used to screen enomnous sequence libraries (up to or 
more than 10" in a single calculation) overcoming the key limitation of experimental library screen! ng 
methods such as directed molecular evolution. There are a wide variety of methods known for 
generating and evaluating sequences. These include, but are not limited to, sequence profBing 
(Bowie and Bsenberg, Science 253(5016): 164-70. (1991)), rotamer library selectfons (Dahiyatand 
Mayo. Protein Sci 5(5): 895-903 (1996); Dahiyat and Mayo, Science 278(5335): 82-7 (1997); 
Desjartais and Handel. Protein Science 4: 2006-2018 (1995); Harbury et al, PNAS USA 92(18): 8408- 
8412 (1995); Kono etal., Proteins: Structure. Function and Genetics 19: 244-255 (1994); Hellinga and 
Richards. PNAS USA 91: 5803-5807 (1994)); and residue pair potentials (Jones, Protein Science 3: 
567-574, (1994)). (see Altschul and Koonin, Trends Biochem Sd 23(11): 444-447. (1998); (see 
Altschul et al.. J. Mol. Biol. 215(3): 403 (1990) and Lockless and Ranganathan, Science 286:295-299 
(1999). Pattern discovery in Biomolecular Data: Tools, Techniques, and Applications; edited by 
Jason T.L. Vteng, Bruce A Shapiro. Dennis Shasha. New Yori<: Oxford University. 1999.) 

Directed evolutwn is a random technique. Currently, there is no comprehensive rational design 
approach that allows eflkaent e}q3toratfon of all possible sequence space. 



SUMMARY OF THE INVENTION 

The present invention provides methods for generating a secondary library of scaffold protein variants 
comprising provkJing a primary library comprising a rank-ordered list or filtered set of scaJfold protein 
primary variant sequences. A list of primary variant poslttons in the primary library is then generated, 
and a pluraUty of the primary variant positions Is then combined to generate a secondary library of 
secondary sequences. 



It is an object of the present invention to provide computattonal methods for prescreening sequence 
libraries to generate and select secondary libraries, which may then be made and evaluated 
experimentally. 



In an additbnal object, the invention provides methods for generating a secondary library of scaffold 
protein variants comprising providing a primary library comprising a rank-ordered list or filtered set of 
scaffold protein primary variant sequences, and generating a probability distribution of amino acM 
residues in a plurality of variant positrons. The pluraBty of the amino add residues is combined to 
generate a secondary library of secondary sequences. These sequences may then be optionally 
synthesized and tested, in a variety of ways, induding multiplexing PGR with pooled oligonucleotWes. 
error prone PCR. gene shuffling, etc. 
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In a further object, the invention provides compositions comprising a plurality of secondary variant 
proteins or nucleic acids encoding the proteins, wherein the plurality comprises all or a subset of the 
secondary library. The invention further provides cells comprising the library, particularly mammalian 
cells. 

In an additional object, the invention provides methods for generating a secondary library of scaffold 
protein variants comprising providing a first library rank-ordered list or filtered set of scaffold protei n 
primary variants, generating a probability distribution of amino acid residues In a plurality of variant 
positions; and synthesizing a plurality of scaffold protein secondary variants comprising a plurality of 
the amino acid residues to form a secondary library. At least one of the secondary variants is 
different from the primary variants. 

it Is a further object of the invention to provide a method for receiving a scaffold protein structure with 
residue positions; selecting a collection of variable residue positions from said residue positions; 
establishing a group of potential rotamers for each of said variable residue positions, and wherein a 
first group for a first variable residue position has a first set of rotamers from at least two 
different amino acid side chains, and wherein a second group for a second variable residue position 
has a second set of rotamers from at least two different amino acid side chains; and, analyzing the 
Interaction of each of said rotamers in each group with all or part of the remainder of said protein to 
generate a set of optimized protein sequences. 

It is a further object of the invention to provide a method for receiving a scaffold protein with residue 
positions; selecting a collection of variable residue positions from said residue positions; establishing 
a group of potential amino acids for each of said variable residue positions, wherein a first group for a 
first variable residue position has a first set of at least two amino acid side chains, and wherein a 
second group for a second variable residue position has a second set of at least two different amino 
acid side chains; and, analyzing the interaction of each of said amino acids with all or part of the 
remainder of said protein to generate a set of optimized protein sequences. 

It is a further object of the invention to provide a method for receiving a scaffold protein with residue 
positions; selecting a set of variable residue positions from said residue positions; establishing a 
group of potential rotamers for each of said variable residue positions; analyzing the interaction of 
each of said rotamers with alt or part of the remainder of said protein to generate a set of optimized 
protein sequences, wherein said analyzing step includes the use of at least one scoring function; and, 
generating a library of said optimized protein sequences. 

It is a further object of the invention to provide a method for receiving a scaffold protein with residue 
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positions; selecting a set of variable residue positions from said residue positions; classifying each 
variable residue position as either a core, surfece or boundary position; establishing a group of 
potential amino acids for each of said variable residue positions, wherein the group for at least one 
variable residue position has at least two different amino acid side chains; and, analyzing the 
interaction of each of said amino acids with all or part of the remainder of said protein to generate a 
set of optimized protein sequences, wherein said analyzing step includes the use of at least one 
scoring function. 

It is a further object of the invention to provide a method for receiving a scaffold protein witii residue 
positions; selecting a set of variable residue positions from said residue positions; establishing a 
group of potential rotamers for each of said variable residue positions, wherein the group for at least 
one variable residue position has rotamers of at least two different amino acid side chains, and 
wherein at least one of said amino acid side chains Is from a hydrophlllc amino add and. analyzing 
the interaction of each of said rotamers with all or part of the remainder of said protein to generate a 
set of opUmized protein sequences, wherein said analyzing step includes the use of at least one 
scoring function. 



It is a further object of the invention to provide a computational method for receiving a scaffold protein 
with residue positions; selecting a collection of variable residue positions from said residue positions; 
providing a sequence alignment of a plurality of related proteins; generating a frequency of 
occurrence for individual amino acids in at least a plurality of positions with said alignments; creating 
a pseudo-energy scoring function using said frequencies; using said pseudo-energy scoring function 
and at least one additional scoring function to generate a set of optimized protein sequences. 

It is a further object of the invention to provide a computational method comprising receiving a 
scaffold protein with residue positions; selecting a collection of variable residue positions from said 
residue positions; providing a sequence alignment of a plurality of related proteins; generating a 
frequency of occurrence for individual amino acids in at least a plurality of positions with said proteins; 
selecting a group of potential amino acids for each of said variable residue positions, wherein a first 
group for a first variable residue position has a first set of at least two amino acid side chains, and 
wherein a second group for a second variable residue position has a second set of at least two 
different amino acid side chains according to their frequency of occurrence; and. analyzing the 
interaction of each of said amino acids at each variable residue position with all or part of the 
remainder of said protein using at least one scoring function to generate a set of optimized protein 
sequences. 



It is a further object of the Invention to provide a method computational method for receiving a scaffold 
protein with residue positions; selecting a collection of variable residue positions from said residue 
positions; providing an amino acid substitution matrix; creating a pseudo-energy scoring function 
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using said matrix; using said pseudo-energy scoring function and at least one additional scoring 
function to generate a set of optimized protein sequences. 

It is a further object of the invention to provide a method for receiving a scaffold protein with residue 
positions; selecting a collection of at least one variable residue position from said residue positions; 
importing a set of coordinates for a scaffold protein, said scaffold protein comprising amino add 
positions; analyzing the interaction of each of said amino adds with all or part of the remainder of said 
protein; utilizing a plurality of scoring functions, at least a first a scoring function having a first weight 
and a second scoring function having a second weight, to generate at least one variable decoy 
sequence; and, comparing the scores from said scoring functions of said variable decoy sequence to 
the scores of a reference state to generate modified weights, wherein each weight is increased if the 
corresponding score of the decoy is higher than the corresponding score of the reference state and 
each weight is decreased if the corresponding score of the decoy is lower than the corresponding 
score of the reference state and, wherein the extent of increase or decrease is based on the relative 
individual and total scores of the decoy and reference states. 

It is a further object of the invention to provide a method for receiving a scaffold protein with residue 
positions; selecting a collection of variable residue positions from said residue positions; importing a 
set of coordinates for a scaffold protein, said scaffold protein comprising amino acid positions; 
generating a variable protein sequence comprising a defined energy state for each amino acid 
position; applying an energy Increase to at least one of said defined energy states for a least one of 
said amino acid positions; and, generating at least one alternate variable protein sequence. 

it is a further object of the invention to provide a method for receiving a scaffold protein with residue 
positions; selecting a collection of variable residue positions from said residue positions; importing a 
set of coordinates for a scaffold protein, said scaffold protein comprising amino acid positions; 
generating a variable protein sequence comprising a defined energy state for each amino acid 
position; applying a probability parameter to at least one of said amino acid positions; and generating 
at least one alternate variable protein sequence. 

It is a further object of the invention to provide a method for receiving a scaffold protein with residue 
positions; selecting a collection of variable residue positions from said residue positions; importing a 
set of coordinates for a scaffold protein, said scaffold protein comprising amino acid positions; 
generating a set of optimized variant protein sequences comprising one or more variant amino acids; 
and, applying a clustering algorithm to cluster said set into a plurality of subsets. 

It is a further object of the invention to provide a method for receiving at least one scaffold protein 
structure with variable residue positions of a target protein; computationally generating a set of 
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primary variant amino acid sequences that adopt a conformation similar to the conformation of said 
target protein; and. Identifying at least one protein sequence that is similar to at least one member of 
said set of primary variants, but is dissimilar to said target protein amino acid sequence. 

It is an additional object of the Invention to provide a method for generating variant protein sequence 
libraries comprising providing populations of at least two double stranded donor fragments 
corresponding to a nucleic acid template; adding polymerase primers capable of hybridizing to end 
regions of each of said population of donor fragments; generating a population of hybrid double 
stranded molecules wherein one strand comprises a 5'-purification tag and the other strand comprises 
a 6'-phosphorylated overhang; enriching for variant strands by removing strands comprising a 5'-biotin 
moiety; annealing said variant strands to form at least two double stranded ligation substrates; and, 
ligating said ligation substrates to form a double stranded ligation product wherein said ligation 
product encodes a variant protein. 



These and other objects of the invention are to provide computational protein design and optimization 
techniques via an objective, quantitative design technique implemented in connection with a general 
purpose computer. 



BRIEF DESCRIPTION OF THE DRA\AnNGS 

Figure 1 depicts a gene assembly scheme. 

Figure 2 illustrates that most protein design simulations do not sufficiently may sequence space. As 
shown in the upper graph, most protein design simulations only map the lowest energy basin; thereby 
omitting other low energy basins that could provide viable sequences for computationally generated 
protein sequences. 

Figure 3 illustrates the point that the alternate low energy basins can represent equally good 
sequences for incorporation into a protein template. This is because the force field representation of 
the energy (i.e.. Ecaic) is not necessarily Identical to the actual energy (i.e., Etrue) associated with a 
native protein structure. 

Figure 4 illustrates the application of taboo for mapping sequence space. The calculated energy 
surface is manipulated based on previous solutions to discourage repeated convergence to the same 
local minimum. 

Figure 5 Illustrates clustering algorithms that may be used in the methods of the present invention. 

Figure 6 depicts an example of energy matrix clustering of designed \AA/Vdomaffi proteins using a 
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single linkage clustering algorithm. • 

Figure 7 depicts the data used to generate Figure 7. 

Figure 8 depicts representative structures from cluster 1 , 3, and 9. 

Figure 9 depicts an example of energy matrix clustering of designed SH3 proteins. 

Figure 10 depicts the superfamily of sequences designed for SH3. As shown In Figure 6, the virtual 
superiiamily of sequences designed using an SH3 backbone structure have significant homology to 
the template sequence and other members of the natural SH3 family. Identities with the native 
sequence are highlighted In dark grey. Functbnal positions are shaded in light grey. Note that 
although the simulations did not include a functional constraint, the nath/e functional residue usually 
appears with low frequency in the alignment 

Figure 11 illustrates coupling pattems in SH3 subfamilies. Interaction-based clustering reveals a 
series of virtual sequence subfamilies that contain various combinations of coupled amino acids 
(highlighted in different shades of grey. Note that some subfamilies differ by amino acids coupled at 7 
positions (medium intensity shading). The amino ack) couplings lead to multiple low energy solutions 
in different sequence subspaces. As a result, some subfamilies have more similarity to the wild type 
sequence than others. 

Figure 12 depicts the synthesis of a full-length gene and all possible mutations by PGR. Overiapping 
oligonucleotides corresponding to the full-length gene (black bar, Step 1) are synthesized, heated and 
annealed. Addition of Pfu DNA polymerase to the annealed oligonucleotides results in the 6' 3' 
synthesis of DNA (Step 2) to produce longer DNA fragments (Step 3). Repeated cycles of heating, 
annealing (Step 4) results in the production of longer DNA, including some full-length molecules. 
These may be selected by a second round of PGR using primers (an^owed) corresponding to the end 
of the full-length gene (Step 5). 

Figure 13 depicts the reduction of the dimensionality of sequence space by PDA™ technology 
screening. From left to right, 1: without PDA"^*^ technology; 2: without PDA'^'^ technology not counting 
Cysteine, Proline, Glycine; 3: with PDA^^ technology using the 1% criterion, modeling free enzyme; 4: 
with PDA™ technology using the 1% criterion, modeling enzyme-substrate complex; 5: with PDA™ 
technology using the 5% criterion nwdeling free enzyme; 6: with PDA^"^ technology using the 5% 
criterion modeling enzyme-substrate complex. 
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Rgure 14 depicts the active site of B. circulans xylanase. Those positions included in the PDA™ 
technology design are shown by their side chain representation. 

Figure 15 depots cefotaxime resistance of E. coli expressing wBd-type (WT) and PDA™ technology. 



Figure 16 depicts a prefened scheme for synthesizing a library of the invention. The wild-type gene, 
or any starting gene, such as the gene for the gtobal minima gene, may be used. Oligonucleotides 
comprising different amino acids at the different variant positions may be used during PGR using 
standard primers. This generally requires fewer oligonucleotides and may result in fewer en-ors. 



Rgure 1 7 depicts and overtapping extension method. At the top of Figure 6 is the template DNA 
showing the locations of the regions to be mutated (black boxes) and the binding sites of the relevant 
primers (arrows). The primers R1 and R2 represent a pool of primers, each containing a different 
mutation; as described herein, this may be done using different ratios of primers if desired. The 
variant position is flanked by regions of homology sufficient to get hybridization. In this example, 
three separate PGR reactions are done for step 1. The first reactbn contains the template plus oligos 
F1 and R1 . The second reaction contains template plus F2 and R2. and the third contains the 
temptete and F3 and R3. The reaction pnxlucts are shown. In Step 2. the products from Step 1 tube 
1 and Step 1 tube 2 are taken. After purification away from the primers, these are added to a fresh 
PGR reaction together with F1 and R4. During the Denaturation phase of the PGR. the overtapping 
regions anneal and the second strand is synthesized. The product is then amplified by the outside 
primers. In Step 3, the purified product ft-om Step 2 is used in a third PGR reaction, together with the 
product of Step 1 , tube 3 and the primers F1 and R3. The final product corresponds to the full-length 
gene and contains the required mutations. 



Figure 18 depicts a ligation of PGR reaction products to synthesize the libraries of the invention. In 
this technique, the primers also contain an endonuclease restriction site (RE), either bfunt. 5' 
overtianging or 3' overtianging. We set up three separate PGR reactions for Step 1. The first 
reaction contains the template plus oligos F1 and R1. The second reaction contains template plus F2 
and R2. and the third contains the template and F3 and R3. The reaction products are shown. In 
Step 2. the products of step 1 are purified and then digested with the appropriate restriction 
endonuclease. The digestion products from Step 2. tube 1 and Step 2. tube 2 and Igate them 
together with DNA ligase (step 3). The products are then amplified in Step 4 using primer F1 and R4. 
The whole process is then repeated by digesting the amplified products, ligating them to the digested 
products of Step 2, tube 3, and then amplifying the final product by primers F1 and R3. It would also 
be possibte to ligate all three PGR products from Step 1 together in one reaction, providing the two 
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restriction sites (RE1 and RE2) were different 

Figure 19 depicts blunt end ligation of PGR products. In this technique, the primers such as F1 and 
R1 do not overlap, but they abut Again three separate PGR reactions are performed. The products 
from tube 1 and tube 2 are ligated, and then amplified with outside primers F1 and R4. This product 
is then ligated with the product from Step 1, tube 3. The final products are then amplified with primers 
F1 and R3. 

Figure 20A and B depicts M13 single stranded template production of mutated PGR products. 
Primerl and Primer2 (each representing a pool of primers corresponding to desired mutations) are 
mixed with the Ml 3 template containing the wild type gene or any starting gene. PGR produces the 
desired product (11) containing the combinations of the desired mutations incorporated in Primerl 
and Primer2. This scheme may be used to produce a gene with mutations, or fragments of a gene 
with mutations that are then linked together via ligation or PGR for example. 

Figure 21 A-E depict examples of some preferred combinations. 

DETAILED DESCRIPTION OF THE INVENTION 
As used herein, the following terms shall have the meaning as described below. 

By "altered phenotvpe" or "changed physiology" or other grammatical equivalents herein is meant that 
the phenotype of the cell containing a variable amino acid sequence (preferably an optimized 
sequence) is altered in some way, preferably in some detectable, observable and/or measurable way. 
Examples of phenotypic changes include, but are not limited to: gross physical changes such as 
changes in cell morphology, cell growth, cell viability, adhesion to substrates or other cells, and 
cellular density; changes in the expression of one or more RNAs, proteins, lipids, hormones, 
cytokines, or other molecules; changes In the equilibrium state (i.e. half -life) or one or more RNAs, 
proteins, lipids, hormones, cytokines, or other molecules; changes in the localization of one or more 
RNAs, proteins, lipids, hormones, cytokines, or other molecules; changes in the bioactivlty or specific 
activity of one or more RNAs, proteins, lipids, hormones, cytokines, receptors, or other molecules; 
changes In phosphorylation; changes in the secretion of ions, cytokines, hormones, growth factors, or 
other molecules; alterations in cellular membrane potential, polarization, integrity or transport; 
changes in infectlvity. susceptibility, latency, adhesion, and uptake of viruses and bacterial pathogens; 
etc. By "capable of altering the phenotype" herein is meant that the library member (e.g. the variable 
amino acid sequence and/or the variable nucleic acid sequence) may change the phenotype of the 
cell in some detectable and/or measurable way. 



9 



wo 03/014325 



PCT/US02/25588 



By "alternate amino acid" as used herein is meant an amino acid slate that differs from the amino acid 
defined by the starting an^no acid sequence in the protein design cycle. As outlined below, this 
starting amino add sequence (e.g, the scaffold protein) may be a wild-type sequence or a variant 
sequence. 

By ."amino add identity' as used herein is meant the identity of an amino add at a specified position; 
e.g. when the position of an amino add is specified, which one of the 20 naturally occurring or non- 
natural analogs is present at that position. 

By "boundary residues" as used herein is meant, residue positions that are not clearly in the protein 
core or on the protein surface. Methods for determining boundary residues are outlined below. The 
solvent accessibility of side chains in boundary positions is determined by the confbmiation and 
identities of the residues surrounding it. in a preferred embodiment, both hydrophobic and polar 
amino adds can be considered as possible replacement residues at boundary positions. 

"candidate bioactive aqenf ' or "candidate drugs" or grammatical equivalents herein Is meant any 
molecule, e.g. proteins (which herein indudes proteins, polypeptides, and peptides), small organic or 
inorganic molecules, polysaccharides, polynudeotides, etc. which are to be tested against a particular 
target Candidate agents encompass numerous chemical classes. In a preferred embodiment, the 
candidate agents are organic molecules, particularly small organic molecules, comprising functional 
groups necessary for structural interaction with proteins, partlculariy hydrogen bonding, and typically 
indude at least an amine, carbonyl, hydroxyl or carboxyl group, preferably at least two of the 
functional chemical groups. The candidate agents often comprise cyclical carbon or heterocydic 
structures and/or aromatic or polyaromatic structures substituted with one or more chemical functional 
groups. A prefenred embodiment Is a protein where the uses include therapeutic, veterinary, 
agricultural, and industrial applications. 

By a "cellular library" herein is meant a plurality of cells wherein generally each cell within the library 
contains at least one member of the library. Ideally each cell contains a single and different library 
member, although as will be appreciated by those in the art. some cells within the library may not 
contain a library nnember and sonne may contain more than one library member. When methods 
other than retroviral infection are used to introduce the library members into a plurality of cells, the 
distribution of library members within the individual cell members of the cellular library may vary 
widely, as it is generally difficult to control the number of nucleic adds which enter a cell during 
electroporation and other transformation methods. Suitable cell types for cellular libraries are induded 
below. In addition, as will be appredated by those in the art. a cellular fibrary generally includes a 
single cell type, although in some embodiments, a cellular library may contain two or more ceil types. 
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By "chemically modified" as used herein is meant to Include modification via chemical reactions as 
well as enzynnatic reactions. The substrates in these reactions generally Include, but are not limited 
to, alkyl groups (including but not limited to straight and branched alkanes, alkenes, and alkynes), aryi 
groups (including but not limited to arenes and heteroaryl), alcohols, ethers, amines, aldehydes, 
ketones, carboxylic ackJs, esters, amides, heterocyclic compounds (Including, but not limited to, 
piperidines, pyrrolidines, purines, pyrimidines, benzodiazepins, and carbohydrates), steroids 
(including but not limited to estrogens, androgens, cortisone, ecodysone, etc.), secondary metabolites 
(including, but not limited to, terpenoids, alkaloids, polyketides, beta-lactams, polyether antibiotics, 
and aminoglycosides), organometallic compounds, lipids, amino acids, and nucleosides. The 
reactions generally include, but are not limited to, hydrolysis, reduction, oxidation, alkylation, aromatic 
substitutions, electrocydizations, dipolar cyclizations, radical anion, radical cation, metal mediated 
couplings, and polymerization. 

By "clustering algorithm" herein is meant an algorithm that may be used to separate a large selection 
or set of computationally generated sequences into subsets that represent various sub-regions of 
sequence space. Clustering algorithms are well known in the art, and representative examples are 
outlined below. 

By "control sequences" or "regulatory sequences" as used herein refers to DNA sequences necessary 
for the expression of a gene in a particular host organism. The control sequences that are suitable for 
prokaryotes, for example, include a promoter, optionally an operator sequence, and a ribosome 
binding site. Eukaryotic cells utilize control sequences including, but not limited to, promoters, 
polyadenylation signals, and enhancers. 

By "core positions" as used herein is meant, positions that are in the interior of a protein or which are 
inaccessible or nearly inaccessible to solvent. Methods for determining which position comprise core 
positions are outiined below. As more fully outlined bebw, in a preferred embodiment, for design 
purposes, only hydrophobic amino acids are considered for incorporation into variable positions at 
core variable positions. As more fully outiined bebw, in an alternate preferred embodiment, polar 
amino acids are considered at core positions only if they form favorable electrostatic or hydrogen 
bond interactions with other polar groups, or If disruption of the scaffold is desired. 

By " coupling" as used herein is meant the non-additive contribution (e.g. synergistic) of two or more 
amino acids to an interaction involving said amino acids. Coupling can be positive (the interaction is 
more favorable than the sum of the individual contributions), neutral, or negative (the interaction Is 
less favorable than the sum of tiie individual contributions). Such coupling typically occurs for amino 
acids located very close in space. 
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^ "decoy state." -decoy stnicture." or °decov seq .i,.nrp- as used herein is meant a protein sequence 
and structure that is different from a specified reference state, and that seroes as a comparison state 
for use in various parameter optimization methods. Decoy structures are more fully described betow. 

By "donor fraqmenf or "donor nudeic acid fragmenf as used herein is meant nucleic acid fragments 
generated from or corresponding to a template nucleic acid molecule. Preferably, the donor 
fragments are generated using nrradified primers and a polymerase, although fragments may be 
generated using enzymatic, chemical or physical cleavage (e.g. shearing) of template nucleic aad 
molecules. Any DN/VRNA polymerase is suitable; however thermophilic polymerases are preferred. 

An -enerqy matrix" is defined for the present purposes as fbltows. A protein design cycle simulation is 
performed to yield a single protein sequence/structure. In the context of this state, all amino adds (in 
all rotemer states) are sampled at each position or at each variable position. Alternatively, less than 
all rotamer states, or less than all amino adds, are sampled at some or all of the positions. Suite bte 
sampling techniques to generate the energies are outlined herein. The context-dependent energy of 
each amino acid is stored. An energy matrix is defined by the listing of the context-dependent energy 
of each amino add at each position of the stmcture. The similarity of two energy matrices (from two 
different simulations) may be defined as the root-mean-squared^leviation of two energy matrices. It 
should be noted that In some cases, energy matrices comprising less than all of the possibte 
interactions can be constructed. 

By " filtered set " herein is meant the optimized protein sequences that are generated using some sort 
of selection criteria. Although in some cases, the set may comprise an arbitrary or random selection 
of a subset of the primary sequences. In a preferred embodiment, the filtered set comprises a rank 
ordered list of sequences. As ouffined herein, this may be done in a variety of ways, induding an 
artjitrary cutoff (for example, the top 10,000 sequences are chosen, or the top 1000 and the bottom 
1000), an energy limitation (e.g. anything with a total energy calculation below X). or when a certain 
number of residue positions have been varied (e.g. the set is complete when 10 variable positfons is 
adiieved. etc). As is outlined more fully below, filtering can be used as all or part of the primary, 
secondary, tertiary, ete. library generation; that is. filtering can be the sole computational analysis or 
part of a larger analysis, at one or more of the steps of the invention. For example, a primary library 
may be computationally generated using PDA and a filtering step applied to define the set for 
secondary library generation, eto. 

By I fixed position' herein is meant, residue posifrons at which the amino acid identity vtrill be held . 
constant in a protein design calculation. In some embodiments, fixed positions may be floated, as 
defined betaw. That is. in some embodimente. an amino acid identity is kept fixed, but its rotemeric 
state is allowed to change. In other embodiments, the amino add identity and rotemeric state are 
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held constant The conformation and amino acid identity may be that observed in the scaffold 
structure or the conformation and/or amino acid identity may be different than that obsen/ed in the 
scaffold structure. 

By "floated position" herein Is meant, a position at which the amino acid confomnatfon but not the 
amino add identity is allowed to vary in a protein design calculation. The floated position may be 
fixed as a non-wild type residue. For example, when known site-directed mutagenesis techniques 
have shown that a particular amino acid is desirable (for example, to eliminate a proteolytic site or 
alter the substrate specificity of an enzyme), the position may be constrained to allow only that amino 
acid. Alternatively, the methods of the present invention may be used to evaluate specific mutations 
de novo. 

By " gene assembly procedures" as used herein is meant either enzymatic or chemical methods of 
joining gene fragments. A wide variety of exemplary methods are included herein and described 
below. 

By " global optimum protein sequence" as used herein is meant an amino acid sequence that best fits 
the mathematical equations of the computational process. As will be appreciated by those in the art, 
a global optimum sequence is the sequence that has the lowest energy or best score of any possible 
sequence In the context of the particular computational analysis utilized . That Is. the global optimum 
sequence depends on the scoring or ranking systems used, and may change with different 
computational parameters. For example, when PDA'^'^ is used, the global optimum will depend on the 
scoring functions utilized, the weighting factors, etc. In addition, tiiere are any number of sequences 
that are not the global minimum but that have low energies or favorable scores referred to herein as 
"optimized sequences", defined below. 

By "labeled" herein is meant that nucleic acids, proteins, candidate agents, antibodies or other 
components of the Invention have at least one element, isotope, or chemical compound attached to 
enable the detection of nucleic acids, proteins and antibodies of the invention. 

By " ligation product " as used herein is meant either tiie single stranded or double stranded nucleic 
acid molecule resulting when at least two ligation substrates are ligated togetiier. 

By " ligation substrate" as used herein is meant either a single or double stranded nucleic acid 
molecule formed by annealing from two complementary donor fragments in which one donor fragment 
has a 5'-phosphorylated overhang and the other firagment has a free 3'-termlnus (see Figure 1). 

By " nucleic acid template" herein is meant a single or double stranded nucleic acid. In a preferred 
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embodiment, the nucleic add template is used to generate donor fragments, defined above. The 
donor fragments may be obtained directly from the nucleic acid template or separately obtained, e.g.. 
by nucleic acid synthesis, fragmentation (e.g. enzymatic, chemical or physical) or amplification 
reactions. A nucleic acid template may comprise an intact gene, or a fragment of a gene encoding 
functional domains of a protein, such as enzymatic domains, regulatory sequences, binding domains, 
etc., as well as smaller gene fragments The template nucleic acid may be from any organism, either 
prokaryotic or eukaryotic. The template sequence may be naturally occuning. a variant, a product of 
a computational step, etc. 

By "nucleoside*; as used herein, includes nucleotides, nucleosides and analogs, including modified 
nucleosides such as amino modified nucleosides and includes non-roturally occurring analog 
stmctures, I.e. the individual units of a peptide nucleic acid, each containing a base, are refen-ed to 
herein as a nucleoside. 

By 'bperablv linked" as used herein means two or more nucleic acids linked together such that the 
desired functionality is achieved. For example, when a first nucleic acid sequence is placed into a 
functional relationship with another nucleic acid sequence. For example, DNA for a presequence or 
secretory leader is operably linked to DNA for a polypeptide if it is expressed as a preprotein that 
participates in the secretion of the polypeptide; a promoter or enhancer is operably linked to a coding 
sequence if it affects the transcriptton of the sequence; or a ribosome binding site is operably linked to 
a coding sequence if it is positioned so as to facilitate translation. Generally, operably linked DNA 
sequences are contiguous, and in the case of a secretory leader, contiguous and in reading phase. 
However, enhancers do not have to be contiguous. Linking can be accomplished by ligation at 
convenient restriction sites. If such sites do not exist, the synthetic oligonucleotide adaptors or linkers 
are used in accordance with conventional pracBce 

By "optimized protein sequence" as used herein is meant a sequence with at least one optimized 
property. For example, in the context of a particular computational analysis, an optimized sequence 
will exhibit a low energy or favorable score. For example, when PDA™ is used, an optimized 
sequence is one which has a lower energy than the energy of the starting scaffold protein. 
Alternatively, an optimized protein sequence may have one or more protein properties, defined below, 
that are desirably different as compared to the starting scaffold protein. An optimized protein 
sequence may or may not be the global optimum sequence, however, an optimized protein sequence 
has at least one amino acid substitution, insertion or deletion as compared to the starting scaffold 
protein used to generate the optimized sequence. 

By a "plurality ofceltel herein is meant roughly from about 10^ cells to 10^, 10® or 10^ witfi from 10^ to 
10° being preferred. 
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By "position" as used herein Is meant a location in the sequence of a protein. Positions are typically 
numbered using the protein numbering scheme described below. In the context of a given scaffold 
protein, each position is associated with the location and/or orientation of its associated backbone 
atoms in three dimensions. Consequently, positions may be described by their secondary structure 
and by whether an amino acid located at that position would be solvent exposed or buried in the 
protein core. 

By "presentation scaffold" or "presentation structure" as used herein is meant a protein structure that 
allows the scaffold protein, generally a peptide, to take on a certain conformation. For example, there 
are a wide variety of "ministructures" known, sometimes refen^ed to as "presentation structures", that 
can confer conformational stability or give a random sequence a conformationally restricted form. 
Proteins interact with each other largely through conformationally constrained domains. Although 
small peptides with fireely rotating amino and carboxyl termini can have potent functions as is known 
in the art, the conversion of such peptide structures into pharmacologic agents is difficult due to the 
inability to predict side-chain positions for peptidomimetic synthesis. Therefore the presentation of 
peptides in conformationally constrained structures will benefit both the later generation of 
pharmaceuticals and will also likely lead to higher affinity interactions of tiie peptide with tiie target 
protein. This fact has been recognized in the combinatorial library generation systems using 
biologically generated short peptides in bacterial phage systems. A number of workers have 
constructed snrtall domain molecules in which one might present randomized peptide structures. 
Thus, synttietic presentation structures, i.e. artificial polypeptides, are capable of presenting a 
randomized peptide as a conformationally-restricted domain. In addition, random peptide structures 
that are not totally random, i.e., tfiat are selected or filtered as described herein may be presented. 
Preferred presentation structures maximize accessibility to the peptide by presenting It on an exterior 
loop. Accordingly, suitable presentation structures include, but are not limited to, minibody structures, 
loops on beta-sheet turns and coiled-coil stem stixictures in which residues not critical to structure are 
randomized, zinc-finger domains, cysteine-linked (disulfide) structures, transglutaminase linked 
sti-uctures, cyclic peptides, B-loop structures, helical banrels or bundles, leucine zipper motifs, etc. 

By "primary library* as used herein is meant a collection of sequences, preferably optimized and 
generally, but not always, in the form of a filtered set, a rank-ordered list (e.g. a scored or sampled 
set), an alignment, a probability distribution table, etc. A primary library is generated as a targeted 
subset of all or a portion of the sequence space for a particular scaffold protein. That is, a primary 
library is generated using any number of techniques, either alone or in combination; to reduce the size 
of the set of sequences likely to take on a particular fold or have a particular protein property. The 
primary library preferably comprises a set of sequences resulting from computation, which may 
include energy calculations and/or statistical or knowledge based approaches. In general, it is 
preferable to have the prinr^ry library be large enough to randomly sample a reasonable sequence 
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space to allow for robust secondary libraries. Thus, primary libraries that range ftom about 50 to 
about 10" are preferred, with from about 1000 to about 10' being parficulariy preferred, and from 
about 1000 to about 100,000 being especially preferred. 

By "probability parameter' as used herein is meant a parameter that governs the rate at which a given 
amino acid or rotemer state is sampled during a simulation. 

B y-protein' as used herein is meant at least two amino acids Hnked together by a peptide bond. As 
used herein, protein includes proteins, oligopeptides, polypeptides and peptides. The peptidyl group 
may comprise naturally occurring amino acids and peptide bonds, or synthetic pepfldomimetfc 
structures, i.e. "analogs", such as peptoids (see Simon et a!.. PNAS USA 89(20):9367 (1992)). The 
amino acids may either be naturally occumng or non-naturaliy occumng. The side chains may be in 
either the (R) or the (S) configuratton. In a preferred embodiment, the amino adds are in the (S) or 
L-conf^uration. 

"protein numbering scheme" herein is meant the manner in which, as is known in the art, the 
residues, or positions, of proteins are generally numbered. The residues, or positions, are generally 
sequentially numbered starting with the N-terminus of the protein. Thus a protein having a methionine 
at its N-tenrunus is said to have a methionine at residue or amino acid position 1 , with the next 
residues as 2. 3, 4, etc. In some embodiments, a set of aligned proteins is numbered together. In 
such cases, insertions relative to the consensus sequence are denoted by adding a letter after the 
number; for example, a one-residue insertion between positions 1 and 3 would produce the 
numbering 1, 2a, 2b, 3. Similarty. deletions relative to the consensus sequence are denoted by 
skipping a number; for example, a one residue deleflon between positions 1 and 3 would produce the 
numbering 1,3. 

By -protein properties" herein Is meant, biological, chemical, and physical properties including, but not 
limited to, enzymatk: activity, spedfkaty (including substrate specificity, kinetic assodation and 
dissociation rates, reaction mechanism, and pH profile), stability (induding thermal stabMly, stabiUty 
as a function of pH or solution conditions, resistance or susceptibility to ubiquitination or proteolytic 
degradation), solubility, aggregation, structural integrity, the creation of new antibody CDRs. generate 
new DNA RNA binding, generate peptide and peptidomimmetic libraries, crystallizability, binding 
affinity and specificity (to one or more molecules including proteins, nucleic acids, polysacdiarides, 
lipids, and small molecules), oligomerization state, dynamic properties (including conformational 
changes, altostery, con-elated motions, flexibility, rigidity, folding rate), subcellular tocalization. ability 
to be secreteid, ability to be displayed on the surface of a cell, posttranslational modification (including 
N- or C-finked glycosylaUon, Rpidation, and phosphorylation), ammenability to synthetic modification 
(including PEGylation, attachment to other molecules or surfaces), and ability to Induce altered 
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phenotype or changed physiology (including cytotoxic activity, Immunogenicity. toxicity, ability to 
signal, ability to stinnulate or inhibit celt proliferation, ability to Induce apoptosis, and ability to treat 
disease). As Is outlined herein, protein properties may be modulated using the techniques of the 
invention. When a biological activity is the property, modulation in this context includes both an 
increase or a decrease In activity. 

By "pseudo enenav" as used herein is meant an energy-lil<e term derived from non-energetic 
infomiation. Such pseudo energies are typically used as a mechanism for combining non-energetic 
information with energy based scoring ftinctions. For example, statistical information arising from 
structural analysis, sequence alignments, or simulation history may be incorporated into a calculation 
by their conversion to pseudo energies. 

By " recency parameter " as used herein means the application of at least one restraint to the most 
recent moves of a simulation (see Modem Heuristic Search Methods, edited by V.J. Rayward-Smith, 
et al., 1996, John Wley & Sons Ltd., hereby expressly incorporated by reference in its entirety). 

By "residue" as used herein is meant an amino acid side chain. A residue may be one of the naturally 
occurring amino acid side chains or a synthetic anatog. 

By "scaffold protein" herein is meant a protein for v\/hich a library of variants is desired. The scaffold 
protein is used as input in the protein design calculations, and often is used to facilitate experimental 
library generation. A scaffold protein may be any protein that has a known structure or for which a 
structure may be calculated, estimated, modeled or determined experimentally. As outlined more fully 
below, the scaffold protein may be a wild-type protein from any organism, a variant, a chimeric 
protein, etc. Preferred embodiments of scaffold proteins are outlined below. 

By "secondary library" as used herein is meant a library of amino acid sequences that is derived from 
a primary library using, a variety of approaches discussed further below, including both experimental 
and computational methods, or combinations thereof. Secondary libraries are generally generated 
experimentally and analyzed for the presence of members possessing desired protein properties. 
The secondary library may be either a subset of the primary library, or contain new Obrary members, 
i.e. sequences that are not found in the primary library. The secondary library typically comprises at 
least one member sequence that is not found in the primary library, and preferably a plurality of such 
sequences, although this Is not required. 

By " selectable gene." "selection gene" or "selectable marker" as used herein is meant any gene that 
enables survival and/or reproduction of the cells that express it. The marker gene may confer 
resistance to a selection agent such as an antibiotic, or may provide a protein required for growth. 



17 



wo 03/014325 



PCT/US02/25588 



By "sequence space" herem Is meant all sequential combinations of amino acids ttiat are possible for 
a defined protein and a defined set of positions thereof. For example, the sequence space for all 
positions of a 100-residue protein is 20^°°, and the sequence space for ten selected positions of a 
protein would be 20^°. if only the twenty naturally occuning amino acids are considered. 

By ."shuffling", as used herein means recombination of one or more protein, DNA, or RNA sequences. 
Shuffling may be done experimentally and/or computationally (e.g. "in silico shuffling"). See for 
example. U.S. Patent 6.319,714; WO 0042559WO 00/42560; and WO 00/42561. 

By "solid supporf or other grammatical equivalents herein is meant any material that may be modified 
to contain discrete individual sites appropriate for the attachment or association of beads, other solid 
support surfaces not in solution, and is amenable to at least one detection method. As will be 
appreciated by fliose in the art, the number of possible supports Is very large. Possible solid supports 
include, but are not limited to, glass and modified or functionalized glass, plastics (including acrylics, 
polystyrene and copolymers of styrene and other materials, polypropylene, polyethylene, 
polybutylene, polyurethanes. Teflon®, etc.). polysaccharides, nylon or nitroceflulose. resins, silica or 
silica-based materials including silicon and modified silicon, carbon, metals, inorganic glasses, 
plastics, optical fiber bundles, and a variety of other polymers. In general, tiie solid supports aflow 
optical detection and do not themselves appreciably fluoresce. 

By "sticky end" as used herein Is meant tiie end of an enzymatically cleaved DNA fragment thai has 
either a 5' or 3* overhang, and has tiie potential to interact favorably witti another sticky end with 
similar properties. 

By Isurface positions" as used herein is meant amino acid positions within a scaffokJ protein (or a 
variable protein) with a significant degree of solvent accessibility. Methods for the determination of 
surface positions are outiined below. In a prefenred embodiment, only polar amino acids are 
considered as possible replacement residues at surface positions in protein design calculations. 

By rtabu search algorithms" as used herein is meant any algorittims fi-om the class of searching 
metiiods in which searching moves are made such that moves already made, or made recently in tfie 
history of the search, are either avoided or disfavored. 

By "tertiary library" as used herein is meant a libraiy that is generated by computational or 
experimental modification or manipufation of a secondary library. 

By "variant protein sequence" as used herein is meant a protein sequence thai differs from anottier 
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protein sequence. In other words a variant protein sequence has at least one amino acid that differs 
from the amino acid defined by the starting amino add sequence in the protein design cycle. As 
outlined below, this starting amino acid sequence (e.g. the scaffold protein) may be a wild-type 
sequence or a variant sequence. 

By "variable residue position" herein is meant a position at which both the amino acid identity and 
confonration are allowed to be altered in a protein design calculation. The amino acid identity to 
which a position may be mutated may be the full set or a subset of the 20 naturally occumng amino 
acids or may be a set of non-naturally occurring amino acids or synthetic analogs. 

By "temperature factor'' as used herein is meant a parameter in an optimization algorithm that 
determines the acceptance criteria for a sampling jump. As will be appreciated by those skilled in the 
art, high temperature factors allow searches across a broad area of sequence space, and low 
temperature factors allow searches over a narrow region of sequence space. See Metropolis et al., J. 
Chem Phys v21. pp 1087, 1953, hereby expressly incorporated by reference. 

By "variant strand " as used herein is meant a nucleic acid strand generated using the gene assembly 
methods outlined herein to differ from the conresponding template nucleic acid sequence by at least 
one nucleotide or its complement. 

All references cited herein are expressly incorporated by reference. 
Introduction 

The present Invention is directed to methods of using computational screening of protein sequence 
libraries (that may comprise up to 10®** or more members) to select smaller libraries of protein 
sequences (that may comprise up to 10^^ members), which may then be used in a number of ways. 
For example, the proteins may actually synthesized and experimentally tested in the desired assay to 
identify proteins that possess desired properties. Sirnilarly, the library may be subjected to additional 
computational manipulation in order to create a new library, which may be experimentally tested. 

As may be appreciated by those skilled in the art, a variety of user interfaces may be utilized in the 
present invention, in a preferred embodiment, the interface is designed to maximize usability and 
efficiency. Furthermore, any or all of the computational methods described below may be automated 
for increased usability and effidency. 

Computational screening to enrich libraries with proteins possessing desired properties 

By computationally screening very large libraries of variant proteins, a greater diversity of protein 
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sequences may be screened (i.e. a larger sampling of sequence space) than is possible using 
experimental methods alone. Consequently, the probability of identifying proteins with desired 
properties is Increased and greater improvements may be realized compared to the results of purely 
experimental methods. 



The number of possible protein sequences grows exponentially with the number of positions that an 
randomized. Generally, only up to 10^^ - 10^^ sequences may be contained in a physical library 
because of experimental and physical constraints (e.g. transformation efficiency, Instrumentation 
limits, the cost of producing large numbers of biopolymers. and, for larger Obraries. the number of 
carbon atoms in the universe, etc.) Often, practical considerations may limit the library size to 10^ o 
fewer. These limits are reached for only 10 amino acid positions. In contrast, using the automated 
protein design techniques outlined below, virtual libraries of protein sequences that are vastly larger 
than experimental libraries may be generated and analyzed: up to 10®^ or more candidate sequena 
may be screened computationally. 



Using experimental methods alone, only a sparse sampling of sequences is possible in the search for 
proteins or peptides with desired properties, lowering the chance of success (both finding any proteins 
that possess the desired properties, and finding proteins that surpass the minimum acceptable 
criteria) and almost certainly missing desirable candidates. Because of the random nature of the 
mutations in experimental libraries, most of the candidates in the library are not suitable (for example, 
a large fraction of sequence space encodes unfolded, misfolded, incompletely folded, partially folded , 
or aggregated proteins), resulting in an enormous waste of the time and resources required to 
produce the library. In effect, when experimental methods alone are used, the screened library is 
composed of a large amount of "wasted sequence space". 

Computational pre-screening may be used to generate and/or enrich libraries of variant proteins that 
possess desired protein properties. An experimental library consisting of the favorable candidates 
found in the virtual library screening may then be generated, resulting in a much more efficient use of 
the time, money and effort required to construct and .screen an experimental library. In effect, when 
computational pre-screening is used the screened library Is composed of primarily productive 
sequence space. As a result, computational pre-screening increases the chances of identifying one 
or more proteins that possess the desired protein properties. 

Computational pre-screening may also be beneficial when the library of mutants is sufficiently small to 
be screened experimentally (that is. a library size of less than 10^^). It reduces the number of mutants 
that must be tested experimentally, thereby reducing the cost and difficulty associated with protein 
engineering and experimental screening. 
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While experimental methods are typicaity limited to 10^ - 10^^ sequences, computational methods 
have the unique ability to screen 10^ sequences or more. However, purely computational methods 
are limited by an incomplete knowledge of the structure-function relationship in proteins. In contrast, 
experimental methods are capable of identifying sequences with desired protein properties, even in 
cases where the causative link between sequence and obsen/ed protein properties is not understood. 
Thus, computational pre-screening followed by experimental screening of the most promising 
constructs combines the best features of computational and experinnental methods. 

Computational screening for target identification 

In a preferred embodiment, the present invention finds use in the screening of random peptide 
libraries for the purpose of target identification. In this application, random peptides are screened for 
the ability to cause a phenotypic alteration. Following identification of the active peptides, their 
interaction partners, which will typically be other proteins, may be determined. These proteins are 
lilcely to be involved in the biochemical pathway associated with a given phenotypic alteration, and 
therefore could potentially serve as new drug targets. This approach Is analogous to the chemical 
genetics methods that have been developed for small molecule libraries (Chen et al. 6:221-235 
(1999). Knodcaertetal. Chem. Biol. 7(6):411-22, (2000)). 

Screening small molecule libraries for compounds that are capable of inducing specific alterations in 
cellular physiology or phenotype has led to the discovery of proteins that function in a variety of 
biochemical and signal transduction pathways. Cyclosporin A (CsA) and FK506, for example, were 
selected in standard pharmaceutical screens for inhibition of T-cell activation. It Is noteworthy that 
while these two drugs bind completely different cellular proteins, cyclophilin and FK506 binding 
protein (FKBP), respectively, the effect of either drug is virtually the same: profound and specific 
suppression of T-cell activation, phenotypically observable In T cells as inhibition of mRNA production 
dependent on transcription factors such as NF-AT and NF-KB. 

Chemical genetics approaches have typically used libraries of small molecules; however, libraries of 
peptides or proteins could be used instead. Computational pre-screening of the peptide libraries 
could be used to maximize the diversity of properties In the library and to select structured peptides 
that are especially likely to bind other molecules with high affinity. 

Computational screening for fold identification 

The present invention also finds use In fold identification. Structural and functional properties of 
protein sequences, such as those deriving from various genome projects, may often be inferred from 
sequence similarity to proteins whose structural and/or functional properties have been characterized. 
One limitation of this approach is that many newly discovered sequences lack sufficient sequence 
similarity witti any of the better characterized proteins. 
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In a preferred embodiment, a three-dimensional database Is created by modifying a known protein 
structure to incorporate particular amino acid residues required for a characteristic property or 
function, as is described in WO 00/23474. expressly incorporated herein by reference. This allows 
the creation of a database that can be used in a manner similar to other "structural alignment" 
programs. That is, by using the protein design cyde systems outlined herein, a variety of amino acid 
sequences that will take on a particular structural fold are generated. These sequences represent a 
set of artificial sequences that will take on a particular conformation. This database may be searched 
against protein databases to identify new proteins having structural similarity witti the known protein. 
Thus, proteins can be identified that make take on a particular fold but do not have enough sequence 
homology to a naturally occurring protein to be chosen using known alignment programs. In some 
cases, this will allow the assignment of putative functional Information as well; for example, by 
identifying proteins witii sti-uctural homology to a particular class of enzyme or ligand, the new protein 
can be assigned similar function. This finds particular use in Identifying proteins that have been 
sequenced but to which no structure and/or function has been assigned. 

In addition, the database could contain additional computatfonally generated sequences that are 
predicted to be compatible witii a given structure and/br function. Computatfonally supplemented 
databases may contain a signifrcantiy greater diversity and total number of sequences tiian databases 
ttiat rely solely on experimental results. Consequently, the fraction of sequences that may be 
classified into a protein family will be larger using a computationally supplemented database than 
using a purely experimental database. Fold identification using PDA™ technotogy and biolnformatics 
tools (e.g. dynamic programming algorittims, BLAST search), may tfien be used to identify new dmg 
targets and antidotes to biological weapons. 

As an example of this concept, the sequencing of new genomes will reveal proteins, stiiictural motifs, 
and domains tiiat are unique to certain genomes. For example, there may be some domains that are 
unique to bacterial or viral genomes and do not exist In eukaryotic genomes. PDA™ technology 
and/or the ottier computational mettiods outiined herein may be used to identify sequences tiiat are 
compatible witii these structures. Bacterial and viral genomes may then be searched to identify 
additional proteins that are likely to fold to tfie structures, but could not be identified as homologs 
using traditional methods. The resulting proteins may serve as novel drug targets that could be used 
to discover new classes of antibiotics and antiviral drugs. 

Approach to library generation 

The invention describes novel mettiods to create secondary libraries derived from very large 
computational mutant libraries. These metiiods altow the rapkJ experimental and/or computational 
testing of large numbers of computationally designed sequences. 
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As more fully outlined below, the invention may take on a wide variety of configurations. In general, 
primary libraries are generated computationally. This may be done in a wide variety of ways, 
including, but not limited to, sequence alignments of related proteins, structural alignments, structural 
prediction models, SCMF methods, or preferably protein design automation™ (PDA™) teclinology 
computational analysis. 

Once the primary library-is generated, it may be manipulated in a variety of ways. In one 
embodiment, a different type of computational analysis may be done; for example, a new type of 
ranking may be performed. In a preferred embodiment, some subset of the primary library is then 
experimentally generated to form a secondary library. Alternatively, some or all of the primary library 
members are recombined to form a secondary library, resulting in a secondary library that contains 
sequences not included in the primary library. Again, this may be done either computationally or 
experimentally or both. 

Accordingly, the present invention provides computational and experin^ntal metiiods for generating 
secondary libraries of scaffold protein variants. 

Oven^iew of PDA™ Technology Methodotogy 

In a preferred embodiment, tiie computational method used to generate the primary library is Protein 
Design Automation™ (PDA™) technology, as is described in U.S.S.N.s 60/061,097, 60/043,464, 
60/054,678, 09/127,926 and PCT US98/07254. ail of which are expressly incorporated herein by 
reference. Briefly, PDA™ technology may be described as follows. A known, generated or 
homologous protein structure is used as the starting point. The residues to be optimized are then 
identified, which may be the entire sequence or subset(s) tiiereof. The side chains of any positions to 
be modified are then removed. The amino acids that will be considered at each position are selected, 
(for example, core residues generally will be selected from the set of hydrophobic residues, surface 
residues generally will be selected from the hydrophilic residues, and boundary residues may be 
eitiier). Each amino add residue may be represented by a discrete set of allowed conformations, 
called rotamers. Interaction energies are calculated between each rescue in a given rotamer and the 
backbone and between each pair of residues in each of their rotamers at different positions. 
Combinaforial search algorithms, typically DEE and Monte Cario, are used to identify the optimum 
amino acid sequence and additbnal low energy sequences which will comprise tiie primary library. 

PDA™ technology, viewed broadly, has four components that may be varied to alter the output (i.e. 
the primary library): generation of the template or templates, choice of amino acid identities and 
conformations considered at each position, the scoring functions used in the process; and the 
optimization strategy. 
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Selection and preparation of the scaffold protein 
Source of Three-dimensfonal Structure 

The scaffold protein may be any protein for which a three dimensional structure {that is, three 
dimensional coordinates for each atom of the protein) is known or may be generated. The three 
dimensional structures of proteins may be detennined using X-ray crystallographic techniques. NMR 
techniques, de nwo modeling, homology modeling, etc. In general, if X-ray structures are used, 
structures at 2 A resolution or better are preferred, but not required. Suitable protein structures 
include, but are not limited to, aD of those found In the Protein Data Base compfled and serviced by 
the Research CoHaboratoiy for Stnictura! Bioinformatics (RCSB. formerly the Brookhaven Natfonal 
Lab). 

Scope of Scaffold 

TTie scaffoW used in protein design calculations may comprise an entire protein or peptide, a subset 
of a protein such as a domain (including functional domains such as enzymatic domains, substrate- 
binding domains, regulatory domains, dimerization domains, etc.). motif, site, or loop. The scaffoki 
protein may comprise more than one protein chain. That is, the scaffold may be an oligomer 
(including but not limited to dimers, trimers, hexamers, 60-mers such as viral coats, and long protein 
chains such as actin filaments) or a multi-protein complex (including but not limited to ligand-receptor 
pairs, antibody-antigen pairs, ribosome complexes, proteosome complexes, transcription complexes, 
chaperone complexes, the spHcesome. molecular motors, focal adhesbn complexes, multi-protein 
signaling complexes, etc.). TTie scaffold may additnnally contain non-protein components, including 
but not limited to smaO molecules, substrates, cofactors, metals, water molecules, prosthetic groups, 
nuclefc acids such as DNA and RNA. sugars, and lipids. 

Source of the scaffold protein 

The scaffold proteins may be from any organism, including prokaryotes and eukaryotes, with proteins 
from bacteria, fungi, vimses, extremophiles such as the archaebacteria. insects, fish, animals 
(particularly mammals and particularly human) and birds all possible. The scaffold protein does not 
necessarily need to be naturally occurring, for example the scaffoki protein couW be a designed 
protein, or a protein selected by a variety of methods including but not limited to directed evolution 
(Farinas at al. Current Opinion in Biotechnotogy 12:545-551 (2001) lUlorawski et al. Biotechnotogy and 
Bioengineering 76:99-107 (2001), Stammer Nature 370(6488); 389-91 (1994) Ness et al. Adv 
Protein. Chem. 55:261-92 (2000)). DNA shuffling (Maxygen, Enchira. Diversa) or ribosome display 
(Hanes et al. Methods in Enzymology 328:404-430 (2000); Hanes and Pluckthun. Proc. Natl. Acad. 
Sci. USA 94:4937-4942 (1997); Roberts and Szostak. Proc. Natl. Acad. Sd USA 94 12297-302 
(1997). 
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Examples of suitable scaffolds 

As will be appreciated by those skilled in the art. any number of scaffold proteins find use in the 
present invention. Suitable proteins include, but are not limited to, industrial and pharmaceutical 
proteins, including ligands, ceil surface receptors, antigens, antibodies, cytokines, hormonps, 
transcription factors, signaling modules, cytoskeletai proteins and enzymes. 

Specifically, preferred scaffoki proteins include, but are not limited to, those with known or predictable 
structures (including variants): 

• cytokines (IHra (+receptor complex), IL-1 (receptor alone), IL-1a, IL-1b (Including variants 
and or receptor complex). 11-2. IL-3. IL-4, IL-5. IL-6, IL-8, IHO. IFN-f , INF-y, IFN-a-2a; IFN- 
a-2B, TNF-a; CD40 ligand (chk). Human Obesity Protein Leptin. Granulocyte Colony- 
Stimulating Factor, Bone Morphogenetic Protein-7, Ciliary Neurotrophic Factor, Granulocyte- 
Macrophage Colony-Stimulating Factor, Monocyte Chemoattractant Protein 1. Macrophage 
Migratfon Inhibitory Factor, Human Glycosylation- Inhibiting Factor, Human Rantes, Human 
Macrophage Inflammatory Protein 1 Beta, human growth hormone, Leukemia Inhibitory 
Factor, Human Melanoma Growth Stimulatory Activity, neutro|:^il activating peptide-2, Cc- 
Chemokine Mcp-3, Platelet Factor M2, Neutrophil Activating Peptide 2. Eotaxin, Stromal Cell - 
Derived Factor-1, Insulin, Insulin-like Growth Factor I, Insulin-like Growth Factor II, 
Transforming Growth Factor B1, Transforming Growth Factor B2, Transforming Growth 
Factor B3, Transforming Growth Factor A, Vascular Endothelial growth factor (VEGF), acidic 
Fibroblast growth factor, basic Fibroblast growth factor, Endothelial growth factor, Nerve 
growth factor, Brain Derived Neurotrophic Factor, Ciliary Neurotrophic Factor, Platelet 
Derived Growth Factor, Human Hepatocyte Growth Factor. Fibroblast Growth Factor 
(including but not limited to alternative splice variants , abundant variants, and the like). Glial 
Cell-Derived Neurotrophic Factor, and haemopoietic receptor cytokines (including but not 
limited to erythropoietin, thrombopoietin, and prolactin), APM1 (including, but not limited to 
adipose most abundant gene transript 1 ), and the like. 

• other extracellular signaling moieties, including, but not limited to, Sonic hedgehog, protein 
hormones such as chorionic gonadotrophin and leutenizing hormone. 

• blood clotting and coagulation factors including, but not limited to, TPA and Factor Vila; 
coagulation factor IX; coagulation factor X ; PROTEIN S protein; Fibrinogen and Thrombin; 
ANTITHROMBIN III; streptokinase and unDkinase. retevase, and the like. 

• transcription factors and other DNA binding proteins, including but not limited to, histones, 
p53; myc; PIT1; NFkB;AP1;JUN; KD domain, homeodomain, heat shock transcription factors, 
Stat, zinc finger proteins (e.g. zif268). 

• Antibodies, antigens, and trojan horse antigens, Including, but not limited to, immunoglobulin 
super family proteins, including but not limited to CD4 and CDS, Fc receptors. T-cell 
receptors. MHC-I, MHC-ll, CD3, and the like. Also, immunoglobulin-like proteins, including 
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but not limited to fibronectin. pkd domain, integrin domains, cadhrin, invasins, cell surface 

receptors witii Ig-Iike domains, and the Bke. Intrabodles. and the like; Anti-HerC neu antibody 

(e.g. Hercepbn); Anti-VEGF; Anti-CD20 (RItuxan), among others. 

intracellular signaling modules, including, but not limited to, kinases, phosphatases. G- 

proteins Phosphatidylinositol 3-kinase (PI3-kinase) kinase. Phosphatidylinositol 4-kinase. wnt 

family members irnduding but not limited to wnt-1 through wnt 15, EF hand proteins including 

calmodulin, troponin C, S100B, calbindtn and D9k; NOTCH; MEK; MARK; ubitquitin and 

ubiquitin Kke proteins, including UBL1, UBL5. UBL3 and UBL4, and the like. 

viral proteins, including, but not limited to, hemagglutinin trimerizaHon domain and HIV Gp41 

ectodomain (fusion domain); viral coat proteins, viral receptors, integrases. proteases. 

reverse transcriptases. 

receptors, including, but not limited to, the extracellular region of human tissue factor 
cytokine^inding region Of Gp1 30, G-CSF receptor, erythropoietin receptor, Rbroblast Growth 
Factor receptor. TNF receptor, IL-1 receptor. IL-1 receptor/IL1ra complex, \L-4 receptor, INF-y 
receptor alpha chain. MHC Class I. MHC Class II . T Cell Receptor, Insulin receptor, insulin 
receptor tyrosine kinase and human growth hormone receptor; Lectins; GPCRs, including but 
not limited to G-Protein coupled receptors; ABC Transporters/ Multidrug resistance proteins; 
Na and K channels; Nuclear Hormone Receptors; Aquaporins; Transporters. RAGE (receptor 
for advanced glycan end points). TRK -A, -B, -C, and the like, and haemopoiedc receptors, 
enzymes including, but not limited to. hydrolases such as proteases/proteinases, 
synthases/synthetasesfligases, decarboxylases/lyases, peroxidases, ATPases' 
carbohydrases. lipases; isonrerases such as racemases. epimerases. tautomerases. or 
mutases; transferases, hydrolases, kinases, reductases/oxidoreductases. hydrogenases. 
polymerases, phophatases, and proteasomes anti-proteasomes. (e.g., MLN341). Suitable 
enzymes include, but limited to. those listed In the Swiss-Prot enzyme database. 
Additional proteins including but not limited to heat shock proteins, ribosomal proteins, 
glycoproteins, motor proteins, transporters, drug resistance proteins, kinetoplasts and 
chaperonins. 
/^timkrobial peptkles 

small proteins including but not limited to metal ligand and disulfide-bridged proteins such as 
metallothbnein. Kunitiz-type inhibBors. crambin. snake and scorpton toxins, and trefoil 
proteins; antimicrobial peptides such as defensins. thoredoixn. fereodoxin, transferetln. and 
the like. 

protein domains and motifs including, but not limited to, SH-2 domains. SH-3 domains. 
Pleckstrin homology domains, WW domains. SAM domains, kinase domains, death domains." 
RING finger domains, Kringle domains, heparin-binding domains, cysteine-rich domains, 
leudne zipper domains, zinc finger domains, nucleotide binding motifs, transmembrane 
heltoes. and helix-turn-helw motifs. Additionally. ATP/GTP-binding site motif A; Ankyrin 
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repeats; fibronectin domain; Frizzled (Iz) domain; GTPase binding domain; C-type lectin 
domain; PDZ domain; 'Homeobox* domain; KrQeppel-associated box (KRAB); Leucine zipper; 
DEAD and DEAH box families; ATP-dependent lielicases; HMG1/2 signature; DNA mismatch 
repair proteins mutL / hexB / PMS1 signature; Thioredoxin family active site; Thioredoxins; 
Annexins repeated domain signature; Clathrin liglit chains signatures; Myotoxins signature; 
Staphylococcal enterotoxins / Streptococcal pyrogenic exotoxins signatures; Serpins 
signature; Cysteine proteases inhibitors signature; Chaperonins; Heat shock; WD domains; 
EGF-like domains; Immunoglobulin domains, Immunoglobulin-like proteins and the like, 
specific protein sites or other subsets of resklues, including but not limited to protease 
cleavage/ recognitksn sites, phosphorylation sites, metal binding sites, and signal sequences. 
Additionally, proteins having post-translational nruxiiflcations include, but are not limited to: 
glycosylation site; Oglycosylation site; Glycosaminoglycan attachment site; Tyrosine suKatbn 
site; cAI\/IP- and cGMP; dependent protein kinase phosphorylation site; Protein kinase C 
phosphorylation site; Casein kinase I! phosphorylation site; Tyrosine kinase phosphorylation 
site; N-myrlstoylatlon site; Amidation site; Aspartic acid and asparagine hydroxylatbn site; 
Vitamin K-dependent carboxylation domain; Phosphopantetheine attachment site; Prokaryotfc 
membrane lipoprotein lipki attachment site; Prokaryotic N-temiinal methylatbn site; Prenyl 
group binding site {CMX box); Intein N- and C-terminal spFicing motif profiles, and the like. 
Proteins involved in motility, including but not limited to chemokines, 8100 family proteins 
(Including but not limited to NRAGE). 
Peptides - defensins 

peptide ligands including, but not limited to, a short region from the HIV-1 envelope 
cytoplasmic domain (shown to block the action of cellular calmodulin), regions of the Fas 
cytoplasmic domain (death-inducing apoptotic or G protein inducing functions), magainin, a 
natural peptide derived from Xenopus (anti-tumor and anti-microbial activity), short peptide 
fragments of a protein kinase C isozyme, &PKC (blocks nuclear translocation of full-length 
UPKC in Xenopus oocytes follov\^ing stimulation), SH-3 target peptides, naturitic peptides 
(AMP, BMP, and CMP), and fibrinopeptides and neuropeptides. 

presentation scaffolds or "ministructures" including, but are not limited to, minibody structures 
(see for example Blanch! et al., J. Mol. Biol. 236(2):649-59 (1994), and references dted 
therein, all of which are incorporated by reference), maquettes (Grosset et al. Biochemistry 
40:5474-5487 (2001)), k)ops on beta-sheet turns and coiled-coil stem structures (see, for 
example, Myszka et al., Bbchem. 33:2362-2373 (1994) and Martin et al., EMBO J. 
13(22):5303-5309 (1994), incorporated by reference), zinc-finger domains, transglutaminase 
linked structures, cyclic peptides. B-loop structures, coiled coils, helical bundles, helical 
hairpins, and beta hairpins. 

Ion channel protein domains, including but not limited to sodium, calcium, potassium, and 
chloride, including their component subunit. Examples of extracellular ligand-gated ion 
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Channels include nAChR receptors. GABA and glycine. 5H-T. MOD-I. P(2X). glutamate. 
NMDA. AMPA, Kainate receptors. GluR-B. ORCC. P2X3. Inward rectifying channels. ROMK.' 
IRK. BIR, and the like. Examples of voltage-gated ion channels. Examples of intracellular 
ligand-gated ion channels. Mechanosensative and cell volume-regulated ion channels and 
the like. 

In addition, a preiened embodiment utilizes scaffold proteins such as random peptWes. That is. there 
is a significant amount of work being done in the area of utilizing random peptides In high throughput 
screening technkiues to Wentify biologically relevant (particularly disease states) proteins. The 
methods of the invention are particularly relevant for computatfonally prescreening random peptide 
libraries to drastically reduce the amount of wet chemistry that must be done, by removing sequences 
that are unlikely to be successful. Different design criteria can be used to produce candidate sets that 
are biased for properties such as charge, solubility, or active site characteristics (polarity, size), are 
biased to have certain amino acids at certain positions or to lake on certain folds. Tliat is. the 
peptkJes (which may be the scaffold protein or the candWate agents, as outlined betow) are 
randomized, either fully randomized or they are biased in their randomizatton. e.g. in 
nucleotide/rBsidue frequency generally or per position. By "randomized- or grammatical equivalents 
herein is meant that each nucleic acid and peptide consists of essentially random nucleotides and 
amino acids, respectively. Thus, any amino acid residue may be incorporated at any position. The 
synthetic process can be designed to generate randomized peptides and/or nucleic acids, to allow the 
fomiatfon of all or most of the possible combinatfons over the length of the nucleic acid, thus forming 
a library of randomized candidate nucleic acids. 

In one embodiment, the library is fully randomized, with no sequence preferences or constants at any 
position. In a preferred embodiment, the library is biased. That is, some positions within the 
sequence are either heki constant, or are selected from a limited number of possibilities. For 
example, in a prefened embodiment, the nucleotides or amino ackJ residues are randomized within a 
defined class, for example, of hydrophobic amino acids, hydrophilic resWues. sterically biased (either 
small or large) residues, towards the creation of cysteines, for cross-linking, prolines for SH-3 
itomains. serines, threonines, tyrosines or histkJines for phosphoiylatton sites, etc.. or to purines, ete. 
In a preferred embodiment, the bias is towards peptides or nuclefc acids that interact with known 
classes of molecules. For example, it is known that much of intracellular signaling is carried out via 
short regions of polypeptides interacting with other polypeptides through small peptide domains. In 
addition, agonists and antagonists of any number of molecules may be used as the basis of biased 
randomizatton of candidate bioactive agente as well. 
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In general, the generation of a prescreened random peptide libraries may be described as follows. 
Any structure, whether a known structure, for example a portion of a known protein, a known peptkie, 
etc., or a synthetic structure, can be used as the backbone for computational screening. For 
example, structures from X-ray crystallographic techniques, NMR techniques, de novo modelling, 
homotogy modelling, etc. may all be used to pick a backbone for which sequences are desired. 
Similarly, a number of molecules or protein domains are suitable as starting points for the generation 
of biased randomized candidate bloactive agents. A large number of small molecule domains are 
known, that confer a common function, structure or affinity. In addition, as is appreciated in the art, 
areas of weak amino ackJ homology may have strong structural homology. A number of these 
molecules, domains, and/or corresponding consensus sequences, are known, including, but are not 
limited to, SH-2 domains^ SH-3 domains, Pleckstrin, death domains, pn^tease cleavage/recognltton 
sites, enzyme inhibitors, enzyme substrates, Traf. etc. Similarly, there are a number of known nucleic 
acid binding proteins containing domains suitable for use in the invention. For example, leucine 
zipper consensus sequences are known. Thus, in general, known peptide ligands can be used as the 
starting scaffold backbone for the generation of the primary library. 

in a prefenred embodiment, the scaffold protein is a variant protein, including, but not limited to, 
mutant proteins comprising one or a plurality of substitutions, inseriiions or deletions, including 
chimeric genes, and genes that have been optimized in any number of ways, including experimentally 
or computationally. 

In a preferred embodiment, the scaffold protein is a. chimeric protein. A chimeric protein (sometimes 
refenred to as a "fusion protein") in this context means a protein that has sequences from at least two 
different sequences operably linked or fused. The chimeric protein may be made using either a single 
linkage point or a plurality of linkage points. In addition, the source of the parent protein sequences 
may be as listed above for scaffokl proteins, e.g. prokaryotes, eukaryotes, including archebacteria 
and viruses, etc. 

As will be appreciated by those in the art, chimeric proteins may be made from different naturally 
occun^ing proteins in a gene family (e.g. one with recognizable sequence or structural homology) or by 
artificially joining two or more distinct genes. For example, the binding domain of a human protein 
may be fused with the activation domain of a mouse gene, etc 

The sequence of the chimeric gene may be been constructed synthetically (e.g. arbitrary or targeted 
portions of two or more genes are crossed over randomly or purposely), experimentally (e.g. through 
homologous recombination or shuffling techniques) or computationally (e.g. using genetic annealing 
programs, "in silico shuffling", alignment programs, etc.). For the purposes of the invention, these 
techniques can be done at the protein or nucleic acid level. 
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In a preferred embodiment, the scafTold protein is actually a product of a computational design cycle 
and/or screening process. That is. a first round of the methods of the invention may produce one or 
more sequences for which further analysis is desired. 

Although several classes of proteins have been stated herein, this should not be construed as an 
exhaustive list, but rather some examples of proteins that may be optimized using the computational 
methodologies outlined herein, including PDA™ technology. 

Preparation of Protein Backbone for Calculations 

The protein scaffold may be modified or altered at the beginning (and optionally, but not preferably in 
the middle or end) of a protein design calculation, or the unaltered scaffold may be used It is also' 
possible to use methods in which the protein scaffold is modified during later steps of a design 
calculation, including during the energy calculation and optimization steps. 

In a preferred embodiment, protein scaffold backbone (comprising, the nitrogen, the carbonyl carbon 
the a-carbon. and the carbonyl oxygen, along with the direction of the vector from the a-carbon to the 
PHarbon) may be altered prior to the computational analysis, for example by varying a set of 
parameters called supersecondary structure parameters. See for example U S Patent Nos 
6.269,312. 6.188.965. and 6.403.312. all of which are herein expressly incorporated by reference 
Altemafively. the protein scaffoH is altered using other methods, such as manually, indu ding directed 
or random perturbations 

Most protein structures contain loop regions that are flexible or conformationally heterogeneous The 
protein backbone may be modified In the toop regions using methods such as molecular dynamics 
simulattons and analysis of databases of known loop structures. In addition, loops may be modiffed in 
order to incorporate new structural or functional properties such as new binding sites. 

in a preferred embodiment, the design cyde is done using a piuralify or set of scaffold proteins That 

the scaffold may be a set of protein structures created by perturbing the starting stnicture This 
may be done using any number of techniques, including molecular dynamics and IMonteCarto 
analysis, that alter the protein structure (induding changing the backbone and side chain torsion 
angles.) Alternatively, an ensemble of structures such as those obtained from NMR may be used as 
the scaffold. These badcbone modifications are particularly useful for enhancing the diversity of 
sequences derived from protein design simulations. Similarly, other useful ensembles indude sets of 
related proteins, sets of related structures, artifidal created ensembtes, etc. 

In a preferred embodiment, once a protein structure badtbone is generated (with alterations, as 
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outlined above), explicit hydrogens are added if not included within the structure . For example, if the 
structure was determined using X-ray crystallography, hydrogens are typically added. 

In a preferred embodiment, energy minimization of the structure is run to relax strain, including strain 
due to van der Waals clashes, unfavorable bond angles, and unfavorable bond lengths. In an 
especially prefenred embodiment, this is done by doing a number of steps of conjugate gradient 
minimization (see Mayo et al., J. Phys. Chem. 94:8897 (1990)) of atomic coordinate positions to 
minim'ee the Dreiding force field with no electrostatics. Generally from 1 0 to 250 steps is preferred, 
with 50 steps being most preferred. 

Identification of Variable, Floated, and Fixed Positions 

In a preferred embodiment all of the residue positions of the protein are variable. This is particularly 
desirable for smaller proteins, although the present methods allow the design of larger proteins as 
well. In an alternate prefen-ed embodiment, only some of the residue positions of the protein are 
variable, and the remainder are fixed or floated. In this embodiment, the variable residues may be at 
least one, or anywhere firom 0.001% to 99.999% of the total number of residues. Thus, for example, It 
may be possible to change only a few (or one) residues, or most of the residues, with all possibilities 
in between. 

In an alternate embodiment, only one or two residue positions are variable and the residue positions 
within a small distance of, for example, 4A to 6A of the variable residue positions are floated. In this 
embodiment, it is possible to conduct separate calculations for different positions and then combine 
the results to yield protein variants with multiple mutations. Using the results from one calculation as 
a starting point for the next calculation one residue position at a time, the optimization procedure may 
be iterative, iteration may be performed until a consistent result is reached. 

In a prefenred embodiment, residues which may be fixed include, but are not limited to, structurally or 
biologically functional residues. For example, residues which are known to be important for biological 
activity, such as the residues which form the active site of an enzyme, the substrate binding site of an 
enzyme, the binding site for a binding partner (llgand/receptor, antigen/antibody, etc.), 
phosphorylation or glycosylation sites, or structurally Important residues, such as cysteines 
participating in disulfide bridges, metal binding sites, critical hydrogen bonding residues, residues 
critical for backbone conformation such as proline or glycine, residues critical for packing interactions, 
etc. may all be fixed or floated. 



31 



wo 03/014325 



PCT/US02/25588 



Similarly, residues which may be chosen as variable residues may be those that confer undesirable 
biological attributes, such as susceptibility to proteolytic degradation, unwanted oligomerization or 
aggregation, glycosylation sites which may lead to unwanted immune responses, unwanted binding 
activity, unwanted allostery, undesirable enzyme activity, etc. 

Alternatively, residues that confer desired protein properties may be specifically targeted for variation. 
In a preferred embodiment, this design strategy may be used to alter properties such as binding 
affinity and specificity and catalytic efficiency and mechanism. A region such as a binding site or 
active site may be defined, for example, to include all residues witiiin a certain distance, for example 4 
- 1 0 A, or preferably 5 A. of the residues that are in van der Waals contact with the substrate or 
ligand. Alternatively, a region such as a binding site or active site may be defined using experimental 
results, for example, a binding site could include all positions at which mutation has been shown to 
affect binding. 

Select Amino Acids to be Considered at Each Position 

A set of amino add side chains is assigned to each variable position. That is, the set of possible 
amino acid side chains tiiat will be considered at each particular position is chosen. In one 
embodiment, variable positions are not classified and all amino acids are considered at each variable 
position. Alternatively, a subset of amino acids are considered at each variable position. Mettiods for 
determining subsets of amino acids include, but are not limited to, those discussed below. Any 
combination of classification methods, including no classification, may be applied to ttie different 
variable positions. 

In a preferred embodiment, all amino acid residues are allowed at each variable residue position 
identified in ihe primary library. That is. once the variable residue positions are identified, a 
secondary library comprising every combination of every amino acid at each variable residue position 
is made. 
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In a preferred embodiment, subsets of amino acids are chosen to maximize coverage. Additional 
amino acids with properties similar to those contained within the primary library may be manually 
added. For example, if the primary library includes three large hydrophobic residues at a given 
position, the user may chose to include additional large hydrophobic residues at that position when 
generating the secondary library. In addition, amino acids in the primary library that do not share 
similar properties with most of the amino acids at a given posifion may be excluded from the 
secondary library. Alternatively, subsets of amino acids may be chosen firom the primary library such 
that a maximal diversity of side chain properties Is sampled at each position. For example, if the 
primary library includes three large hydrophobic residues at a given position, the user may chose to 
include only one of them in the secondary library, in combination with other amino acids that are not 
large and hydrophobic. 

In a prefenred embodiment, each variable posifion Is classified as either a core, surface or boundary 
residue position. The classification of residue positions as core, surface or boundary may be done in 
several ways, as will be appreciated by those in the art. In a prefenred embodiment, the classification 
is done via a visual scan of the original protein scaffold and assigning a classification based on a 
subjective evaluation of one skilled In the art of protein modeling. Alternatively, a preferred 
embodiment, called RESCLASS, utilizes an assessment of the orientation of the Ca-Cp vectors 
relative to a solvent accessible surface computed using only the template Co atoms, as outlined in 
U.S. Patent Nos. 6,269.312, 6.188,965, and 6,403,312, and expressly herein Incorporated by 
reference. Alternatively, a surface area calculation may be done. In an alternate embodiment, the 
results of the RESCLASS calculation are used in conjunction with the results of a surface area 
calculation in order to classify residue positions. 

A core residue will generally be selected ftt>m a set of hydrophobic residues consisting of alanine, 
valine, isoleucine, leucine, phenylalanine, tyrosine, tryptophan, and methionine (in some 
embodiments, methionine may be removed from the set). Similarly, surface positions are generally 
selected from a set of hydrophilic residues consisting of alanine, serine, threonine, aspartic acid, 
asparagine, glutamine. glutamic acid, arginine, lysine and histidine. Finally, boundary positions are 
generally chosen from alanine, serine, threonine, aspartic acid, asparagine. glutamine. glutamic acid, 
arginine, lysine histidine, valine, isoleucine, leucine, phenylalanine, tyrosin e, tryptophan, and 
methionine. 
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In a preferred embodinrjent, proline, cysteine and glycine are not included In the list of possit>le amino 
acid side cliains. and thus the rotamers for these side chains are not used. However, in an alternate 
preferred embodiment, when the variable residue position has a <p angle (that is, the dihedral angle 
defined by 1) the carbonyl carbon of the preceding amino acid; 2) the nitrogen atom of the current 
residue; 3) the a-carbon of the current residue; and 4) the carbonyl carbon of the current residue) 
greater than 0 degrees, the position is set to glycine to minimize backbone strain. In an alternate 
enr4)odiment. q^steine is considered at positions where disulfide bonds are desired. In another 
aKemate embodiment, proline is considered at positions whose backbone conlbrmatton is altowabte 
for proline. 



As wfll be appreciated by those in the art, there is a computational benefit to classifying the residue 
positions, as it decreases the combinatorial complexity of the problem. It should also be noted that 
there may fee situations where alternative classification approaches will be applied or where the sets 
of core, boundary and surface residues are altered from those described above; for example, under 
some circumstances, one or more amino ackJs is either added or subtracted from the set of allowed 
amino acids. For example, hydrophobfc resMues may be included at solvent exposed posiUons in 
order to confer desired oligomerization or ligand binding activity, and polar residues may be Included 
in the core of the protein in order to construct an active site or a binding site. Similarly, i n one 
embodiment, only residues capable of fomiing N-capping interactions are included at the position 
immediately preceding each helix, and amino adds that interact unfavorably with the helix dipole are 
subtracted from the set of polar residues at the three positions at the beginning and end of each helix. 

In a preferred embodiment, the set of amino adds allowed at each positton is detenrtned using 
sequence or structure alignment methods. For example, the set of amino acMs allowed at each 
position may comprise the set of amino acids that is observed at that positron in the alignment, or the 
set of amino ackJs that is observed most frequently in the alignment 

In another preferred embodiment, the set of amino acids allowed at each positron comprises the set of 
amino acWs that are known to interact with a partfcular dass of molecules or to serve a specific 
iunction. Possible sete indude, but are not limited to, residues that may Ggate or coordinate to certain 
metels (such as zinc, copper, iron, and nrolybdenum), residues that may undergo posttranslational 
modification (such as phosphorylation, glycosylaton, prenylatfon. and Bpklation), and residues that 
are amenable to synthetic modification. Synthetic modifications indude, but are not limited to, 
alkylation or acylation which indudes but is not limited to PEGylation, biotinylation. fluorophore 
conjugation, acetylation, oxidative or reductive homo- or heterooligomerization, native ligation, 
conjugation to synthefic mono- and oligosaccharides, and covalent or non-covalent attachment to a 
solid support (e.g. glass beads, glass slides, or 96-well plates). Sites of synthetic modiffcations 
indude. but are not Bmited to. the amide N-H. the amino add side chains, the amino or carboxyl 
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terminus of the protein, or any of the various posttranslational modifications. 

In a preferred embodiment, the set of allowed amino adds includes one or more non-natural or 
noncanonical amino acids. Synthetic modifications of the non-natural or non-canonical amino acids 
are also viable. In addition to the modifications listed above, these synthetic transformations include, 
but are not limited to intra- and intermolecular metal mediated couplings such as the Heck reaction or 
Suzuki coupling and conjugation through shiff base formation. In a prefenred embodiment, the set of 
allowed amino acids includes more than one charge state for some or all of the acidic or basic 
residues (that is, arginine, lysine, histidine. glutamic ackJ, aspartic acid, cysteine, and tyrosine). 

Select the Set of Rotamers That Will Be Used to Model Each Residue Type 
In a preferred embodiment, a set of discrete side chain confomnations, called rotamers, are 
considered for each amino acid. Thus, a set of rotamers will be considered at each variable and 
floated position. Rotamers may be obtained from published rotamer libraries (see Lovel et al., 
Proteins: Structure Function and Genetics 40:389-408 (2000) Dunbrack and Cohen Protein Science 
6:1661-1681 (1997); DeMaeyer et al., Folding and Design 2:53-66 (1997); Tuffery et al. J. Biomol. 
Struct. Dyn. 8:1267-1289 (1991), Ponder and Richards, J. Mol. Biol. 193:775-791 (1987)), from 
molecular mechanics or ab initio calculations, and using other methods. In a preferred embodiment, a 
flexible rotamer model is used (see Mendes et al., Proteins: Structure, Function, and Genetics 
37:530-543 (1999)) Similarly, artificially generated rotamers may be used, or augment the set chosen 
for each amino acid and/or variable position. In a prefenred embodiment, at least one conformation 
that is not low in energy is included in the list of rotamers. In an alternative embodiment, the identity of 
each amino acid, rather than specific conformational states of each amino ackl, are used, i.e., use of 
rotamers is not essential. 

Generating ranks or lists of possible sequences 

In essence, any computational methods that may resuR in either the relative ranking of the possible 
sequences of a protein or a list of suitable sequences may be used to generate a primary library. As 
will be appreciated by those in the art, any of the methods described herein or known in the art may 
be used. Each method may be used alone, or in combination with other methods. In a preferred 
embodiment, knowledge-based and statistical methods are used. Altematively, methods that rely on 
energy calculations may also be used. Protein design methods use various criteria to screen 
sequences, resulting in sequences that are likely to possess desired properties. The design criteria 
may be altered to generate primary libraries that are likely to contain proteins possessing a different 
set of desired properties. 

Knowledge-based and Statistical Methods 

In a preferred embodiment, sequence and/or stmctural alignment programs may be used to generate 
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primary Dbraries. For example, various algnmerrt methods may be used to create sequence 
alignments of proteins related to the target structure (see for erampte Altschul et al.. J. Mol. Biol. 
215(3): 403 (1990). incorporated by reference). Sequences may be related at the level of primary, 
secondary, or terUary structure. Alternatively, sequences may be related by function or activlfy. 
These sequence alignments are then examined to determine the observed sequence variations. 
These sequence variations are tabulated to define a primary library, or used to bias the convergence 
of a protein design algorithm. 



As is known in the art. sequence aignmente can be analyzed usaig statistical methods to calculate 
the sequence diversity at any position in the alignment, and the occurrence frequency or probability of 
each amino acid at a position. In the simplest embodiment, these occurrence frequencies are 
calculated by counting the number of times an amino acid is observed at an alignment position, then 
dividing by the total number of sequences in the alignment. In other embodiments, the contribution of 
each sequence, position or amino acid to the counting procedure is weighted by a variety of possible 
mechanisms. For example, sequences may be weighted towards or away from a wild type sequence, 
towards a human sequence, etc. 



Furthermore, the sequence alignments may be analyzed to produce the probability of observing two 
residues simultaneously at two positions. These probabilities may serve as a measure of the strength 
of coupling between residues. In one embodiment, the probabilities may then be used to favor 
selection of sequences ttiat maintain conserved residue pairs and disfavor selection of sequences 
that contain pairs that are seldom or never observed In sequence homologs. 

As is known in the art there are a number of sequence-based alignment programs; including for 
example, SmittvWaterman searches. Needleman-WUnsch. Double Afiine Smith-W&terman, frame 
search, Gribskov/GCG profile search, Gribskov/GCG profile scan, profile frame search, Bucher 
generalized profiles. Hidden Markov models, Hfiame. Double Frame, Blast, Psi-Blast, Clustal, 
GeneWise, and FASTA. 



The source of the sequences may vary widely, and include taking sequences from one or more of the 
known databases, including, but not limited to, SCOP (Hubbard, et al., Nucleic Acids Res 27(1 ): 254- 
256. (1999)); PFAM (Bateman, etal., Nucleic Acids Res 27(1): 260-262. (1999) 
http://www.sanger.ac.uk/PfamO; TI6RFAM (http://www.tigr.org/TIGRFAMs); VAST (Gibrat, et al.. Cun- 
Opin Struct Biol 6(3): 377-385. (1996)); CATH (Orengo. etal.. Structure 5(8): 1093-1108. (1997)); 
PhD Predictor (http://www.embl-heidelberg.de/predictprotein/predictprotein.html); Prosite (Hofmann, 
et al.. Nucleic Acids Res 27(1): 215-219. (1999) http://www.expasy.ch/prosite/): SwissProf 
(http:/Awww.expasy.ch/sprotO; PIR (http7/www.mips.biochem.mpg.de/prqi/protseqdb/): GenBank 
(http:/Aiwww.ncbi.nlm.nih.gov/GenbankO; Entrez (http://www.ncbi.nlm.nih.gov/entrezO; RefSeq 
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{http:/ArtftArtw.ncbi.nlm.nih.gov/LocusLink/refseq.html); EMBL Nucleotide Sequence Database 
(http://www.ebi.ac.uk/embl/); DDBJ (http://www.ddbj.nig.ac.jp/); PDB (www.rcsb.org) and BIND 
(Bader, et al., Nucleic Acids Res 29(1 ): 242-245(2001 ) http://www.bind.ca/). In addition, sequences 
may be obtained from genome and SNP databases of organisms including, but not limited to, human, 
mouse, worm, fly, plants, fiingi, bacteria, and viruses. These may include public databases, for 
example The Genome Database of The Human Genome Project (http://gdbwww.gdb.org/), or private 
databases, for example those of Celera Genomics Corporation (http7/www.celera.com/) or Incyte 
Genomics (http://www.incyte.com/). 

In a preferred embodiment, ttie contribution of each aligned sequence to the frequency statistics is 
weighted according to its diversity weighting relative to other sequences in ttie alignment A common 
strategy for accomplishing this is the sequence weighting system recommended by Henikoff and 
Henikoff (see Henikoff S, Henikoff JG. Amino acid substitution matrices, Adv Protein Chem. 2000; 
54:73-97. Review. PMID: 10829225 and Henikoff S, Henikoff JG. Position-based sequence weights. J 
Mol Biol. 1994 Nov 4; 243(4): 574-8.PMID: 7966282). each are herein expressly incorporated by 
reference. 

In a preferred embodiment, only sequences within a preset level of homology to the template 
sequence are included in the alignment (> 60% identity, > 70% identity, etc.) 

In a preferred embodiment, the contribution of each sequence to the statistics is dependent on its 
extent of similarity to the target sequence, such that sequences with higher similarity to the target 
sequence are weighted more highly. Examples of similarity measures include, but are not limited to, 
sequence identity, BLOSUM similarity score. PAM matrix similarity score, and Blast score. 

In a preferred embodiment, the contribution of each sequence to the statistics is dependent on its 
known physical or functional properties. These properties include, but are not limited to, thermal and 
chemical stability, contribution to activity, solubility, etc. For example, when optimizing the target 
sequence for solubility, those sequences in an alignment with high solubility levels will contribute 
more heavily to the cak:ulated firequencies. 

In a preferred embodiment, each of the weighted or unweighted alignment frequencies is converted 
directly to a pseudo-energy as -log (fa). Thus, amino acids with higher frequency are assigned lower 
(more favorable) pseudo energies. If a frequency is zero, a constant positive pseudo energy may be 
applied. 

In a preferred embodiment, each of the final alignment frequencies (fa) is divided by the observed 
frequency (fo) of occurrence of each amino acid in all proteins. The log of this ratio, known to those in 
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the art as the log-odds ratio, tog(W, reflects the extent of natural selection for/against each ambto 
acid at each position in the proteir). Positive numl)ers reflect positive selection while negative 
numbers reflect negative selection. These log-odds ratios may ttien be used as pseudo energy terms 
within a PDA™ technology simulation. In situations where lower energies are favorable, the negative 
log-odds, -log(fA). is a more appropriate pseudo energy term, if a frequency is zero, a constant 
positive energy may be applied. 

In a prefierred embodiment, no pseudo energies are created. Rattier, the position-specific alignment 
information is used directty to generate the Bst of possible amino acids at a variable residue position In 
a PDA™ technology simulation. Lehnrann M. VWyss M. Engineering proteins for thermoslabnity: the 
use of sequence aOgnments versus rational design and directed evolution.Curr Opin Biotechnol. 2001 
Aug; 12(4): 371-5. Review; Lehmann M, Pasamontes L, Lassen SF, Wyss M. The consensus concept 
for ttiemiostability engineering of proteins. Biochim Biophys Acta. 2000 De c 29; 1 543(2): 408-415. 
Review; Rath A, Davidson AR. The design of a hyperstable mutant of the Abpl p SH3 domain by 
sequence alignment analysis.Protein Sd. 2000 Dec;9(12):2457-69; Lehmann M, Kostrewa D. Wyss 
M, Brugger R. D'Arcy A, Pasamontes L, van Loon AP. From DNA sequence to improved functionality: 
using protein sequence comparisons to rapidly design a flTermostable consensus phytase. Protein 
Eng. 2000 Jan;13(1):49-57; Desjarlais JR. Berg JM. Use of a zinc-finger consensus sequence 
framework and specificity rules todesign specific DNA binding proteins. Pfoc Nati Acad Sd U S A. 
1993 Mar 15;90(6):2256-60; Desjarlais JR. Berg JM. Redesigning the DNA-binding specificity of a 
zinc finger protein; a database-guided approach. Proteins. 1992 Feb;12(2):101-4; HenikoffS, Henikoff 
JG. Amino add substitution nratirlces. Adv Protein Chem. 2000; 54:73-97. Review. PWIID: 10829225; 
Henikoff S, Henikoff JG. Position-based sequence weights. J Mol Biol. 1 994 Nov 4 ; 243(4):574- 
8.PMID: 7966282. 



Similarly, structural alignment of stmcturally related proteins may be done to generate sequence 
alignments. There are a wide variety of such structural alignment programs known. See for example 
VAST from the NCBI (http://www.ncbi.nlm.nih.gov: 80/Stmcture/VAST/vast.shtml) ; SSAP (Orengo and 
Taylor, Mettwds Enzymol 266(617-635 (1996)) SARF2 (Alexandrov. Protein Eng 9(9): 727-732. 
(1996)) CE (Shindyalov and Boume, Protein Eng 11(9): 739-747. (1998)); (Orengo et al., Sh^icture 
5(8): 1093-108 (1997); Dali (Holm et al., Nudefc Acid Res. 26(1): 316-9 (1998). all of which are 
incorporated by reference). These sttudurally-generated sequence alignments may ttien be 
examined to determine ttie obsen/ed sequence variations. 

tn a prefen-ed embodiment, residue pair potentials may be used to score sequences (Miyazawa et al.. 
Macromolecules 18(3):534-552 (1985) Jones, Protein Science 3: 567-574. (1994); PROSA (Heindlich 
et al.. J. MoL BfoL 216:167-180 (1990); THREADER (Jones et al.. Nature 358:86-89 (1992). expressly 
incorporated by reference) during computational screening. 
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In a preferred embodiment sequence profile scores (see Bowie et al., Science 253(5016): 164-70 
(1991), incorporated by reference) and/or potentials of mean force (see Hendlich et al., J. Mol. Biol. 
216(1): 167-180 (1990). also incorporated by reference) are calculated to score sequences. 
Weighting using these methods determines the structural homology between the sequence and the 
three-dimensional structure of a reference sequence. These methods assess the match between a 
sequence and a three-dimensional protein structure and hence may act to screen sequences for 
fidelity to the protein structure. In particular, U.S. Patent Nos. 6,269,312. 6.188,965. and 6,403,312, 
and herein expressly incorporated by reference, describe a method termed "Protein Design 
Automation", or PDA""^ technology, that utilizes a number of scoring functbns to evaluate sequence 
stability. 

Primary libraries may be generated by predicting tertiary structure from sequence, and then selecting 
sequences that are compatible with the predicted tertiary structure. There are a number of tertiary 
structure prediction methods, including, but not limited to, threading (Bryant and Altschul, Curr Opin 
Struct Biol 5(2): 236-244. (1995)), Profile 3D (Bowie, et al., Methods Enzymol 266(598-616 (1996); 
MONSSTER (Skolnick, et al.. J Mol Biol 265(2): 21 7-241 . (1 997); Rosetta (Simons, et al., Proteins 
37(S3): 171-176 (1999); PSI-BLAST (Altschul and Koonin, Trends Biochem Sci 23(11): 444^7. 
(1998)); Impala (Schaffer, et al., Bioinformatics 15(12): 1000-1011. (1999)); HMMER (McClure, etal.. 
Proc Int Conf Intell Syst Mol Biol 4(155-164 (1996)); Clustal W (http://www.ebi.ac.uk/clustalw/) ; ), 
helix-coil transition theory (Munoz and Serrano, Biopolymers 41:495, 1997), neural networks, tocal 
structure alignment and others (e.g., see in Selbig et al., Bioinformatics 15:1039, 1999). 

In an alternate embodiment, the primary library consists of all sequences whose binary pattern, or 
arrangement of hydrophobic and polar residues, is predicted to be compatible with formation of the 
desired protein structure (Kamtekar et al., Science 262(5140): 1680-5 (1993). In an alternate 
embodiment two profile methods (Gribskov et al. PNAS 84:4355-4358 (1987) and Fischer and 
Eisenberg, Protein Sci. 5:947-955 (1996), Rice and Elsenberg J. Mol. Biol. 267:1026-1038(1997)). all 
of which are expressly incorporated by reference) are used to generate the primary library. 

In a further embodiment, a knowledge-based amino acid substitution matrbc can be used to guide the 
convergence of a protein design cycle. Examples of such nnatrices include, but are not limited to: 
BLOSUM matrices (e.g. 62, 90, etc.), PAM matrices (e.g. 250, etc.). and Dayhoff matrices. 

Energy Calculation Methods 

Force field calculations that may be used to optimize the conformation of a sequence within a 
computational method, such as molecular dynamics and rotamer placement methods, or to generate 
de novo optimized sequences as outlined herein. These methods can be used in any step of the 
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methods of the invention, including their use to generate a primary or secondary tibrary. 

Force fields include, but are not limited to, ab initio or quantum mechanical force fields, semi-empirical 
force fields, and molecular mechanics force fields. Examples of force fields include OPLS-AA 
(Jorgensen, etal., J. Am. Chem. Soc. (1996), v 118, pp 11225-11236; Joigensen, W.L; BOSS, 
Version 4.1; Yale University: New Haven. CT (1999)); OPLS (Jorgensen. et al., J. Am. Chem. Soc. 
(1988), V 110, pp 1657ff; Jorgensen. etal., J Am. Chem. Soc. (1990), v 112, pp 4768ff); UNRES 
(United Residue Forcefield; Liwo. et a!., Protein Science (1993). v 2. pp1697-1714; Liwo, et al., 
Protein Science (1993). v2, pp1715.1731; Lrwo. etal.. J. Comp. Chem. (1997). v 18, pp849-873; 
Liwo, et al., J. Comp. Chem. (1997), v 18, pp874-884; Liwo, et al., J. Comp. Chem. (1998), v 19, 
pp259-276; Forcefield for Protein Structure Prediction (Liwo, et al.. Proc. Natl. Acad. Sd. USA (1999). 
V 96, pp5482-5485); ECEPP/3 (Liwo et ai., J Protein Chem 1994 May; 13(4): 375-80); AMBER 1.1 
force field (Weiner, et al.. J. Am. Chem. Soc. v106, pp765-784); AMBER 3.0 force field (U.C. Singh et 
al., Proc. Natl. Acad. Sci. USA 82:755-759); CHARMM and CHARMM22 (Brooks, et al.. J. Comp. 
Chem. v4, pp 187-217); cvffS.O (Dauber-Osguthorpe. et al..(1988) Proteins: Structure. Function and 
Genetics, v4,pp31-47); cfTOI (Maple, et al., J. Comp. Chem. v15. 162-182); also, the DISCOVER (cvff 
and cfrai) and AMBER forcefields are used In the INSIGHT molecular modeling package 
(Biosym/MSI, San Diego California) and HARMM is used in the QUANTA molecular modeling 
package (Biosym/MSI, San Diego California). HF, UHF, MCSCF, CI. MPx. MNDO, AMI. and MINDO 
are techniques known to those skilled in the art and which may be used to perform computational site 
directed mutagenesis for protein design, (see Szabo et al. Modern Quantum Chemistry: Introduction 
to Advanced Electronic Structure Theory. Macmillan, New York. (c1982) and Hehre. Ab Initio 
Molecular Orbital Theory. Wiley. New York (c1986) ail of which are expressly incorporated by 
reference.) 

In a prefen-ed embodiment, the scaffold protein is an enzyme and highly accurate electrostatfc models 
may be used for enzyme active site residue scogng to improve enzyme active site libraries (see 
Warshel. Computer Modeling of Chemical Reactions in Enzymes and Solutions, Wiley & Sons, New 
York. 1991, hereby expressly incorporated by reference). These accurate models may assess the 
relative energies of sequences with high precisfon. but are computationally Intensive. Highly accurate 
electrostatrc models may also be used In the design of binding sites. 

Furthermore, scoring functions may be used to screen for sequences that would create metal or co- 
factor binding sites in the protein (Hellinga. Fold Des. 3(1): R1-8 (1998), hereby expressly 
incorporated by reference). Similarly, scoring functions may be used to screen for sequences that 
woukJ create disulfide bonds in the protein. 

In a preferred embodiment, rotamer library selection methods are used to generate the primary 
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library. (Dahiyat and Mayo» Protein Sci 5(5): 895-903 (1996); Dahiyat and Mayo» Science 278(5335): 
82-7 (1997); Desjarlais and Handel, Protein Science 4: 2006-2018 (1995); Harbury et al. PNAS USA 
92(18): 8408-8412 (1995); Kono et al., Proteins: Structure, Function and Genetics 19: 244-255 
(1994); Heinnga and Richards, PNAS USA 91 : 5803-5807 (1 994) . 

In a preferred embodiment, a sequence prediction algorithm (SPA) is used to design proteins that are 
compatible with a l<nown protein backbone stnjcture as is described in Raha, K., et al. (2000) Protein 
Sci., 9: 1106-1119, U.S.S.N. 09/877,695; USSN to be detemnined for a continuation-in-part 
application filed on February 6, 2002, entitled APPARATUS AND METHOD FOR DESIGNING 
PROTEINS AND PROTEIN LIBRARIES, with John R. Desjarlais as inventor, expressly incorporated 
herein by reference. 

In an altemate embodiment, other inverse folding methods such as those described by Simons et al. 
(Proteins, 34:535-543, 1999), Levitt and Gerstein (PNAS USA, 95:5913-5920, 1998). Godzik et al., 
PNAS. V89, PP 12098-102; Godzik and Skolnick (PNAS USA, 89:12098-102. 1992), Godzik et al. (J. 
Mol. Biol. 227:227-38, 1 992) may be used. 

In an altemate embodiment, molecular dynamics calculations may be used to computationally screen 
sequences by individually cak^latlng mutant sequence scores and compiling a rank ordered list 

In addition, other computational methods such as those described by Koehl and Levitt (J. Mot. BioL 
293:1161-1181 (1999); J. Mol. Biol. 293:1183-1193(1999); expressly incorporated by reference) may 
be used to create a primary library. 

PDA™ Technology Calculations 

In an especially preferred embodiment, the primary library is generated and processed as outlined in 
U.S. Patent Nos. 6,269,312, 6,188,965, and 6,403,312, and are herein expressly incorporated by 
reference. This processing step entails analyzing interactions of the rotamers with each other and 
with the protein backbone to generate optimized protein sequences. Simplistically, the processing 
initially comprises the use of a number of scoring functions to calculate energies of interactions of the 
rotamers. with the backbone and with other rotamers. Prefen^ed PDA™ technology scoring functions 
Include, but are not limited to, a van der Waals potential scoring function, a hydrogen bond potential 
scoring function, an atomic solvation scoring function, a secondary structure propensity scoring 
function and an electrostatic scoring function. As is further described below, at least one scoring 
function is used to score each variable or floated position, although the scoring functions may differ 
depending on the positbn classification or other consideratk)ns. 

As will be appreciated by those skilled in the art, a variety of force fields that may be used in the 
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PDA™ technology calculations. These include, but are not limited to. those listed previously. As 
outlined in U.S. Patent Nos. 6.269.312, 6,188.965, and 6,403,312. which are herein expressly 
incorporated by reference, any combination of the preferred scoring functions, either alone or in 
combination, may be used. For example. In an alternate embodiment, rotamer intemal energies are 
included. In additional embodimente, energies or scores that are a function of the confonnatlon 
and/or identity of three or wore amino acids are included. 

In further embodiments, additional terms are included to influence the energy of each rotamer state, 
including but not fimited to. reference energies, psuedo energies based on rotamer statistics, and 
sequence biases derived from multiple sequence alignments. Because sequence alignment 
information and rational methods have demonstrated utility for protein optimization, \he invention is an 
improvement via its combination of information from both methods. Sequence alignment Information 
atone may sometimes be misleading because of unfavorable couplings between amino acids that 
occur commonly in a multiple sequence alignment Rational methods alone, may have limitations, for 
example, are subject to systematic enters due to improper parameterization of force field components 
and weights. 

In a preferred embodiment, the scoring functions may be altered. Additfonal scoring functions may be 
used. Additional scoring functions include, but are not limited to torsional potentials, entropy 
potentials, additional solvation models including contact models, solvent exclusion models (see 
Lazaridis and Karplus. Proteins 35(2): 133-52 (1999)), and the like; and models for immunogenidty, 
(see U.S.S.N.s 09^903,378, 10/039.170, and PCT/US02/00165, herein expressly incorporated by 
reference) such as functions derived from data on binding of peptides to MHC (Major 
Histocompatibflity Complex), that may be used to Ideritify potentially immunogenic sequences. Such 
additional scoring functions may be used alone, or as functions for processing the library after it is 
initially scored. 

Altered scoring functions may also be obtained from analysis of experimental data. For example, if 
the presence of certain residues at certain positions are conrelated with tiie presence of desired 
protein properties, a scoring function may be generated which favor these certain residues. 

In addition, other metiiods may be used to "train" scoring functions by comparing designed sequences 
and their properties to natural sequences and their properties. That is, tiie relative importance, or 
weight, given to individual scoring functions can be optimized in a variety of ways. Although a variety 
of useful scoring functions exist tiiat represent van der Waals. electrostatics, solvation, and other 
terms, an important aspect of a force field is ttie contribution (or weight) of each scoring function to 
the total score. In a preferred embodiment, computational sequence screening may be used to 
identify force field parameters such that properties of natural proteins are mimicked in computationally 
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instances it may be desirable to include all sequences when a defined number of variable positions 
are used. It is usually preferable for the primary library to be small enough that a reasonable fraction 
of the sequence space of a particular sequence may be sampled, allowing for robust generation of 
secondary libraries. Thus, primary libraries that range from about 50 to 10^^ are preferred, with from 
1000 to 10^ being particularly prefen-ed, and from 1000 to 100,000 being especially preferred. Thus, 
In one preferred embodiment, the primary library excludes from 1% to about 90-95% of possible 
sequence space sequences, with exclusion of at least 1%, 2%, 5%, 10%. 20%, 40%, 50% and 70% 
being preferred. Altematively, the library may include 1 in 10^, 1 in 10^ 1 in 10^**. I ln 10^,1 in 10^, 
lin 10™and1in10®^ 

A variety of approaches may be used to select a set of sequences for the primary library, including 
structure-based methods such as PDA™ technology sequence-based methods, or combinations as 
outlined herein. In addition, as noted herein, any method used to generate a primary or secondary 
library may be used as the other step. 

It should also be noted that while these methods are described in conjunction with limiting the size of 
the primary library, these same techniques may be used to formulate a cutoff for inclusion in the 
secondary and tertiary libraries as well. 

The set of protein sequences in the primary and secondary libraries are generally, but not always, 
significantly different from the wild-type sequence from which the backbone was taken, although in 
some cases the primary or secondary library may contain the wild- type sequence. That is, the range 
of optimized protein sequences is dependent upon many factors including the size of the protein, 
properties desired, etc. However, for example, comprises between 0.001% and 100% variant amino 
acids, with about at least 90%. 70%. 50%. 30%. 10% variant amino acids being prefenred. 

In a preferred embodiment, the primary library sequences are obtained from a rank-ordered list or 
filtered set generated using an algorithm such as Monte Carlo, B&B, or SCMF. For example, the top 
10^ or the top 10^ sequences in the rank-ordered list or filtered set may comprise the primary library. 
Alternatively, all sequences scoring within a certain range of the optimum sequence may be used. 
For example, all sequences within 10 kcal/mol of the optimum sequence could be used as the 
primary library. In addition, as outlined below, any cut of a rank-ordered list or a filtered set may be 
used depending on the conditions, use and additional methodologies of the resulting set; for example, 
the top X number of sequences may be used, or the top X and the bottom Y number of sequences, 
for example when a wider range of sequence space is to be explored or when clustering is used. This 
method has the advantage of using a direct measure of fidelity to a three-dimensional structure to 
determine inclusion. 
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Alternatively, the total number of sequences defined by the recombination of all mutations may be 
used as a cutoff criterion for the primary sequence library. Preferred values for the total number of 
recombined sequences range from 100 to 10^", particulariy preferred values range from 1000 to 10" 
especially preferred values range from 1000 to 10^- /MtemaBvely. a cutoff may be enforced v^hen a 
predetermined number of mutations per position is reached. As a rank-ordered (or unordered) or 
filtered set sequence list is lengthened and the fibrary is enlarged . the number of mutations per 
position wrtll typically increase. Altematively. the first occurrence in the list of predefined undesirable 
residues may be used as a cutoff criterion. For example, the first hydiophfllc residue occurring in a 
core position could limit the set of sequences Included in the primary library. Alternatively when 
multiple related structures are used for the scaffold, the set of optimal sequences for each stmcture 
may be used to make ttie primary library. 

In addition, in some embodiments, sequences that do not make the cutoff are included in the primary 
library. This may be desirable in some situations, for instance to evaluate the primary library 
generation method, to sen/e as controls or comparisons, or to sample additional sequence space 
For example, in a preferred embodiment the wild-type sequence is Included, even if it dkl not make 
Uie cutoff. 



As IS further outiined betow. it should also be noted that different primary libraries may be combined 
For example, positions in a protein Uiat show a great deal of mutational diversity in computational 
screening may be fixed as outiined betow and a different primary library regenerated. A rank-ordered 
list or filtered set of the same lengtti as the first wouW now show diversity at positions tiiat were 
largely conserved in ttie first library. The variants from a first primary library may be combined witt, 
the vanants from a second primary library to provWe a combined Hbraiy at lower computational cost 
than creating a very long rank-ordered list or filtered set. This approach may be particulariy useful to 
sample sequence diversity in both highly mutatable and highly conserved positions. In addition 
primary libraries may be generated by combining ttie results of two or more calculations to form one 
primary library. 

CLUSTERIMG 

Clustering algoritiims may be useful for classifying sequences derived by protein design algorittims 
into representative groups. Clustering can serve a wide variety of purposes. For example, sets of 
sequences that are dose in sequence space can be distinguished from other sets, and thus 
recombination can be confined wittiin sets. That is, sequences tiiat share a local minima may be 
recombined. to allow better results, ratiier ttian recombine sequences from two local minima ttiat may 
have qurte different sequences. Thus, for example, a primary library can be clustered around local 
minima ("clustered sets of sequences'), recombination or secondary library generation is within each 
Clustered set. and tiien each -clustered- secondary library is added to fom. the secondary library 
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genus. 

Clustering algorithms require two key components. First is a metric for comparing ti^e similarity of two 
entities. Measures of similarity include, but are not limited to sequence identity, sequence similarity, 
and energetic similarity. Second, clustering algorithms require an algorithm to separate the entities 
Into groups based on relative similarities. Many types of clustering algorithms exist, the most simple 
and commonly used are single-linkage, complete linkage, and average linkage methods (see Figure 
5). These are often applied hierarchically, such that the relationships between entities may be 
described with a tree structure. 

Preferably, clustering algorithms including but not limited to, single linkage clustering algorithms, 
complete linkage clustering algorithms, and average linkage clustering algorithms are used to analyze 
the results from computational protein cycles described herein. Clustering algorithms may be used to 
form subsets using computationally generated energy matrices to measure energetic similarity (see 
Figure 6). Alternatively, clustering algorithms may be used to form subsets directly from a set of 
optimized protein sequences. 

In a preferred embodiment, a single-linkage clustering algorithm is used to form subsets from 
computationally generated energy matrk^s. An example of the use of a single-linkage clustering 
algorithm to form subsets from a computationally generated energy matrix is shown in Figures 5, 6. 
and 7. 

In altemative embodiments, a single linkage clustering algorithm is used to form subsets directly from 
a set of optimized protein sequences whereby the measure of similarity between two sequences is 
the extent of sequence identity. Alternatively, the measure of similarity between two sequences may 
be based on a standard sequence similarity comparison. As will be appreciated by those skilled in 
the art. similarity scores Include but are not limited to BLOSUM similarity score, Dayhoff similarity 
score, RAM similarity score, etc. Specific examples of the aforementioned similarity scores include but 
are not limited to BLOSUM tables, 62 and 90; RAM tables: 250, etc., among others. In a preferred 
embodiment, subsets of designed protein sequences derived by clustering or related methods may be 
used to define multiple primary or secondary libraries. 

In an altemate embodiment, sets of sequences that may be recomblned productively are defined as 
those that minimize disruption of sets of interacting or correlated residues. Identification of sets of 
interacting residues may be carried out by a number of ways, e.g. by using known pattern recognition 
methods, comparing frequencies of occurrence of mutations or by analyzing the calculated energy of 
interaction among the residues (for example. If the energy of interaction is high, the positions are said 
to be correlated or interacting). These correlations may be positional correlations (e.g. variable 
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residue positions 1 and 2 always change togetlier or never change together) or sequence cofrelations 
(e.g. if there is a residue A at position 1 . there is always residue D at position 2). In addition, 
programs used to search for consensus motifs may Ije used. See: Lockless and Ranganathan. 
Science 286:295-299 (1999). Pattern discovery in Biomolecular Data: Tools. Techniques, and 
Applications, edited by Jason T.L. Wang, Bruce A Shapiro. Dennis Shasha. New York: Oxfbrtl 
University. 1999; Andrews. Harry C. Introduction to mathematical techniques in patter recognition; 
New York. VWIey-lnterscience (1 972); ApplicaBons of Pattern Recognition; Editor. KS. Fu. Boca ' 
Raton, Fla. CRC Press. 1982; Genetic Algorithms for Pattern Recognition; edited by Sankar K. Pal. 
Paul P. W^ng. Boca Raton: CRC Press. c1996; Pandya. Abhijit S.. Pattem recognitfon with Neural 
networks in C++/AbhiiitS. Pandya. Robert B. Macy. Boca Raton. Fla.: CRC Press. 1996; Handbook 
of pattem recognition and computer vision / edited by C.H. Chen, LF. Pau, P.S,P. Wang. 2™* ed. 
Signapore;.River Edge, N.J. : World Scientific, c1999; and Friedman, Introduction to Pattem 
Recognition : Statistical. Structural, Neural, and Fuzzy Logic Approaches; River Edge, N.J. : World 
Scientific. c1 999. Series Title: Serien a machine perception and artificial intelligence; vol. 32. A II 
references cited herein are expressly Incorporated by reference. 

GENERATION OF SECONDARY LIBRARIES 

As described herein, there are a wide variety of methods to generate secondary libraries from primary 
libraries. The first is a selection step, where some set of primary sequences are chosen to form the 
secondary library. The second is a computational step, again generally including a selection step, 
where some subset of the primary library is chosen and then subjected to further computational ' 
analysis, including both protein design cycles as weH as techniques such as "in silico" shuffling 
(recombination). The third is an experimental step, where some subset of the primary library is 
chosen and then recombined experimentally to form a secondary library. 

SELECTING SEQUENCES FOR THE SECONDARY LIBRARY 

In a preferred embodiment, the primary library of the scaffoW protein is used to generate a secondary 
library. The secondary library may then be generated and tested experimentally or subjected to 
further computational manipulation. A variety of approaches, including but not limited to those 
described below, may be used to select sequences for the secondary library. Each approach may be 
used alone, or any combinafion of approaches may be used. As will be appreciated by ttiose in the 
art, the secondary library may be either a subset of the primary Hbrary, or contain new library 
members, i.e. sequences that are not found in the primary library. That is. in general, the variant 
positions and/or amino acid residues in the variant positions may be recombined in any number of 
ways to form a new library that exploits ttie sequence variations found in the primary library. In such 
embodiments, the secondary library will contain sequences ttiat were not included in the primary 
library. In all cases, if the secondary library is generated experimentally, it may optionally comprise 
one or more "error' sequences, which result from experinrental errors, as well as one or more 
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sequences generated intentionally. That is, additbnal variability can be added to the secondary (or, 
In fact, to the primary library as well), either experinientally (e.g. through the use of error-prone PGR 
in secondary library sequences) or computationally (adding an "in sllico" variant generation step to 
sample more sequence space). In the latter case, it is possible to introduce this additional level of 
variability in a random fashion (as used herein random includes variation introduced in a controlled 
manner or an uncontrolled manner) or in a directed fashion. For example, directed variability may be 
Introduced by adding certain residues from a particular sequence, e.g. the human sequence. 

Selecting a subset of the primary library 

As described herein, there are a wide variety of techniques that can be used to generate a secondary 
library. In a prefenred embodiment, a subset of the primary library is used as the secondary library. 
This subset can be chosen in a variety of ways, as outlined herein. For example, similar to the 
primary library cut-off, an arbitrary numerical cut-off can be applied: the top X number of sequences 
forms the basis of the secondary library (or the top X number and the bottom Y number, or any 
sequences in the top X number plus anything within Z energy of the wild- type sequence, etc. ). As will 
be appreciated by those in the art, there are a wide variety of relatively simple numerical cutoffs that 
can be applied. 

In a preferred embodiment, ail amino acid residues are allowed at each variable residue position 
Identified in the primary library. That is, once the variable residue positions are identified, a 
secondary library comprising every combination of every amino acid at each variable residue position 
is made. 

In a preferred embodiment, subsets of amino acids are chosen to maxinruze coverage. Additional 
amino acids with properties similar to those contained within the primary library may be manually 
added. For example. If tiie primary library includes three large hydrophobic residues at a given 
position, the user may chose to include additional large hydrophobic residues at that position when 
generating the secondary library. In addition, amino acids in the primary library that do not share 
similar properties with most of tiie amino acids at a given position may be excluded from tiie 
secondary library. Alternatively, subsets of amino acids may be chosen from the primary library such 
that a maximal diversity of side chain properties Is sampled at each position. For example, if the 
primary library includes three large hydrophobic residues at a given position, the user may chose to 
include only one of them In tiie secondary library, in combination with other amino acids that are not 
large and hydrophobic. 

In a prefen-ed embodiment, the primary library may be analyzed to determine which amino acid 
positions in the scaffold protein have a high mutational frequency, and which positions have a low 
mutation firequency. The secondary library may be generated by varying the amino acids at tiie 
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positions that have high numbers of mutattons. while Iceeping constant the positions that do not have 
mutations above a certain frequency. For example, if a position has less than 20% and more 
preferably less than 10% mutations, it may be held mvariant 

In a preferred embodiment, the secondary library is generated from a probability distribution table. As 
outlined herein, there are a variety of methods of generating a probability distribution table, including 
using PDA™ technology output, the result^ of other energy calculation methods, (e.g. SCiWF). and/or 
the results of knowledge- or sequence*ased methods, all described previously. In addition, the 
probability distribution may be used to generate information entropy scores for each position, as a 
measure of the mutational frequency observed in the library. In this embodiment, the frequency of 
each amino acid residue at each variable residue position in the list is identified. Frequencies may be 
thresholded, wherein any variant frequency lower than a cutoff is set to zero. This cutoff Is preferably 
1%, 2%, 5%. 10% or 20%, with 10% being particuiariy preferred. These frequencies may be buHt into 
the secondary library, so that the frequency at which each amino acid is present in the primary library 
is equal, within experimental error, to the firequency at which that amino acid will be present in the 
secondary library. 

Recombination of Some or All Primary Libra r y Sequences to Generate a Secondan/ Library 
In an alternate embodiment, variable residue positions may be recombined to generate novel 
sequences to form a secondary library. Thus, the secondary library comprises at least one member 
sequence and preferably a plurality of such member sequences not found in the primary library. 
Recombination may be performed experimentally and/or computationally using a variety of 
approaches. For example, a list of naturally occurring sequences may be used to calculate all 
possible recombinant sequences, with an optional rank ordering or filtering step. Alternatively, once a 
primary library is generated, one could rank order only those recombinations that occur at cross-over 
points with at least a threshold of identity over a given window (for example. 100% identic over a 
contiguous 18 nucleotide sequence, or 80% identity over a contiguous 24 nucleotide sequence). 
Alternatively, the homology could be considered at the DNA level, by computationally translaUng the 
amino adds to their respective DNAcodons. Different codon usages could be considered. A 
prefened embodiment considers only recombinations with crossover points that have DNA sequence 
identify suffident for hybridizatton. 

In some embodiments, all possible recombinant sequences are experimentally generated and tested. 
Altematively. in a preferred embodiment, the recombinant sequences are scored computationally and 
a subset of these sequences are experimentally generated and tested. Computational screening of 
the set of recombinant sequences may be used to reduce the library to an experimentally tractable 
size and/or to enrich the library in sequences predfcted to possess desired properties. The 
recombinant sequences may be analyzed using methods induding, but not restricted to, those 
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methods used to generate and analyze primary litxary sequences, and by considering the role of 
clusters of interacting residues, as discussed below. 

In a preferred embodiment, the secondary library in generated by using any of the techniques outlined 
for primary library generation (SPA, PDA™, taboo, clustering, "in sllico" recombination, etc.) on the 
primary library that has been chosen. Particular combinations of computational analyses for primary 
and secondary libraries are outlined below. 

In a preferred embodiment, the secondary library is generated experimentally, using any number of 
the techniques outlined below, including gene assembly procedures. 

It is possible that some recombinant sequences will be inviable. that is, they will fail to fold, aggregate , 
possess other undesired properties, or lacks desired properties. In certain cases, some algorithms 
will generate a plurality of local minima, the combination of which may lead to unsatisfactory 
sequences. 

However, computational screening approaches may be used to differentiate and bias or select for 
viable constructs from inviable constructs. For example. If recombining all library members is 
predicted to yield an excessive number of unviable sequences, subsets of a library could be 
recombined instead. Strategies for identifying sets of sequences ttiat may be productively 
recombined include, but are not limited to, clustering based on sequence identity or similarity, 
clustering based on similarity of the energy matrix, and identification of sets of interacting residues. 

As will be appreciated by those in the art and outiined herein, probability distribution tables can be 
generated in a variety of ways. In addition to the methods outlined herein, self -consistent mean field 
(SCMF) methods can be used in tiie direct generation of probabOity tables. SCMF is a deterministic 
computational method tiiat uses a mean field description of rotamer interactions to calculate energies. 
A probability table generated in this way can be used to create secondary libraries as described 
herein. SCMF can be used in tiiree ways: the frequencies of amino acids and rotamers for each 
amino acid are listed at each position; the probabilities are determined directiy from SCMF (see 
Delarue et la. Pac. Symp. Biocomput. 109-21 (1997), expressly incorporated by reference). In 
addition, highly variable positions and non-variable positions can be identified. Alternatively, another 
metiiod is used to determine what sequence is jumped to during a search of sequence space; SCMF 
is used to obtain an accurate energy for that sequence; tiiis energy is then used to rank it and create 
a rank-ordered list of sequences (similar to a Monte Carlo sequence list). A probability table showing 
the frequencies of amino acids at each position can then be calculated from this list (Koehl et al., J. 
MoL Biol. 239:249 (1994); Koehl et al., Nat. Struc. Biol. 2:163 (1995); Koehl et al.. Cun-. Opin. Struct. 
Biol. 6:222 (1996); Koehl et al., J. MoL Bio. 293:1183 (1999); Koehl et al., J. Mol. Biol. 293:1161 
(1999); Lee J. Mol. Biol. 236:918 (1994); and Vasquez Biopolymers 36:63-70 (1995); all of which are 
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expressly incorporated by reference. Other forcefields that can be used in similar methods are 
outlined above. 

In addition, as outlined herein, a preferred method of generating a probability distribution table is 
through the use of sequence alignment programs. In addition, the probability table can be obtained 
by a combination of sequence alignments and computational approaches. For example, one can add 
amino acids found In the alignment of homologous sequences to the result of the computation. 
Preferable one can add the wild type amino add Identity to the probability toble if it is not found in the 
computation. 

Generation of Tertiary Libraries 

In a prefenred embodiment, a variety of additional steps may be done to one or more secondary 
libraries; for example, further computational processing may occur, secondary libraries may be 
recombined, or subsets of different secondary libraries may be combined. 

In a prefenred embodiment a tertiary library can be generated from combining secondary libraries. 
For example, a probability distribution table from a secondary library can be generated and 
recombined, whether computationally or experimentally, as outlined herein. A PDA secondary library 
may be combined with a sequence alignment secondary library, and either recombined (again, 
computationally or experimentally) or just the cutoffs from each joined to make a new tertiary library. 
The top sequences from several libraries can be recombined. Primary and secondary libraries can 
similariy be combined. Sequences from the top of a library can be combined with sequences from the 
bottom of the library to more broadly sample sequence space, or only sequences distant from the top 
of the library can be combined. Primary and/or secondary libraries that analyzed different parts of a 
protein can be combined to a tertiary library that treats the combined parts of the protein. These 
combinations can be done to analyze large proteins, especially large multidomain proteins or 
complete protoesomes. 

In a preferred embodiment, a tertiary library can be generated using con-elations in the secondary 
library. That is, a residue at a first variable position may be correlated to a residue at second variable 
position (or correlated to residues at additional positions as well). For example, two variable positions 
may sterically or electrostatically interact, such that if the first residue is X, the second residue must 
be Y, This may be either a positive or negative correlation. This con-elation, or "cluster" of residues, 
may be both detected and used in a variety of ways. (For the generation of correlations, see the 
earlier cited art). 

In addition, primary and secondary libraries can be combined to form new libraries; these can be 
random combinations or the libraries, combining tiie "top" sequences, or weighting the combinations 
(positions or residues from the first library are scored higher than those of the second library). 
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Additional variability can be added to the terfiary library as well), either experimentally (e.g. through 
the use of error-prone PGR in tertiary library sequences) or computationally (adding an "in silico" 
variant generation step to sample more sequence space). In the latter case, It is possible to introduce 
this additional level of variability In a random fashion (as used herein random Includes variation 
Introduced in a controlled manner or an uncontrolled manner) or in a directed fashion. For example, 
directed variability may be introduced by adding certain residues from a particular sequence, e.g. the 
human sequence. 

In a prefen^ed embodiment, when two computational steps are used (e.g. a PDA™ step to generate a 
primary library and in sifico shuffling or a probability table to ge nerate a secondary library), the 
experimental generation of the secondary library can result in a tertiary library, that is. a library that 
contains members not found in the secondary library. Alternatively, the tertiary library may just be a 
subset of the secondary library as outlined above. 

In a prefen-ed embodiment, a secondary library may be computationally remanlpulated to form an 
additional secondary library (sometimes referred to herein as "tertiary libraries"). For example, any of 
ttie secondary library sequences may be chosen for a second round of PDA™ technology 
calculations, by freezing or fixing some or all of the changed positions in the first secondary library. 
Alternatively, only changes seen in the last probability distribution table would be allowed. 
Alternatively, the stringency of the probability table may be altered, either by increasing or decreasing 
the cutoff for inclusion. 

In a preferred embodiment, the sequence information derived from experimental screening of a 
secondary library could be used to guide the design for the tertiary library. In this way, the library 
generation is an iterative process. In a prefenred embodiment, the tertiary library could be derived by 
computationally screening the secondary library for desired protein properties as previously 
mentioned. 

Experimentally Making the Lfbraiv 

Once a library is generated using any of the methods outlined herein or combinations thereof, the 
library (or a tertiary, quaternary, etc. library) Is made any number of techniques, Including using gene 
assembly procedures. Accordingly, the present invention provides methods for making protein 
libraries in any of a variety of different ways. 

Chemical synthesis of proteins 
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In a preferred embodiment, different protein members of tfie secondary library may be chemically 
synthesized. This is particularly useful when the designed proteins are short, preferably less than 150 
amino adds in length, with less than 100 amino adds being preferred, and less than 50 amino adds 
being particularly preferred, although as is known in the art. longer proteins may be made diemically 

or enzymatically. 

These amino add sequences could then be joined together via chemical ligation to form larger 
proteins as needed (see Yan. L and Dawson. P.E. J. Am. Chem. Soc. 123 (2001) 526-533. and 
Dawson. P.E. and Kent. S.B.H. Ann. Rev. Biochem. 69, (2000) 923-960), hereby expressly 
incorporated by reference. Furthermore, peptides corresponding to sequences from different library 
members could be shuffled or randomly ligated together to form a secondary library. For example, 
one or more peptides with different amino acid sequences firom the N-terminai region of the protein 
could be rigated to one or more peptides with different amino acid sequences from the C-terminal 
region of the protein. Such an assembly could be repeated for several further rounds of synthesis. 
Using such a method, a secondary library could be chemically synthesized. 

In a preferred embodiment, proteins could be constructed by chemically synthesis of peptides and 
fomied by ligation of the peptides using intein technology (Evans et al. (1999) J. Biol. Chem. 274, 
18359-18363; Evans etal. (1999) J. Biol. Chem. 274, 3923-3926; Mathys etal. (1999) Gene 231* 1- 
13; Evans et al. (1998) Protein Sci. 7.2256-2264; Southworth et al. Biotechniques 27. 1 10-120). 



Generating nucleic aci ds that encode single members of a library 

In a preferred embodiment, particularly for longer proteins or proteins for which large samples are 
desired, the secondary library sequences are used to create nucleic adds sudi as DNA which encode 
the member sequences and whidi may then be cloned into host cells, expressed and assayed, if 
desired. Thus, nudeic acids, and particularly DNA. may be made which encodes each member 
protein sequence. This is done using well-known procedures. See Maniatis and current protocols, 
(see Current Protocols In Molecular Biology, Wiley & Sons, and Molecular Ctoning - A Laboratory 
Manual - 3"" Ed. . CoW Spring Harbor laboratory Press, New York (2001)). The choice of codons. 
suitable expression vectors and suitable host cells will vary depending on a number of factors, and 
may be easily optimized as needed. 

Gene Assembhr Procedures 

As will be appredated by those in the art. the generation of exact sequences for a library comprising a 
large number of sequences (despite the fact that the set number is much smaller than the original set) 
is still potentially expensive and time consuming. Accordingly, in a preferred embodiment, there are a 
variety of gene assembly tediniques that may be used to generate the secondary or higher order 
libraries of the present invention. As discussed herein, these experimentally generated libraries 
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generally recombine sequences within the library, resulting in sequences present in the original library 
as well as recombined combinations of those sequences. 

Gene Assembly Using Pooled Qligonucieottdes 

In a preferred embodiment, multiple amplification reactions with pooled oligonucleotides are done, as 
is generally depicted in Figure 12, comprising variant protein sequences created by the assembly of 
gene fragments generated from a nudeic acid template. This generally involves generating variant 
protein sequences created by the assembly of gene fragments generated from a nucleic acid 
template. They can be full length "overlapping" oligonucleotides, or primers. In one embodiment, 
overlapping oligonucleotides are synthesized which correspond to the full-length gene. As may be 
appredated by one skilled In the art, these oligonucleotides nnay represent all of the different amino 
acids at each variant position or subsets. Once these oligonucleotides are made, they are 
reassembled into a set of variable sequences in any number of ways, outlined below. While the 
reactions described below focus on PGR as the amplification techniques, others are induded as is 
generally outlined below. 

In general, the invention may take on a wide variety of configuratbns. For example, libraries of 
nucleic acids encoding all or a subset of possible proteins are generated by assembling nudeic add 
fragments. Preferably, the gene fragments are linked together using an enzymatic or non-enzymatic 
method for the ligation of gene fragments. For example, for each gene fragment, a pair of donor 
fragments is generated such that the sense strand from one donor fragment complements the 
antisense strand of the other donor fragment and creates a 5'-phosphorylated overhang when the two 
strands are hybridized under conditions that allow for the formation of a double stranded molecules. 
The 5' phosphorytated overhang Is located at one of the 5' ends of the resulting double stranded 
molecule to allow ligation to a free 3'- terminus of an adjacent gene fragment. In some embodiments, 
5'-phosphorylated overhangs are generated at both ends, preferably with unique sequences to 
prevent setf-ligation. 

Chemically synthesized oligonucleotides are used as primers for the generation of donor fragments. 
For each pair of donor fragments, one primer is labeled at the 5*-end with a purification tag. The 
purification tag may be a his, myc. flag, or HA tag or a fusion protein may be used instead, for 
example gst, thioredoxin, nusA, among others known in the art. Preferably, the purification tag is 
biotin. The other primer is designed to bind to the other member of the donor fragment pair to create 
a 5 -phosphorylated overhang, from about 1 to 20 or more base pairs in length. 

In a preferred embodiment, at least one of the populations of nucleic acid fragments comprise variant 
sequences that result in the formation of a variant nucleic acid sequence. In a further embodiment, 
both the 5 'phosphorylated primer and at least one of the populations of nucleb acid firagments are 
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used to generate variant nucleic add sequences. In a preferred embodiment, ligation substrates are 
formed from at least two different donor fragment pairs. The donor fragment pairs may be generated 
from the same template or from different templates. 

In a preferred embodiment, the ligation product is generated using the following steps: (1) generating 
at least two donor fragments from a template molecule using primer dependent DNA polymerization 
wherein one strand comprises a purification tag and the other strand comprises a 5'-phosphorylated 
• overhang; (2) removing strands tagged with a purifrcation tag using a suitable capture molecule; (3) 
annealing the remaining 5'-phosphorylated strand to form first and second ligation substrates; and. (4) 
ligating said first and second ligation substrates after annealing strands with 6* phosphorylated 
overhangs to generate nucleic acid molecules encoding variant proteins, (see Kneidinger, Graininger 
and Messner, Biotechniques 30: 249-252 (2001); Au. Yang, Yand, Lo, and Kao; Biochem Blophys 
Res Comm 248: 200-203 (1998)). Each of the above-cited references are herein expressly 
incorporated by reference. This method is more fully described in U.S. Pat No. 6.1 10.668 and 
W09815567. 

In a prefen-ed embodiment, the donor fragments are generated using modified primers and a 
polymerase. The nucleic acid template may be single stranded (i.e. Ml 3 DNA) or double stranded 
(i.e., plasmid, genomic, or cDNA). The overall design of the primers will depend on the linkage 
scheme between the donor fragments. For example, (Rgure 20) Illustrates the controlled linkage 
between two neighboring fragments A and B. Initially for each gene fragment, a pair of donor 
fragments is generated (DFA1/DFA2 and DFB1/DFB2). The donor fragment pairs are designed such 
that the sense strand from one donor fragment, DFA1 or DFB1 , complements the antisense sfrand of 
the other donor fragment, DFA2 or DFB2. and creates a S'nahosphorylated overhang on the hybrid 
product of the con-esponding two strands. The overhang is located on the side where two 
neighboring gene fragments are to be joined. The sequence of the overhang is a sequence that 
belongs either to the 3'-end of fragment A or the S'-end of fragment B (in Figure 1 0. it belongs to B 
FIX). The sfrands not used to form the sticky end hybrid molecule are removed using a purification 
tag. 

In a preferred embodiment, the strands not used to form the sticky end are removed using 
biotin/streptavidin capture technology as is known in the art In an alternative embodiment, a 5'- 
phosphorylated primer is incorporated on the strand to be removed, followed by digestion of this 
strand with lambda exonuclease. Subsequent 5'-phosphorylation of the remaining sfrand will allow 
formation of a hybrid molecule with a phosphorylated overhang. 

In a prefenred embodiment, equimolar amounts of the con^esponding single sfrands of the donor 
fragments are combined under conditions suitable to renature double stranded molecules (A/A' and 
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B/B*), with a 5 -phosphorylated overhang. Preferably, these double stranded molecules, also referred 
to herein as ligatbn substrates are joined using enzymatic or non enzymatic ligation to fbnm a nucleic 
acid ligation product that encodes a protein variant Altematively, the ligation substrate is not ligated, 
but instead is used as a source of donor fragments and the process repeated. 

The following U.S. patents are incorporated herein in their entirety: U.S. Patent No. 6,188.965; U.S. 
Patent No. 6,269,312; and U.S. Patent No. 6,403,312. The following U.S. patent applications are 
incorporated herein in their entirety: U.S.S.N. 09/927.790, filed August 10, 2001 and U.S.S.N. 
10/101 ,499. filed March 1 8. 2002. 

In a preferred embodiment, the oligonucleotides are pooled in equal proportions and multiple PCR 
reactions are performed to create full length sequences containing the combinations of mutations 
defined by the secondary library. In addition, this may be done using methods that introduce 
addib'onal variations, such as error-prone amplification (e.g. PCR) methods or by intentionally 
introducing other variables. 

In a preferred embodiment, the different oligonucleotides are added in relative amounts 
con-esponding to either a probability distribution table or to an arbitrary or computatbnally derived 
formula. The multiple PCR reactions thus result in full length sequences with the desired 
combinations of mutations in the desired proportions. 

The total number of oligonucleotides needed is a function of the number of positions being mutated 
and the number of nnutations being considered at these positions: 

(number of oligos for constant positions) + M1 + M2 + M3 + Mn = (total number of ollgos required), 
where Mn is the number of mutations considered at position n in the sequence. 
In a preferred embodiment, each overlapping oligonucleotide comprises only one position to be 
varied; in alternate embodiments, the variant positions are too cbse together to allow this and multiple 
variants per oligonucleotide are used to allow complete recombination of ail the possibilities. That is, 
each oligo can contain the codon for a single position being mutated^ or for more than one position 
being mutated. The multiple positions being mutated must be close In sequence to prevent the oligo 
length from being impractical. For multiple mutating positions on an oligonucleotide, particular 
combinations of mutations can be included or excluded in the library by including or excluding the 
oligonucleotide encoding that combination. For example, as discussed herein, there may be 
con-elatlons between variable regions; that is, when position X is a certain residue, position Y must (or 
must not) be a particular residue. These sets of variable positions are sometimes referred to herein 
as a "cluster". When the clusters are comprised of residues close together, and thus can reside on 
one oligonuclotide primer, the clusters can be set to the "good" correlations, and eliminate the bad 
combinations that may decrease the effectiveness of the library. However, if the residues of the 
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Cluster are far apart in sequence, and thus wOl reside on diffiBrent oligonuctatides for synthesis, it may 
be desirable to either set the residues to the "good" correlation, or eliminate them as variable residues 
entirely. In an altemative embodiment, the library may be generated in several steps, so that the 
cluster mutations only appear together. This procedure, i.e.. the procedure of dentifying mutation 
clusters and either placing them on the same oligonucleotides or eliminating them from the library or 
library generation In several steps preserving clusters, can considerably enrich the experimental 
library with properly folded protein. Identification of clustere can be carried out by a number of ways. 
e.g. by using known pattern recognition methods, comparisons of fifequencies of occurrence of 
mutattons or by using energy analysis of the sequences to be experimentally generated (Ibr example, 
if the energy of interaction is high, the positions are correlated), these correlations may be positional ' 
correlations (e.g. variable positions 1 and 2 always change together or never change together) or 
sequence correlations (e.g. if there is a residue A at position 1. there is always residue B at position 
2). See: Pattern discovery in Biomolecular Data: Tools. Techniques, and Applications; edited by 
Jason T.L Mteng. Bruce A. Shapiro. Dennis Shasha. New York: Oxford Unviersity. 1 999; Andrews. 
Harry C. Introduction to mathematical techniques in patter recognition; New York, Wiley -Interscience 
[1972]; Applications of Pattern Recognition; Editbr. KS. Fu. Boca Raton. Fla. CRC Press, 1982; 
Genetic Algorithms for Pattern Recognttton; edited by Sankar K Pal. Paul P. \Nang. Boca Raton : 
CRC Press. c1996; Pandya. Abhijit S.. Partem recognition with Neuial networlcs in C++/Abhijit S. 
Pandya. Robert B. Macy. Boca Raton. Fla.: CRC Press. 1996; Handbook of pattern recognition and 
computer vision / edited by C.H. Chen. LP. Pau. PS.P. Wang. 2nd ed. Signapore ; River Edge. N.J. 
: World Scientific. c1999; Friedman. Introduction to Pattern Recognition : Statistical. Structural. 
Neural, and Fuzzy Logic Approaches ; River Edge. N.J. : Wortd Scientific. c1999. Series title: Serien a 
machine perception and artifldal intelligence; vol. 32; all of wh ich are expressly incorporated by 
reference. In addition programs used to search for consensus motifs can be used as weB. 

In addition, correlations and shuffling can be fixed or optimized by altering the design of the 
ongonucleotkJes; that is. by deciding where the oligonucleotides (primers) start and stop (e.g. where 
the sequences are "cuT). The start and stop sites of oligos can be set to maximize the number of 
clusters that appear in single oligonucleotides, thereby enriching the library with higher scoring 
sequences. Different oligonucleotkies stert and stop site opttons can be compuiatnnally modeled 
and ranked according to number of clusters that are represented on single oligos. or the percentage 
of the resulting sequences consistent with the predicted libarary of sequences. 

The total number of oligonucleotides required increases when multiple mutable posittons are encoded 
by a singte oligonucleotide. The annealed regions are the ones that remain constant. i.e. have the 
sequence of the reference sequence. 

Oligonucleotkies with insertfons or deletions of codons can be used to create a library expressing 
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different length proteins. In particular computational sequence screening for insertions or deletions 
can result in secondary libraries defining different length proteins, which can be expressed by a library 
of pooled oligonucleotide of different lengths. 

Preferably, an individual gene that serves as the template nucleic acid is obtained from at least two 
different species. In this embodiment, the gene from one species is ctoned into a vector to produce a 
template molecule comprising single stranded nucleic acid molecules. The DNAfrom the second 
species is cleaved into firagments. The resulting fragments are added to the template molecule under 
conditbns that permit the fragments to anneal to the template molecule. Unhybridized termini are 
enzymaticaliy removed. Gaps between hybridized fragments are filled using an appropriate enzyme, 
such as a polymerase and nicks sealed using a ligase. The chimeric gene can be amplified using 
suitable primers or other techniques that are well known to those of skill in the art. 

In a preferred embodiment, sequences derived from introns are used to mediate specific cleavage 
and ligation of discontinuous nucleic acid molecules to create libraries of novel genes and gene 
products as described in U.S. Patent Nos. 5,498,531, and 5,780,272, both of which are hereby 
expressly incorporated by refierence in their entirety. In one embodiment a library of ribonucleic acids 
encoding a novel gene product or novel gene products is created by mixing spficing constructs 
comprising an exon and 3' and 5' intron fragments. See U.S. Patent No. 5,498,531. 

In another embodiment, DNA sequence libraries are created by mixing DNA/RNA hybrid molecules 
that contain intron derived sequences that are used to mediate specific cleavage and ligation of the 
DNA/RNA hybrid molecules such that the DNA sequences are covalentiy linked to form novel DNA 
sequences as described in U.S. Patent No. 6,150.141, WO 00/40715 and WO 00/17342, all of which 
are hereby expressly incorporated by reference in their entirety. 

In a preferred embodiment, the secondary library is done by shuffling the family (e.g. a set of 
variants); that is, some set of the top sequences (if a rank-ordered list is used) can be shuffled, either 
with or without error-prone PGR. "Shuffling" in this context means a recombination of related 
sequences, generally in a either a targeted or random way. It can Include "shuffling" as defined and 
exemplified in U.S. Patent Nos. 5.830.721; 5,811.238; 5,605.793; 5,837,458 and PCT US/19256, all 
of which are expressly incorporated by reference in their entirety. This set of sequences can also be 
an artificial set; for example, from a probability table (for example generated using SCMF) or a Monte 
Carlo set Similarly, the ''iamWf can be the top 10 and the bottom 10 sequences, the top 100 
sequences, etc. This may also be done using enror-prone PGR. 

Thus, in a preferred embodiment, in sllico shuffling is done using the computational methods 
described therein. That is, starting with either two libraries or two sequences, random recombinations 
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Of the sequences can be generated and evaluated computationally, and then experimental libraries 
generated. 

PGR with pooled oHaos 

Use of pooled oligos for synthetic shuffling is more fully described In U.S. Pat No. 6.368.861 (see 
also US6423542; US6376246: US6368861; US6319714; W00042561A3; WO0042561A2; 
WO0042560A3; WO0042560A2; WO0042559A1; WO0018906C2; WO0018906A3; and 
WO0018906A2.) 

In a preferred embodiment, PGR using a wild type gene or other gene may be used, as Is 
schematically depicted in Figure 16. In this embodiment, a starting gene is used: the gene may the 
wild-type gene, the gene encoding the global optimized sequence, or any other sequence of the list. 
In this embodiment, oligonucleotides are used that correspond to the variant positions and contain the 
different amino acids of the secondary library. PGR Is done using PGR primers at the temiini, as is 
known in the art PGR provides many benefits namely, fewer oligonucleotides, may result in fewer 
errors, and if the wild type gene is used, it need not be synthesized. An alternative method for 
creating members of the library, are llgase chain reaction-based methods, (see Chalmers and 
Cumow, Biotechniques 30 (2001) 249-252). which in herein expressly incorporated by reference. 

In a preferred embodiment, these oligonucleotides are pooled in equal proportions and multiple PGR 
reactions are performed to create fuD-length sequences containing the combinations of mutations 
defined by the secondary library. In a prefenred embodiment, the different oligonucleotides are added 
in relative amounts, e.g. in amounts con-esponding to a probability distribution table, an alignment, or 
other parameters. The multiple PGR reactions thus result in full-length sequences with the desired 
combinations of mutations in the desired proportions. 

Number of mutations per oilqo 

In a preferred embodiment, each overiapping oligonucleotide comprises at least one or more 
positions to be varied and zero or more positions that are not varied. As may be appreciated by one 
skilled in the art, the distance between multiple variants may affect the completeness of 
recombination of all possible library members. That is. each oligo may contain the codon for a single 
position being mutated, or for more than one position being mutated. For multiple mutating positions 
on an oligonucleotide, particular combinations of mutations may be included or excluded in the library 
by including or excluding the oligonucleotide encoding that combination. The total number of 
or^onucleotides required increases when multiple mutable positions are encoded by a single 
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oligonucleotide. The annealed regions are the ones that remain constant, l.e. have the sequence of 
the reference sequence. 

Random codons 

In some cases, oligos with random mutatbns may be used. That is, any amino acid may be 
represented at a codon position. As known by those skilled in the art, subsets of random codons may 
be used, where the bias is for or against specific amino acids. By judicial design, certain amino ackls 
may be favored or excluded from the set of possible mutations. 

Multiple DNA libraries may be synthesized that code for different subsets of amino acids at certain 
positions, allowing generation of the amino acid diversity desired without having to fully randomize the 
codon and thereby waste sequences in the library on stop codons, frameshifts, undesired amino 
acids, etc. This may be done by creating a library that at each position to be randomized is only 
randomized at one or two of the positions of the triplet where the positlon(s) left constant are those 
that the amino ackJs to be considered at this position have in common. Multiple DNA libraries rnay be 
created to insure that all amino acids desired at each position exist in the aggregate library. 
Alternatively, shuffling, as is generally known In the art, may be done with multiple libraries. 
Alternatively, the random peptide libraries may be done using the frequency tabulation and 
experinnental generation methods Including, multiplexed PGR, shuffiing, and the like. 

Error-prone PGR 

In a preferred embodiment, error-prone amplification methods (e.g. error prone PGR) is done to 
generate additional members of the secondary library, or the whole library. See U.S. Patent Nos. 
5,605,793, 6,811,238, and 5,830,721, all of which are hereby incorporated by reference. This may be 
done on the optimal sequence or on top members of the library, or some other artificial set or family. 
Enror prone PGR is then performed on the optimal sequence gene in the presence of oligonucleotides 
that code for the mutations at the variable residue positions of the secondary library (bias 
oligonucleotides). The addition of the oligonucleotides will c reate a bias favoring the incorporatbn of 
the mutations in the secondary library. Alternatively, only oligonucleotides for certain mutations may 
be used to bias the library. 

In addition to enror-prone PGR, mutations could be introduced in specific regions using minor 
modifications to several other methods, either in vitro or In vivo, including but not limited to "DNA 
shuffling" (see WO 00/42561 A3: WO 01/70947 A3;), exon shuffling (see US 8365 377 B1; Kolkman & 
Stammer (2001) Nature Biotechnology 19, 423-428), family shuffling (see Grameri etal. (1998) 
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Nature 391, 288-291; US 6376246 B1), RACHITT™ (Coco et al. (2001) Nature Bfotechnology 19, 
354- 359; WO 02/06469 A2). STEP and random priming of in vitro recombination ( see Zhao et al., 
(1998) Nature Biotechnology 16. 258-261; Shao et al (1998) Nucleic Acids Research 26, 681-683; 
exonudease mediated gene assembly (US 6352842 B1, US 6361974 81), Gene Site Saturation 
Mutagenesis™ (US 6358709 B1). Gene Reassembly^ (US 6368709B1) and SCRATCHY (Lutz et 
al.(2001), PNAS 98. 11248-11253). DNA fragmentation methods (Kikuchietal.. Gene 236. 159-167). 
single-stranded DNA shuffling (Kikuchi et al, (2000) Gene 243, 133-137). Although these methods 
are Intended to introduce random mutations throughout the gene, those skilled In the art will 
appreciate that specific regions (those defined by computational methods such as PDA™ technology: 
see WO 01/75767) of the gene couW be mutated, whilst others could be left untouched, either by 
isolating and combining the mutated region with the unmodified region (for example, by cassette 
mutagenesis; see WO 01/75767 A2; Kim & Mass, (2000) Biotechnkjues 28. 196-198; Lank) & Jeltsch 

(1998) Biotechniques 25, 958- 965; Ge & Rudolph (1997) Biotechniques 22. 28- 30; Ho et al., (1989) 
Gene 77, 51059), or via in vitro or in vivo recombination (see for example see WO 02/10183 Al and 
Abecassis et al.. (2000) Nucleic Acids Research 28. e88 for examples). All of the above-cited 
references are hereby expressly incorporated by reference. In addition, it should be noted that the 
compulational equivalents of all of these methods can be used as a computatfonal step to generate 
primary and/or secondary Obraries. That is. "in silico' shuffling of a primary library rank -ordered list 
may be further "shuffled" using experimental procedures. 

Additional methods for gene construction 

The creation of members of the secondary library may be performed by several other methods, 
including, but not limited to, classical site-directed mutagenesis, e.g. Quickchange commercially 
avanable from Stratagene. cassette mutagenesis as well as other amplification techniques. Cassette 
mutagenesis could include the creation of DNA molecules firom restriction digestion fragments using 
nucleic ackJ ligation, and includes the random ligation of restriction fragments (see Kikuchi et al., 

(1999) , Gene 236, 159-167). Additionally, cassette mutagenesis could also be achieved using 
randomly-cleaved nucleic acids (see Kikuchi et al.. (1999), Gene 236, 133-137), by PCR-ligation PGR 
mutagenesis (see for example Ali & Steinkasserer (1995), Biotechniques 18, 746-750), by seamless 
gene engineering using RNA- and DNA- overhang cloning (see Roc & Doc; Coljee et al.. (2000) 
Nature Bfotechnology 18, 789-791). by ligatfon mediated gene construction (U.S.S.N. 60/311,545), by 
homotogous or non-homologous random recombination (see US6.368.861; US6423642; US6376246; 
US6368861; US6319714; WO0042661A3; WO0042661A2; WO0042560A3; WO0042660A2; 
WO0042559A1; WO0018906C2; WO0018906A3; and WO0018906A2), or in vivo using 
recombination between flanking sequences (see WO 02/10183 A1 and Abecassis et al., (2000) 
Nucleic Acids Research 28. e88 for examples). In addition, regions of the gene could be mutated in E 
coli lacking correct mismatch repair mechanisms, (e.g. E.coli XUnutS strain commercially available 
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from Stratagene), or by using phage display techniques to evofve a library (e.g. Long-McGie et al., 

(2000) . Biotechnol Bioeng 68. 121-125). 

In addition to the PGR methods outlined herein, there are other amplification and gene synthesis 
methods that can be used. For example, the library genes may be "stitched" together using pools of 
oligonucleotides with polymerases (and optionally or solely) ligases. These resulting variable 
sequences can then be amplified using any number of amplification techniques, including, but not 
limited to. polymerase chain reaction (PGR), strand displacement amplification (SDA), nucleic acid 
sequence based amplification (NASBA), ligation chain reaction (LCR) and transcription mediated 
amplification (TMA). In addition, there are a number of variations of PGR which may also find use in 
the invention, including "quantitative competitive PGR" or "QG-PGR", "arbitrarily primed PGR" or "AP- 
PGR" . "immuno-PGR", "Alu-PGR", "PGR single strand conformational polymorphism" or "PGR- 
SSCP'. "reverse transcriptase PGR" or "RT-PGR". "biotin capture PGR", "vectorette PGR", "panhandle 
PGR", and "PGR select cDNA subtration", among others. Furthermore, by incorporating the 17 
polymerase initiator into one or more oligonucleotides, IVT amplification can be done. 

Experimental Itflodification of Libraries to Generate Further Libraries 

It will be appreciated by those skilled in the art that many of the methods used to construct the 
secondary libraries can be used in further modifications. For example, cassette mutagenesis could 
include the creation of DNA molecules from restriction digestion fi-agments using nucleic acid ligation, 
and includes the random ligation of restriction fragments (see KIkuchi et al., (1999), Gene 236, 159- 
167). Additionally, cassette mutagenesis could also be achieved using randomly -cleaved nuclerc 
ackJs (see Kikuchi et al., (1999), Gene 236. 133-137), by PCR-ligation PGR mutagenesis (see Ail & 
Steinkasserer (1995), Blotechniques 18. 746-750), by seamless gene engineering using RNA- and 
DNA- overhang cloning (Roc & Doc; Goljee et al., (2000) Nature Biotechnology 18, 789-791), by 
ligation mediated gene construction (U. S.S.N. 60/311,545), by homologous or non-homologous 
random recombination (see US6,368,861: US6423542; US6376246; US6368861; US6319714; 
WO0042561A3; WO0042561 A2; WO0042560A3; WO0042560A2; WO0042559A1; WO0018906G2: 
WO0018906A3; and WO0018906A2). 

Tertiary libraries could be created from secondary libraries using any of the techniques outlined herein 
or one or more of the following, either in a step-wise fashion or in combination: DNA shuffling ( see 
WO 00/42561 A3; WO 01/70947 A3;), exon shuffling (see US 6365 377 B1; Kolkman & Stemmer 

(2001) Nature Biotechnology 19, 423-428), Family Shuffling (see Grameri et al. (1998) Nature 391, 
288-291; US 6376246 B1), RAGHITTtm (see Coco et al. (2001) Nature Biotechnology 19. 354- 359; 
WO 02/06469 A2), STEP and random priming of in vitro recombination (see Zhao et a!., (1998) 
Nature Biotechnology 16, 258-261; Shao et al (1998) Nucleic Acids Research 26, 681-683; 
exonuclease mediated gene assembly (see US 6352842 B1. US 6361974 B1), Gene Site Saturation 
MutagenesisTM (see US 6358709 B1), Gene Reassembly^" (see US 6358709B1) and SGRATGHY 
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(see Lutzetal.(2001). PNAS98. 11248-11253), DNA fragmentation methods (see Kikuchi et al.. 
Gene 236. 159-167), single-stranded DNA shuffling (see Kikuchi et al., (2000) Gene 243, 133-137), ii 
vitro or in vivo recombination (see WO 02/10183 Al and Abecassis et al.. (2000) Nucleic Acids 
Research 28, e88 for examples). Additionally, in vivo mutagenesis could be perfornned in strains of 
£co// that lack conrect DNA mismatch repair mechanisms, e.g. E.coli XLmutS strain commerdally 
available from Stratagene, or by using phage display techniques to evolve a library (e.g. Long-McGie 
etaL. (2000), Biotechnol Bloeng 68, 121-125). 

Preferred Combinations 

In general, as more fully outlined below, the invention can take on a wide variety of configurations. In 
general, primary libraries, e.g. libraries of all or a subset of possible proteins are generated 
computationally. This can be done in a wide variety of ways, including sequence alignments of 
related proteins, structural alignments, stmctural prediction models, databases, or (preferably) protein 
design automation computational analysis. Similarly, primary libraries can be generated via 
sequence screening using a set of scaffoW structures that are created by perturbing the starting 
structure (using any number of techniques such as molecular dynamics, Monte Carlo analysis) to 
make changes to the protein (including backbone and sidechain torsion angle changes). Optimal 
sequences can be selected for each starting structures (or, some set of the top sequences) to make 
primary libraries. 

Some of these techniques result in the list of sequences in the primary library being "scored", or 
"ranked" on the basis of some particular criteria. In some embodiments, lists of sequences that are 
generated without ranking can then be ranked using techniques as outlined below. 

In a preferred embodiment some subset of the primary library is then experimentally generated to 
form a secondary library. Alternatively, some or all of the primary library members are recombined to 
form a secondary library, e.g. with new members. Again, this may be done either computationally or 
experimentally or both. 

Alternatively, once the primary library is generated, it can be manipulated in a variety of ways. In one 
embodiment, a different type of computational analysis can be done; for example, a new type of 
ranking may be done. Alternatively, and the primary library can be recombined. e.g. residues at 
different positions mixed to forni a new, secondary library. Again, this can be done either 
computationally or experimentally, or both. 
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As will be appredated by those in the art, there are a number of specific combinations that can be 
used with the methods of the present invention. Examples of some preferred combinations are 
shown in Figures 21 A-E. 

Expression Systems 

The library proteins of the present invention are produced by culturing a host cell transformed with 
nucleic acid, preferably an expression vector, containing nucieic acid encoding a library protein, under 
the appropriate conditions to induce or cause expression of the library protein. The conditions 
appropriate for library protein expression will vary with the choice of the expression vector and the 
host cell, and will be easily ascertained by one skilled in the art through routine experimentation. For 
example, the use of constitutive promoters in the expression vector will require optimizing the growth 
and proliferation of the host cell, while the use of an inducible promoter requires the appropriate 
growth conditions for induction. In addition, in some embodiments, the timing of the harvest is 
important. For example, the baculoviral systems used In insect cell expression are lytic viruses, and 
thus harvest time selection can be crucial for product yield. 

Examples of expression systems 

As will be appreciated by those in the art, the type of cells used in the present invention can vary 
widely. The lists that follow are applicable both to the source of scaffold proteins as well as to host 
cells in which to produce the variant libraries. A wide variety of appropriate host cells can be used, 
including yeast, bacteria, archaebacteria, fungi, and insect, plant and animal cells, including 
mammalian cells. Of particular interest are Drosophila melanogaster cells, Saccharomyces cerevisiae 
and other yeasts, E, coli, Bacillus subtilis. Streptococcus cremoris, Streptococcus IMdans, pED 
(commercially available from Novagen), pBAD and pCNDA (commercially available from Invitrogen), 
pEGEX (commercially available from Amersham Biosciences), pQE (commercially available from 
Qiagen), SF9 cells, C129 ceils, 293 cells, Neurospora, BHK, CHO. COS, and HeLa ceils, fibroblasts, 
Schwanoma cell lines, immortalized mammalian myeloid and lymphoid eel! lines, Jurkat cells, mast 
cells and other endocrine and exocrine cells, and neuronal cells. See the ATCC cell line catalog, 
hereby expressly incorporated by reference. In one embodiment, the cells may be genetically 
engineered, that is, contain exogenous nucleic add, for example, to contain target molecules. 

In a preferred embodiment, the library proteins are expressed in mammalian expression systems, 
including systems in which the expression constructs are introduced into the mammalian cells using 
virus such as retrovirus or adenovirus. Any mammalian cells may be used, with mouse, rat, primate 
and human cells being particularly prefen-ed, although as will be appreciated by those in the art, 
modifications of the system by pseudotyping allows all eukaryotic cells to be used, preferably higher 
eukaryotes. Accordingly, suitable marnmalian cell types include, but are not limited to, tumor cells of 
all types (particularly melanoma, myeloid leukemia, carcinomas of the lung, breast, ovaries, colon. 
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lodney. prostate, pancreas and testes), cardiomyocytes. endotheGal cells, epithelial cells, lymphocytes 
(T-cells and B cells) . mast cells, eosinophils, vascular intimal cells, hepatocytes. leukocytes including 
mononuclear leukocytes, stem cells such as haemopoetic. neural, skin. lung, kidney, fiver and 
myocyte stem cells (for use in screening for differentiation and de^lifferentiation factors), osteoclasts, 
chondrocytes and other connective tissue cells, keratinocytes, melanocytes, liver cells, kidney cells ' 
and adipocytes. Suitable cells also include known research cells, including, but not limited to. Jurkat 
T cells. NIH3T3 cells. CHO. Cos. etc. Again. scaffbU proteins may be obtained from these sources 
as well. 



In a preferred embodiment, library proteins are expressed in bacterial systems, including bacteria in 
which the expression constructs are introduced into the bacteria using phage. Bacterial expressfon 
systems are well known iij the art. and include Badttus subWis. £ coB. Stfeptococcus cnmoris. and 
Streptococcus IMdans 

In an alternate embodiment, library proteins are produced in insect ceils, including but not limited to 
Dmsophila melanogaster S2 cells, as well as cells derived from members of the order Lepidoptera 
which includes all butterflies and moths, such as the silkmofh Bombyx mon and the alphalpha tooper 
Autographa califomica. Lepklopteran insects are host organisms for some members of a family of 
virus, known as baculoviruses (more than 400 known species), that inlect a variety of arthropods 
(see U.S. 6,090,584). 

In an alternate embodiment, fibrary proteins are produced in insect cells. Tlie library can be 
transfected into SF9 Spodoptera fmgiperda insect ceUs to generate baculovirus which are used to 
infect SF21 or High Rve commercially available fiom Invitrogen. insect cells for high level protein 
productfon. Also, transfections into the DrosopMa SchneWer S2 cells will express proteins. 

In a preferred embodiment, library protein is produced in yeast cells. Yeast expression systems are 
well known in the art. and Include expression vectors for Sacchammyces cerevisiae, Candida 
albicans and C. maltosa. Hansenula pofymorpha, Kluyveromyces fragHis and K. lactis. Pichia 
guHlerimondii and P. pastoris, Schizosacchammyces pombe, and Yarrowia lipotytica. 

In one embodiment the library proteins are expressed in vitro using cell free translation systems. 
Several commercial sources are available for this including but not limited to Roche RapW Translatfon 
System. Promega TnT system, Novagen's EcoPro system. Ambton's ProteinScipt-Pro system in 
vibv translation systems derived firom both prokaryotic (e.g. E. coli) and eukaryotic (e.g. Wheat germ 
Rabbit reticulocytes) cells are available and can be chosen based on the expression levels and 
functfonal properties of the protein of interest. Both linear (as derived from a PGR amplificatton) and 
circular (as in plasmkJ) DNA molecules are suitable for such expresston as long as they contain the 
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gene encoding the protein operably linked to an appropriate promoter. Other features of the nnolecule 
that are important for opttmal expression in either the bacterial or eukaryotic cells (including the 
ribosome binding site etc) are also included in these constructs. The proteins can again be 
expressed individually or in suitable size pools consisting of multiple library members. The main 
advantage offered by these in vitro systems is their speed and ability to produce soluble proteins. In 
addition the protein being synthesized can be selectively labeled if needed for subsequent functional 
analysis. 

Transformation and transfection methods 

The methods of introducing exogenous nucleic acid into host cells is well known in the art, and will 
vary with the host cell used. Techniques include dextran-mediated transfection, calcium phosphate 
precipitation, calcium chloride treatment, potybrene mediated transfection, protoplast fusion, 
electroporation, viral or phage infection, encapsulation of the polynucleotide(s) in liposomes, and 
direct microinjection of the DNA into nuclei. In the case of mammalian ceils, transfection may be 
either transient or stable. 

Expression Vectors 

A variety of expression vectors may be utilized to express the library proteins. The expression 
vectors are constructed to be compatible with the host ceil type. Expressbn vectors may comprise 
self-replicating extrachromosomal vectors or vectors which Integrate into a host genome. Expression 
vectors typically comprise a library member, any fusion constructs, control or regulatory sequences, 
selectable markers, and/or additional elements. 

Preferred bacterial expression vectors include but are not limited to pET. pBAD. bluescript, pUC. 
pQE. pGEX. pMAL, and the like. 

Preferred yeast expression vectors include pPICZ, pPIC3.5K, and pHIL-SI commercially available 
from Invitrogen. 

Expression vectors for the transformation of insect cells, and in particular, baculo virus -based 
expression vectors, are well known in the art and are described e.g.. in O'Reilly et al.. Baculovirus 
Expression Vectors: A Laboratory Manual (New York: Oxford University Press, 1994). 

A preferred mammalian expression vector system is a retroviral vector system such as is generally 
described in Mann et al., Cell, 33:153-9 (1993); Pear et al.. Proc. Natl. Acad. Sci. U.S.A.. 
90(18):8392-6 (1993); KItamura et al., Proc. Natl. Acad. Sci. U.S.A., 92:9146-50 (1995); Kinsella et 
al., Human Gene Therapy. 7:1405-13; Hofmann eta!., Proc. Natl. Acad. Sci. U.S.A., 93:5185-90; 
Choate et al., Human Gene Therapy. 7:2247 (1996); PCT/US97/01019 and PCT/US97/01Q48. and 
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references cited therein, all of which are hereby expressly incorporated by reference, 
inctusio n of control or reaulatofv sequences 

Generally, expression vectors include transcriptional and translational regulatory nucleic add 
sequences which are operably linked to the nucleic add sequence encoding the library protein. 

The transcriptional and translational regulatory nucleic acid sequences will generally be appropriate to 
the host cell used to express the library protein, as will be appredated by those in the art For 
example, transcriptional and translational regulatory sequences from £ coli are preferably used to 
express proteins in £ colL 

Transcriptional and translational regulatory sequences may include, but are not limited to. promoter 
sequences, ribosomal binding sites, transcriptional start and stop sequences, translational start and 
stop sequences, and enhancer or activator sequences. In a prefen-ed embodiment, the regulatory 
sequences comprise a promoter and transcriptional and translational start and stop sequences. 

A suitable promoter is any nucleic ackJ sequence capable of binding RNA polymerase and initiating 
the downstream (3') transcription of the coding sequence of library protein into mRNA Promoter 
sequences may be constitutive or indudble. The promoters may be naturally occurring promoters, 
hybrid or synthetic promoters. 

A suitable bacterial promoter has a transcriptbn initiation region which is usually placed proximal to 
the 5' end of the coding sequence. The transcription initiation region typically includes an RNA 
polymerase binding site and a transcription initiation site. In £ co//. the ribosome-blnding site is called 
the Shine-Dalgarno (SD) sequence and includes an initiation codon and a sequence 3-9 nucleotides 
in length located 3-11 nucleotides upstream of the initiation codon. Promoter sequences for 
metabolic pathway enzymes are commonly utilized. Examples include promoter sequences derived 
from sugar metabolizing enzymes, such as galactose, lactose and maltose, and sequences derived 
from biosynthetic enzymes such as tryptophan. Promoters from bacteriophage, such as the T7 
promoter, nnay also be used. In addition, synthetic promoters and hybrid promoters are also useful; 
for example, the fac promoter is a hybrid of the trp and lac promoter sequences. 

Prefen-ed yeast promoter sequences indude the indudble GALI.IO promoter, the promoters from 
alcohol dehydrogenase, enolase, glucokinase. glucose-6-phosphate isomerase, glyceraldehyde-3- 
phosphate-dehydrogenase, hexokinase, phosphofructokinase. 3-phosphoglycerate mutese. pyruvate 
kinase, and the acid phosphatase gene. 

A suitable mammalian promoter will have a transcription initiating region, which is usually placed 
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proximal to the 5' end of the coding sequence, and a TATA box, usually located 25-30 base pairs 
upstream of the transcription initiation site. The TATA box is thought to direct RNA polymerase II to 
begin RNA synthesis at the correct site. A mammalian promoter will also contain an upstream 
promoter element (enhancer element), typically located within 100 to 200 base pairs upstream of the 
TATA box. Typically, transcription termination and polyadenylation sequences recognized by 
mammalian cells are regulatory regions located 3' to the translation stop codon and thus, together 
with the promoter elements, flank the coding sequence. The 3' terminus of the mature mRNA Is 
formed by site-specific post-translational cleavage and polyadenylation. Examples of transcription 
terminator and polyadenylation signals include those derived from SV40. An upstream promoter 
element determines the rate at which transcription is initiated and can act in either orientation. Of 
particular use as mammalian promoters are the promoters from mammalian viral genes, since the 
viral genes are often highly expressed and have a broad host range. Examples include the SV40 
earty promoter, mouse mammary tumor vims LTR promoter, adenovims major late promoter, herpes 
simplex virus promoter, and the CMV promoter. 

Inclusion of a selectable marker 

In addition, in a preferred embodiment, the expression vector contains a selectbn gene or marker to 
allow the selection of transformed host cells containing the expression vector. Selection genes are 
well known in the art and will vary with the host cell used. 

For example, a bacterial expression vector may include a selectable marker gene to allow for the 
selection of bacterial strains that have been transformed. Suitable selection genes include genes 
which render the bacteria resistant to drugs such as ampicillin, chloramphenicol, erythromycin, 
kanamycin, neomycin and tetracycline. 

Yeast selectable maricers include the biosynthetic genes ADE2, HIS4, LEU2, and TRP1 when used in 
the context of auxotrophe strains; ALG7, which confers resistance to tunlcamycin; the neomycin 
phosphotransferase gene, which confers resistance to G418; and the CUP1 gene, which allows 
yeast to grow in the presence of copper ions. 

Suitable mammalian selection markers include, but are not limited to, those that confer resistance to 
neomycin (or its analog G418), blastbidin S, histinkiol D, bleomycin, puromycin, hygromydn B, and 
other drugs. Selectable mariners confenring survivability in a specific media include, but are not limited 
to Blasticidin S Deaminase. Neomycin phophotranserase II, Hygromydn B phosphotranserase, 
Puromycin N-acetyl transferase, Bleomycin resistance protein (or Zeocin resistance protein, 
Phleomycin resistance protein, or phleomycin/zeocin binding protein), hypoxanthine guanosine 
phosphoribosyl transferase (HPRT), Thymidylate synthase, xanthine-guanine phosphoridosyl 
transferase, and the like. 
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tnctusion of additional elements 

In addition, the expression vector may comprise additional elennents. In a preferred embodiment, the 
vector contains a fusion protein, as discussed below. In another embodiment, the expression vector 
may have two replication systems, thus allowing it to be maintained in two organisms, for example in 
mammalian or insect cells for expression and in a prokaryotic host for cloning and amplification. 
Furthemnore, for integrating expression vectors, the expression vector contains at least one sequence 
homologous to the host cell genome, and preferably two homologous sequences which flanic the 
expression construct The integrating vector may be directed to a specific locus in the host cell by 
selecting the appropriate homologous sequence for inclusion in the vector. Such vectors may include 
cre-lox recombination sites, or affR. aftB. affP, and attL sites. Constructs for integrating vectors and 
appropriate selection and screening protocols are well known in the art and are described in e.g., 
Mansour et al., Ceil, 51:603 (1988) and Murray, Gene Transfer and Expression Protocols, Methods in 
Molecular Biology, Vol. 7(C!iflon: Humana Press, 1991). In a preferred embodiment, the expression 
vector contains a RNA splfcing sequence upstream or downstream of the gene to be expressed in 
order to increase the level of gene expression.. (See Barret etal., Nucleic Adds Res. 1991; Grooset 
aL, mi Cell. Bk>l. 1987; and Budiman etal.. Mol. Ceil. Btol. 1988.) 

Fusion Constructs 

The library protein may also be made as a fusion protein, using techniques well known in the art. For 
example, fusion partners such as targeting sequences can be used which allow the localization of the 
library members into a subcellular or extracellular compartment of the cell. Purification tags may be 
fused with a library, allowing the purification or isolation of the library protein. Rescue sequences can 
be used to enable the recovery of the nucleic acids encoding them. Other fusion sequences are 
possible, such as fusions which enable utilization of a screening or selection technology. 

Targeting or signal sequences 

The expression vector may also include a signal peptide sequence tiiat directs library protein and any 
assodated fusions to a desired cellular tocation or to ttie extracellular media. Suitable targeting 
sequences indude, but are not limited to, binding sequences capable of causing binding of the 
expression product to a predetermined molecule or dass of molecules while retaining bioactivity of 
the expression product, (for example by using enzyme inhibitor or substi-ate sequences to target a 
dass of relevant enzymes); sequences signalling selective degradation, of itself or co-bound proteins; 
and signal sequences capable of constitutively localizing the candidate expression products to a 
predetemnihed cellular locale, induding a) subcellular locations such as tiie Go^i. endoplasmic 
reticulum, nudeus. nudeoli. nudear membrane, mitochondria, chloroplast. secretory vesicles, 
lysosome, and cellular membrane; and b) extracellular locations via a secretory signal. Target 
sequences also may be used in conjunction with cell surface display technology as discussed below. 
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Particularly preferred is localization to either subcellular locations or to the outside of the cell via 
secretion. For example some targeting sequences enable secretion of library protein in bacteria. The 
signal sequence typically encodes a signal peptide comprised of hydrophobic amino acids which 
direct the secretion of the protein from the cell, as is well known in the art. This method may be useful 
for gram-positive bacteria or gram-negative bacteria. The protein can be either secreted into the 
growth media or into the periplasmic space, located between the inner and outer membrane of the 
cell. 

Purification tags 

In a preferred embodiment, the library member comprises a purification tag operably linked to the rest 
of the library peptide or protein. A purification tag is a sequence which may be used to purify or 
isolate the candidate agent, for detection, for Immunoprecipitation, for FACS (fluorescence- activated 
ceil sorting), or for other reasons. Thus, for example, purification tags include purification sequences 
such as polyhistidine, including but not limited to Hise, or other tag for use with Immobilized Metal 
Affinity Chromatography (IMAC) systems (e.g. Ni*^ affinity columns), GSTfijsions, MBP fusfons, 
Strep-tag, the BSP biotinylation target sequence of the bacterial enzyme BIrA, and epitope tags which 
are targeted by antibodies. Suitable epitope tags include but are not limited to c-myc (for use with the 
commercially available 9E10 antibody), fiag tag, and the like. 

Rescue fusions 

A rescue fusion is a fusion protein which enables recovery of the nucleic acid encoding the library 
protein. In a preferred embodiment, such a rescue fusion would enable screening or selection of 
library members. Such fusion proteins may include but are not limited to, rep proteins, viral VPg 
proteins, transcription factors Including but not limited to zinc fingers, RNA and DNA binding proteins, 
and the like. Attachment can be covalent or noncovalent 

Alternatively, the rescue sequence may be a unique oligonucieotide sequence that serves as a probe 
target site to allow the quick and easy isolation of the retroviral construct, via PCR. related 
techniques, or hybridization. 

In an alternate embodiment, rescue sequences could also be based upon in vivo recombination 
systems, such as the cre-lox system, the Invitrogen Gateway system, forced recombination systems 
in yeast, mammalian, plant, bacteria or fungal cells (see WO 02/10183 A1), or phage display systems. 

In an alternate embodiment, display technologies are utilized. For example, in phage display (see 
Kay, BK et al, eds. Phage display of peptides and proteins: a laboratory manual (Academic Press. 
San Diego, CA, 1996); Lowman HB. Bass SH, Simpson N, Wells JA (1991) Selecting high-affinity 
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binding proteins by monovalent phage display. Bioechennistry 30:10832-10838; Smith GP (1985) 
Filamentous fusion phage: novel expression vectors that display cloned antigens on the virion 
surface. Science 228:1315-1317.) Bbrary proteins can be fused to the gene III protein. Cell surfece 
display (Witrrup KD. Protein engineering by celi-surlace display. Curr. Opin. Biotechnology 2001, 
12:395-399.) may also be useful for screening. This includes but is not limited to display on bacteria 
(see Georgiou G. Poetschlte HL, Stathopoulos C. Francisco JA, Practical applications of engineering 
gram-negative bacterial cell surfaces. Trends Biotechnol. 1993 Jan;1 1(1):6-10; Georgiou G, 
Stethopoulos C. Daugherty PS. Nayak AR. Iverson BL, and Curtiss RR (1997) Display of 
heterologous proteins on the surface of microorganisms: from the screening of combinatorial libraries 
to live recombinant vaccines. Nature Biotechnol. 15, 29-34; Lee JS. Shin KS. Pan JG. Wm CJ. 
Surface-displayed viral antigens on Salmonella earner vaccine. Nature Biotechnology, 2000 , 18:845- 
648; Jun et al, 1998). yeast (see Boder ET, Witlrup KD: Yeast surface display for screening 
combinatorial polypeptde libraries. Nat Biotechnol 1997, 15:553-557 Boder ET and Wrttrup KD. Yeast 
surface display for directed evolution of protein expresston, affinity, and stability. IWethods Enzymol 
2000, 328:430-44.). and mammalian cells (see Whitehom EA, Tate E, Yanofslqr SD. Kochersperger 
L. Davis A. Mortensen RB. Yonkovlch S. Bell K. Dower WJ. and Barrett RW 1995. A generic method 
for expression and use of "fagged" soluble verstons of cell surface receptors. Bio/lechnology, 13. 
1215-1219.). 



Additional fusions that allovy for screening or selection 

In an alternate embodiment, a protein fragment complementation assay is used (see Johnsson N & 
Varshavsky A. Split UbiqulBn as a sensor of protein interacHons in vivo. 1994 Proc Natl Acad Sci 
USA. 91: 10340-10344; Peltetier JN. Campbell-Valois FX, Michnick SW. Oligomerization domain- 
directed reassembV of active dihydroiblate reductase firom rationally designed fragments. 1998. Proc 
Natl Acad Sci USA 95:12141-12146.) Other fusion methods which may allow screening include but 
are not limited to periplasmic expression and cytometric screening (see Chen G. Hayhurst A, Thomas 
JG, Harvey BR, Iverson BL. Georgiou G: Isolation of high-affinity ligand-blnding proteins by 
periplasmic expresston with cytometric screening (PECS). Nat Biotechnol 2001. 19: 537-542.). and 
the yeast two hybrid screen (see ReWs S. Song O: A novel genetic system to detect protein-protein 
interactions. Nature 1989, 340:245-246.) 

Other fusions 

Additional fusion partners may also be utilized. For example, library protein may be made as a fusion 
protein to increase expression, increase solubility, confer stability or protection from degradation, 
and/or confer other properties. For example, when raising monoclonal antibodies to a small epitope, 
the library protein may be fused to a earner protein to form an immunogen. According to 
Varshavsky's N-End Rule, susceptibBity to ubiquHinatton and subsequent degredatfon can be 
minimized by the incorporation of glycines after the initiation nnethtonine (MG or MGG), thus 
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conferring bng half-life in tlie cytoplasm. Simitarly, adding two prolines to flie C-terminus confers 
resistance to carboxypeptidase action. 

Linkers 

Linker sequences may be used to connect tlie library protein to its fusion partner or tag. The linker 
sequence will generally comprise a small number of amino acids, typk:ally less than ten. However, 
longer linkers may also be used. As will be appreciated by those skilled in the art, any of a wide 
variety of sequences may be used as linkers. Typically, linker sequences are selected to be flexible 
and resistant to degradation. A common linker sequence comprises the amino acid sequence 
GGGGS. The preferred linker between a protein and C-terminal PP tag consists of two glycines. 

Labels 

In one embodiment, the library nucleic acids, proteins and antibodies of the invention are labeled. In 
general, labels fall into three classes: a) Immune labels, which may be an epitope incorporated as a 
fusion constructs nfiay which Is recognized by an antibody as discussed above, isotopic labels, whrch 
may be radioactive or heavy Isotopes, and c) small molecule labels which may include fluorescent 
and cotorlmetric dyes or molecules such as blotin which enable the use of other labeling techniques. 
Labels may be incorporated into the compound at any position and may be incorporated in vivo during 
protein or peptide expression or in vitro. 

Protein Purification 

In a preferred embodiment, the library protein is purified or isolated after expression. Library proteins 
may be isolated or purified in a variety of ways known to those skilled in the art depending on what 
other components are present in the sample. The degree of purificatton necessary will vary 
depending on the use of the library protein. In some instances no purification will be necessary. For 
example in one embodiment, if library proteins are secreted, screening or selection can take place 
directly from the media. 

Standard purification methods include electrophoretic, molecular, immunological and chromatographic 
techniques, including ion exchange, hydrophobic, affinity, size exclusion chromatography, and 
reversed-phase HPLC chromatography, as well as precipitation, dialysis, and chromatofocusing 
techniques. Purification can often be facilitated by the inclusion of purification tag, as described 
above. For example, the library protein may be purified using glutathione resin if a GST fusion is 
employed, Immobilized Metal Affinity Chromatography (IMAC) if a His or other tag is empbyed, or 
immobilized anti-flag antibody if a flag tag is used. Ultrafiltration and diafiltration techniques, in 
conjunction with protein concentration, are also useful. For general guidance in suitable purification 
techniques, (see Scopes, R., Protein Purification: Principles and Practice 3'^^ Ed., Springer-Verlag, NY 
(1994).), hereby expressly incorporated by reference. 
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In a preferred embodiment, the libraries are used In any number of display techniques. For example 
the libraries may be displayed using phage or enveloped virus systems, bacterial systems, yeast two 
hybrid systems or mammalian systems. 

In a prBfierred embodiment, the libraries are displayed using a phage or enveloped virus system. For 
example, a library of viruses, each canying a distinct peptide sequence as part of the coat protein 
can be produced by inserting random oligonucleotides sequences into the coding sequence of viral 
coat or envelope proteins. Several different viral systems have been used to display peptides as 
described in Smith. G.P. (1985) Science. 228:1315-1317; Santini. C. et al.. (1998) J. Mol. Biol.. 
282:125-135; Sternberg. N. and Hoess, R.H. (1995) Proc. Natl. Acad. Sci. USA. 92:1609-1613-' 
IWaruyama. I.N.. elal. (1994) Proc. NaO. Acad. Sci. USA, 91:8273-8277; Dunn. I.S.. (1995) J. Mol 
Bfol.. 248:497-506; Rosenberg. A., etal. (1996) Innovations 6:1-6); Ren. Z.J., etal. (1996) Protein 
Sci.. 5:1833-1843; Bimov. V.P., etal. (1995) Virus Genes 10:173-177; Dulbecco. R.. U.S Patent No 
4.593.002; Ladner, R.C... et al.. U.S. Patent No. 5.837.500; Ladner, R.C.. et al.. U.S. Patent No 
5.223.409; Dower, et al.. U.S. Patent No. 6.427.908; Russell et ai;. U.S. Patent No. 5.723.287; Li U.S. 
Patent No. 6,190.856; and the application entitled "METHODS AND COMPOSITIONS FOR THE 
CONSTRUCTION AND USE OF ENVELOPE VIRUSES AS DISPLAY PARTICLES", filed August2 
2001 . serial number not yet assigned, ail of which are expressly incorporated by reference. 

In a prefen-ed embodiment, the libraries are displayed on the surface of a bacterial cell as is described 
in WO 97/37025. which is expressly incorporated by reference in Hs entirety. In this embodiment 
surface anchoring vectors are provided for the surface expression of genes encoding proteins of* 
interest At a minima, the vector includes a gene encoding an ice nucleation protein, a secretion 
signal a targeting signal and a gene of interest Preferably, the bacterial host is a gram negative 
bacterium belonging to the genera Escherichia. Acetobacter, Pseudomonas. Xanthomonas. Erwinia 
and Xymomonas. Advantages to using the ice nucleation protein as the surface anchoring protein ara 
the high level of expression of the ice nudeatton protein on the surface of the bacterial ceil and its 
stable expression during the stationary phase of bacterial cell growth. 

In a preferred embodiment the libraries ara displayed using yeast two hybrid systems as is described 
.n Fields and Song (1989) Nature 340:245. which is expressly incorporated herein by reference 
Yeast-based two-hybrid systems utilize chimeric genes and detect protein-protein interactions via the 
activabon of reporter-gene expression. Reporter-gene expression occurs as a result of reconstitution 
of a functional transcription factor caused by the association of fusion proteins encoded by the 
Chimeric genes. Preferably, the yeast two-hybrid system commercially available from Clontech is 
used to screen libraries Ibr proteins that interact with a candidate proteins. See generally. Ausubel et 
aL. Current Protocols in Molecular Biology. John Wiley & Sons, pp.13.14.1.13.14.14. which is 
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expressly Incorporated herein by reference. 

In a preferred embodiment, the libraries are displayed using mammalian systems. For example, a 
cell-based display can be used to display large cDNA libraries in nfiammalian cells as descrit)ed in 
Nolan, et a!., U.S. Patent No. 6,153,380; Shioda , et al. U.S. Patent No. 6,251,676, both of which are 
expressly incorporated herein by reference. 

Screening of Libraries 
High-throughput screening technology 

Fully robotic or microfluidic systems include automated liquid-, particle-, celK and organism-handling 
including high throughput pipetting to perform all steps of experimental library generation, protein 
expression, and library screening. This Includes liquid, particle, cell, and organism manipulatbns such - 
as aspiration, dispensing, mixing, diluting, washing, accurate volumetric transfers; retrieving, and 
discarding of pipette tips; and repetitive pipetting of identical volumes for multiple deliveries from a 
single sample aspiration. These manipulations are cross-contamination-free liquid, particle, cell, and 
organism transfers. This instrument performs automated replication of microplate samples to filters, 
membranes, and/or daughter plates, high-density transfers, full-plate serial dilutions, and high 
capacity operation. 

In addition, as will also be appreciated by those in the art, biochlps may be part of the HTS system 
utilizing any number of components such as biosensor chips with protein arrays to measure protein- 
protein interactions or DNA-sensor chips to measure protein-DNA interactions. Microfluidic chip 
arrays (e.g., those commercially available from Caliper) may also be utilized in the context of 
automated HTS screening. 

The automated HTS system used can include a computer workstation comprising a microprocessor 
programmed to manipulate a device selected from the group consisting of a thermocycler, a 
multichannel pipetter, a sainpie handler, a plate handler, a gel loading system, an automated 
transformation system, a gene sequencer, a colony picker, a bead picker, a cell sorter, an incubator, a 
liglit microscope, a fluorescence microscope, a spectrofluorimeter, a spectrophotometer, a 
luminometer, a CCD camera and combinations thereof. 

in vivo screening 

In a prefenred embodiment, the library is screened using in vivo assay systems, including cell-based, 
tissue-based, or whole-organism assay systems. Cells, tissues, or organisms may be exposed to 
individual library members or pools containing several library members. Alternatively, host cells can 
be transformed or transfected with DNA encoding the library proteins and analyzed for phenotypic 
alterations. 
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To screen the library, experimental systems are developed in which the activity for the library protein 
of interest is coupled to an observable property. Typical observable properties include changes in 
absorbance. fluorescence, or luminescence. Screens may also monitor changes in properties such 
as cell morphology or viability. 

For example, cen death or viability can be measured using dyes or immuno^ochemical reagents 
(e.g. Caspase staining assay for apoptosis. Alamar blue for cell vitality) that specificalV recognize 
either viable or inviable cells. 

In an alternate cell death or viabflity assay, the cells are transformed or transfected writh a receptor or 
binding partner protein responsive to the ligand represented by the library. The receptor may be 
coupled toa signaling pathway that causes cell death, allows cell survival, or triggers expression of a 
reporter gene. These readout modalities can be measured using dyes or immuno-cytochemical 
reagentsthat Indicate cell death, cell vitality (e.g. Caspase staining a ssay for apoptosis. Alamar blue 
for cell vitality). 

Alternatively, readout can be via a reporter construct Reporter constructs may be proteins that are 
intrinsically fluorescent or colored, or proteins that modify the spectral properties of a substrate or 
binding partner. Common reporter constructs include the proteins luclferase. green fluorescent ' 
protein, and beta-gabctosldase. 

The assays described can also be perfomied by measuring morphological changes of the cells as a 
response to the presence of a library variant These morphotogical changes can be registered using 
microscopic image analysis systems (e.g. Cellomics ArrayScan technology) such as those now 
available commercially. 

in vitro screening 

In a preferred embodiment, different physical and functional properties of the library members are 
screened in an in vitro assay. Properties of library members that may be screened include, but are 
not limited to. various aspects of stability (including pH. thermal, oxidative/reductive and solvent 
stebility). solubility, affinity, activity and specificity. Multiple properties can be screened 
simultaneously (e.g. substrate spedfidty in organic solvents, receptor-ligand binding at tow pH) or 
individually. 

Protein properties can be assayed and detected in a wide variety of ways. Typical readouts include- 
but are not limited to, chromogenic. fluorescent, luminescent, or isotopic signals. These detection ' 
modalities are utBized in several assay methods including, but not limited to. FRET (fluorescence 
resonance energy transfer) and BRET(bioluminescence resonance energy transfer) based assays. 
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AlphaScreen (Amplified Luminescent Proximity Homogeneous Assay). SPA (scintillation proximity 
assay), ELISA (enzynr^inked immunosorbent assays), BIACORE (surface plasmon resonance), or 
enzymatic assays. !n vitro screening may or may not utilize a protein fusion or a label. 

Selection of Libraries 

In an alternatively prefen-ed embodiment, a selection method is used to select for desired library 
members. This is generally done on the basis of desired phenotypic properties, e.g. the protein 
properties defined herein. This is enabled by any nnethod which couples phenotype and genotype, 
l.e. protein function with the nucleic add that codes for rt. In some cases this will be a "trans" effect 
rather than a "cis" effecL In this way, isolation of library protein variants simultaneously enables 
isolation of its coding nucleic acid. Once isolated, the gene or genes encoding library protein can be 
purified ("rescued") and/or amplified. This process of isolation and amplification can be repeated, 
allowing favorable protein variants in the library to be enriched. Nucleic acid sequencing of the 
selected library members ultimately allows for identification of library members witii desired 
properties. 

Isolation of library protein can be accomplished by a number of metiiods. In some embodiments, only 
cells containing library protein variants with desired protein properties are allowed to survive or 
replicate. In alternate embodiments, the library protein and its genetic material are obtained by 
binding the library protein to anotiier protein, RN A aptamer, or otiier nrK>lecuIe. 

In one embodiment, the selection method is based on the use of specific fusion constmcts. For 
example, if phage display is used, the library members are fused to the phage gene 111 protein. 

In one embodiment selection is accomplished using a rescue fusion sequence, which forms a 
covalent or noncovalent link between the library member (phenotype) and the nucleic ackJ that 
encodes the library member (genotype). For example, in a preferred embodiment the rescue fusion 
protein binds to a specific sequence on the expression vector (see U. S.S.N. 09/642,574; 
PCT/USOO/22906; U. S.S.N. 10/023,208; PCT/U S0 1/49058; U.S.S.N. 09/792,630; U.S.S.N. 
10/080,376; PCT/US02/04852; U.S.S.N. 09/792,626; PCT/US02/04853; U.S.S.N. 10/082.671; 
U.S.S.N. 09/953,351; PCTUS01/28702; U.S.S.N. 10/097,100; and PCTyUS02/07466), and envelope 
virus (see U.S.S.N. 09/922,503 and PCT/US01^4535). 

In an alternate embodiment, selection is accomplished using a display technology including, but not 
limited to phage display, In which the library members are fused to a protein such as the phage gene 
Hi protein, (see Kay, BK et al, eds. Phage display of peptides and proteins: a laboratory manual 
(Academic Press, San Diego, CA, 1996); Lowman HB, Bass SH, Simpson N, Wells JA (1991) 
Selecting high-affinity binding proteins by monovalent phage display. Bioechemistry 30:10832-10838; 



85 



wo 03/014325 PCT/US02/25588 



Smith GP Filamentous fusion phage: novel expression vectors that display ctaned ant^ens on the 
virion surface. (1985) Science 228:1315-1317.) and its derivatives such as selective phage Infection 
(see Malmborg AC. Soderiind E, Frost L. Borrebaeck CASelective phage infection mediated by 
epitope expression on F pilus. (1997) J Mol Biol 273:544-551.). selectively infective phage (see 
Krebber C, Spada S. Desplanq D, Krebber A Ge L, Pluckthun A Selectively infective phage (SIP): a 
mechanistic dissection of a novel method to select for protein-ligand interactions. (1 997) J Mol Biol 
268:61 9-630.), and delayed infecth/lty panning (see Benhar I, Azriel R, Nahary L. Shaky S. 
Berdichevsky Y. Tamarkin A. Wels W(2000) Highly efficient selection of phage antibodies mediated 
by display of antigen as Ijjp-OmpA' fusfons on live bacteria. J Mol Biol 301:893-904.), Other display 
technotogies. which could be used, include but are not limited to cell surface display (see Witrnip KD, 
Protein engineering by cell-surface display. Cun-. Opin. Biotechnology 2001, 12:395-399) such as 
display on bacteria (see Georgiou G, Poetschke HL, Stathopoutos C. Francisco JA. Practical 
applications of engineering gram-negative bacterial cell surfaces. Trends Biotechnol. 1993 
Jan;11(1):6-10; Georgiou G, Stathopoulos C. Daugherty PS. Nayak AR. Iverson BL. and Curtiss RR 
(1997) Display of heterologous proteftis on the surface of mkafoorganisms: from the screening of 
combinatorial Hbraries to live reconibinant vaccines. Nahire Btotechnol. 15, 29-34; Lee JS. Shin KS, 
Pan JG, Kim CJ. Suriiace-displayed viral antigens on Salmonella canier vaccine. Nature 
Biotechnology, 2000, 18:645-648; Jun HC, Lebeault JM. Pan JG. Surface display of Zymomonas 
mobilislevansucmse by using the ice-nucleation protein of Pseudomonas syringae. Nat Biotechnol 
1998, 16:576-80.), yeast (see Boder ET and Wittrup KD. Yeast surface display for directed evolution 
of protein expression, affinity, and sfability. Methods Enzymol 2000. 328:430-44.; Boder ET. Wittrup 
KD: Yeast surface display for screening combinatorial polypeptide libraries. Nat Biotechnol 1997, 
15:553-557), and mammalfan cells (see Whitehom EA, Tate E, Yanofeky SD, Kochersperger L, Davis 
A, Mortensen RB, Yonkovich S. Bel K. Dower WJ. and Barrett RW (1995). A generic method for 
expression and use of "tagged" soluble versions of cell surface receptors. Bfci/technotogy, 13, 1215- 
1219.), as weH as in vitro display technologies such as polysome display (see Mattheakis LC, Bhatt 
RR, Dower WJ. Proc. Natl Acad Sci USA 1994, 91: 9022-9026; Hanes J and Pluckthun A Proc Natl 
Acad Sd USA 1997, 94:4937-4942.), ribosome display (see Hanes J and Pluckthun A Proc Natl Acad 
Sci USA 1997, 94:4937-4942). mRNAdispfay (Roberts RWand Szostak JW Proc Natl Acad Sci USA 
1997, 94, 12297-12302; Nemoto N, Miyamoto-Sato E, Husimi Y, Yanagawa H FEBS Lett 1997. 
414:405-408), and ribosome-inactlvatlon display system (see Zhou J. Fujifa S. Vterashina M. Baba. T, 
Taira K J Am Chem Soc (2002). 124, 538-543.) 

In an altemate embodiment, in vitro selection methods that do not rely on display technologies are 
used. These methods include but are not limited to periplasmic expression and cytometric screening 
(see Chen G. Hayhurst A. Thomas JG, Han/ey BR, Iverson BL, Georgiou G: Isolation of high-affinity 
ligand-binding proteins by peripfasmic expression with cytometric screening (PECS). Nat Biotechnol 
2001, 19: 537-542), protein fragment complementation assay (see Johnsson N & Varshavsl^ A Split 
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Ubiquitin as a sensor of protein interactions In vivo. (1 994) Proc Natl Acad Sci USA, 91:1 0340- 
10344.) and the yeast two hybrid screen (see Fields S, Song O: A novel genetic system to detect 
protein-protein interactions. Nature 1989, 340:245-246.) used in selection mode (see VIsintin M, Tse 
E, Axelson H, Rabbltts TH, Cattaneo A: Selection of antibodies for intracellular function using a two- 
hybrid in vivo system. Proc Natl Acad Sci USA 1999. 96: 11723-11728.). 

In an alternative embodiment, in vh/o selection can occur if expression of the library protein imparts 
some growth, reproduction, or survival advantage to the cell. For example, if host cells transformed 
with a library comprising variants of an essential enzyme are grown in the presence of the 
con-esponding substrate; only clones with a functional variant of the enzyme will survive. 
Altematively, an advantage may be conferred if the library member comprises a growth or survival 
factor and the host cell expresses the appropriate receptor. 

Additional Characterization 

In a preferred embodiment, a library member or members isolated using some screening or selection 
method are further characterized. The library member(s) may be subjected to further biological, 
physical, structural, kinetic, and thermodynamic analysis. Thus, for example, a selected library 
variant may be subjected to physical-chemical characterization using gel electrophoresis, reversed- 
phase HPLC, SEC-HPLC, mass spectrometry (MS) including but not limited to LC-MS, LC-MS 
peptide mapping and the like, ultraviolet absorbance spectroscopy, fluorescence spectroscopy, 
drcular dichroism spectroscopy, isothermal titration calorimetry, differential scanning calorimetry, 
surfece plasmon resonance, analytical ultra-centrifugation, proteolysis, and cross-linking. Structural 
analysis employing X-ray crystailographic techniques and nuclear magnetic resonance spectroscopy 
are also useful. As is known to those skilled in the art. several of the above methods can also be 
used to detemnine the. kinetics and thermodynamics of binding and enzymatic reactions. The 
biological properties of one or more library members, including pharmacokinetics and toxicity, can 
also be characterized in cell, tissue, and whole organism experiments. 

Expression vectors 

Using the nuclec acids of the present invention, which encode library members, a variety of 
expression vectors are made. The expression vectors may be either self-replicating 
extrachromosomal vectors or vectors which integrate into a host genome. 

Nucleic acki is operably linked when it is placed into a functional relationship with another nucleic acid 
sequence. For example, DNA for a presequence or secretory leader is operably linked to DNA for a 
polypeptide if it is expressed as a preprotein that participates in the secretion of the polypeptide; a 
promoter or enhancer is operably linked to a coding sequence if it affects the transcription of the 
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sequence; or a ribosome binding site is operably finked to a coding sequence if it is positioned so as 
to fadOtate translation However, enhancers do not have to be contiguous. 



Inclusion of control or regulatory sequences 

Generally, these expression vectors Include transcriptional and translational regulatory nucleic acid 
operably linked to the nucleic acid encoding the library protein. 

The transcriptional and translational regulatory nucleic acid will generally be appropriate to the host 
ceB used to express tiie library protein, as will be appredated by those in the art; for example, 
transcriptional and translational regulatory nudeic acid sequences from Bacillus are preferably used 
to express the library protein in Bacfflfos. Numerous types of appropriate expression vedors. and 
suitable regutetory sequences are known in the art for a variety of host ceBs. 

In general, the transcriptional and translational regulatory sequences may include, but are not limited 
to. promoter sequences, ribosomal binding sites, transcriptional start and stop sequences, 
translational start and slop sequences, and enhancer or activator sequences. In a prefen^ 
embodiment, the regulatory sequences indude a promoter and transcriptional start and stop 
sequences. 

Promoter sequences Include constitutive and inducible promoter sequences. TTie promoters may be 
naturally occumng promoters, hybrid or synttietic promoters. Hybrid promoters, whfch combine 
elements of more than one promoter, are also known in the art. and are useful in the present 
invention. 



Inclusion of a selectable marl(er(s) 

In addition, in a preferred embodiment the expression vector contains one or more selectable genes 
or parts of selectable mariner genes to allow ttie selectton of transfonried host cells containing the 
expression vector, and particularty in the case of mammalian cells, ensures ttie stability of the vector 
since cells which do not contain tiie vector wfll generally die. Selection genes are well known in the ' 
art and will vary with ttie host cell used. 

The bacterial expression vector may also include at least one selectable mari<er gene(s) to allow for 
ttie selection of bacterial strains that have been transformed. Suitable selectable gene(s) or parts of 
selectable martcer genes, indude genes, which render ttie baderia resistant to drugs such as 
ampidllin. ditoramphenfcol. erytiiromydn. kanamycin. neomydn and tetracydine. Selectable markers 
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also include blosynthetic genes, such as those In the histldine. tryptophan and leucine btosynthetic 
pathways. 

Inclusion of additional elements 

In a preferred embodiment, the expression vector contains a RNA splicing sequence upstream or 
downstream of the gene to be expressed in order to increase the level of gene expression. See 
Barret etal., Nucleic Acids Res. 1991; Groosetal., Mol. Cell, Biol. 1987; and Budiman etal., Mol. 
Cell. Biol. 1988. 

In addition, the expression vector way comprise additional elements. For example, the expression 
vector may have two replication systems, thus allowing it to be maintained in two organisms, for 
example in mammalian or insect cells for expression and In a prokaryotic host for cloning and 
amplification. Furthermore, for integrating expression vectors, the expression vector contains at least 
one sequence homologous to the host cell genome, and preferably two homologous sequences which 
flank the expression construct. The integrating vector may be directed to a specific locus in the host 
cell by selecting the appropriate homologous sequence for inclusion in the vector. Such vectors may 
include crs-lox recombination sites, or attR. affi, aftP. and aftL sites. Constructs for integrating 
vectors and appropriate selection and screening protocols are well known in the art and are described 
in e.g., Mansour et al.. Cell, 51:503 (1988) and Murray, Gene Transfer and Expression Protocols, 
Methods in Molecular Biology, Vol 7 (Clifton: Humana Press, 1991). 

Constructs 

Targeting or signal sequences 

The expression vector may also include a signal peptide sequence that provides for secretion of the 
library protein in bacteria. The signal sequence typically encodes a signal peptide comprised of 
hydrophobic amino acids which direct the secretion of the protein firom the cell, as is well known In the 
art. The protein is either secreted into the growth media (gram-positive bacteria) or into the 
periplasmic space, tocated between the inner and outer membrane of the cell (gram-negative 
bacteria). 

Thus, suitable targeting sequences Include, but are not limited to, binding sequences capable of 
causing binding of the expression product to a predetermined molecule or class of molecules while 
retaining bioactivity of the expression product, (for example by using enzyme inhibitor or substrate 
sequences to target a class of relevant enzymes); sequences signaling selective degradation, of itself 
or co-bound proteins; and signal sequences capable of constitutlvely localizing the candidate 
expressfon products to a predetermined cellular locale, including a) subcellular locations such as the 
Golgi. endoplasmic reticulum, nucleus, nucleoli, nuclear membrane, mitochondria, chbroplast, 
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secretory vesicles, lysosome, and cellular membrane; and b) extracellular locations via a secretory 
signal. Particularly preferred is localization to either subcellular locafions or to the outside of the ceH 
via secretion. 



iD(Purificatton>Taos 

In a prefenned embodiment, the library member comprises a rescue sequence operably linked to the 
rest of the peptide or protein. A rescue sequence is a sequence which may be used to purify or 
isolate either the candidate agent or the nucleic acid encoding it Thus, for example, peptide rescue 
sequences include purification sequences such as polyhistidlnes, including but not limited to the Hise. 
and the like or other tag for use with Ni*^ affinity columns and epitope tags for detection, 
immunopredpitation or FACS (fluorescence-activated cell sorting). Suitable epitope tags include c- 
myc (for use witti the commercially available 9E10 antibody), ttie BSP biotinylatbn target sequence of 
the bacterial enzyme BirA, flu tags. lacZ, and GST. 

A rescue sequence could also be a nucleic acid sequence operably linked to an epitope in a 
covalently attached protein, or a protein that specifically recognizes the nucleic acid. Such sequences 
include, but are not limited to. most sequence specific RNA and DNA binding proteins, preferably 
those ttiat recognize specific sequences or structures, and the like. 

Alternatively, the rescue sequence may be a unique oUgonudeotade sequence tiiat serves as a probe 
target site to allow the quick and easy isolation of the construct, via PCR. related technkjues. or 
hybridization. 

In a preferred embodiment, rescue sequences could also be based upon in vivo recombination 
systems, such as the cre-tox system, ttie Invifr^ogen Gateway™ system, forced recombination 
systems in yeast mammafian, plant bacteria or fungal cells (for example WO 02/10183 A1). or phage 
display systems. 



Fusion constructs 

The library protein may also be made as a fosion protein, using techniques well known in tiie art 
Thus, for example, fertile creation of monoctonal antibodies, if ttie desired epitope is small, tiie library 
protein may be fosed to a earner protein to form an immunogen. Altematively, ttie Bbrary protein may 
be made as a fusion protein to increase expression, or for other reasons. For example, when ttie 
library protein is a library peptide, ttie nucleic acid encoding ttie peptide may be linked to ottier nucleic 
acid for expression purposes. Similarty, ottier fusion partners may be used, such as targeting 
sequences whrch albw ttie kjcalization of the flbrary members into a subcellular or extiacellular 
compartonent of the cell, rescue sequences or purification tegs which altow the purification or isolation 
of eittier the library protein or ttie nuciek: ackls encoding ttiem; stability sequences, which confer 
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stability or protection from degradation to the library protein or the nucleic acid encoding it, for 
example resistance to proteolytic degradation, or combinations of these, as well as Wnker sequences 
as needed. 

In a preferred embodiment, the fusion partner is a stability sequence to confer stability to the library 
member or the nucleic acid encoding il Thus, for example, peptides may be stabilized by the 
incorporation of glycines after the initiation methionine (MG or MGGO), for protection of the peptide to 
ubiqultination as per Varshavsky's N-End Rule, thus confemng long half-life in the cytoplasm. 
Similarly, two prolines at the C- terminus impart peptides that are largely resistant to carboxypeptidase 
action. The presence of two glycines prior to the prolines impart both flexibility and prevent structure 
initiating events in the di-proline to be propagated into the candidate peptide structure. Thus, 
preferred stability sequences are as follows: MG(X)„GGPP, where X is any amino acid and n is an 
Integer of at feast four. 

Labeling (isotopic, fluorescent affinity) 

In one embodiment, the library nucleic acids, proteins and antibodies of the invention are labeled. By 
"labeled" herein is meant that nucleic acids, proteins and antibodies of the invention have at least one 
element, isotope or chemical compound attached to enable the detection of nucleic acids, proteins 
and antibodies of the invention. In general, labels fall Into three classes: a) isotopic labels, which may 
be radioactive or heavy isotopes; b) affinity labels, which may be antibodies or antigens; and c) 
colored or fluorescent dyes. The labels may be Incorporated into the compound at any position. 

Expression systems 

The library proteins of the present invention are produced by culturing a host cell transformed with 
nucleic acid, preferably an expression vector, containing nucleic acid encoding an library protein, 
under the appropriate conditions to induce or cause expression of the library protein. As outlined 
below, the libraries may be the basis of a variety of display techniques, including, but not limited to, 
phage and other viral display technologies, yeast, bacterial, and mammalian display technologies. 
The conditions appropriate for library protein expression will vary virith the choice of the expression 
vector and the host cell, and will be easily ascertained by one skilled in the art through routine 
experimentation. For example, the use of constitutive promoters in the expresston vector will require 
optimizing the growth and proliferation of the host cell, while the use of an inducible promoter requires 
the appropriate growth conditions for induction. In addition, in some embodiments, the timing of the 
harvest is important. For example, the baculoviral systems used in insect cell expression are lytic 
viruses, and thus harvest time selection may be crucial for product yield. 
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As Will be appredated by those in the art, the type of cells used in the present invention may vary 
widely. Basically, a wide variety of appropriate host cells may be used, including yeast, bacteria, 
archaebacteria, fungi, and insect and animal cells, including mammalian cells. Of particular interest 
are Drosophila melanogaster ceWs, Saccharomyces cerevisiae and other yeasts, E coh\ Bacillus 
subtilis, SF9 ceHs, C129 cells, 293 cells, Neurospora, BHK, CHO, COS, and HeLa cells, fibroblasts, 
Schwanoma cell lines, immortalized mammalian myeloid and lymphoid cell lines. Juri<at cells, mast 
cells and other endocrine and exocrine cells, and neuronal cells. See the ATCC cell line catalog, 
hereby expressly Incorporated by reference. In addition, the expression of the secondary libraries in 
phage display systems, such as are well known in the art. are particularty preferred, especially when 
the secondary library comprises random peptides. In one embodinnent, the cells may be genefically 
engineered, that is, contain exogenous nudeic acid, for example, to contain target molecules. 

Mammalian expression systems 

In a preferred embodiment, the library proteins are expressed In mammalian cells. Any mammalian 
ceils may be used, with nwuse. rat, primate and human ceils being particularly prefen-ed. although as 
will be appreciated by those in the art, modifications of the system by pseudotyping allows all 
eukaryotic cells to be used, preferably higher eukaryotes. As is more fully described below, a screen 
will be set up such that the cells exhibit a selectable phenotype in the presence of a random library 
member. As is more fully described below, cell types implicated in a wide variety of disease 
conditions are particularly useful, so long as a suitable screen may be designed to allow the selection 
of cells that exhibit an altered phenotype as a consequence of the presence of a library member 
within the cell. 

Accordingly, suitable mammalian cell types include, but are not limited to, tumor cells of all types 
(particularly melanoma, myeloid leukemia, carcinomas of the lung, breast, ovaries, colon, kidney, 
prostate, pancreas and testes), cardiomyocytes, endothePial cells, epithelial cells, lymphocytes (T-cell 
and B cell) . mast cells, eosinophils, vascular intimal cells, hepatocytes, leukocytes including 
mononuclear leukocytes, stem cells such as haemopoietic, neural, skin, lung, kidney, liver and 
myocyte stem cells (for use in screening for differentiation and de-differentiation factors), osteoclasts, 
chondrocytes and other connective tissue cells, kerab'nocytes, melanocytes, liver cells, kidney cells, 
and adipocytes. Suitable cells also include known research cells, including, but not limited to, Juriot 
T cells, NIH3T3 cells. CHO, COS, eta See the ATCC cell line catalog, hereby expressly incorporated 
by reference. 
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Mammalian expression systems are also known In the art, and Include retroviral systems. A 
mammalian promoter is any DMA sequence capable of binding mammalian RNA polymerase and 
initiating the downstream (3*) transcription of a coding sequence for library protein into mRNA, A 
promoter will have a transcription-initiating region, which is usually placed proximal to the 5* end of the 
coding sequence, and a TATA box, usually located 26-30 base pairs upstream of the transcription 
initiation site. The TATA box is thought to direct RNA polymerase II to begin RNA synthesis at the 
correct site. A mammalian promoter will also contain an upstream promoter element (enhancer 
element), typically located within 100 to 200 base pairs upstream of the TATA box. An upstream 
pronrK)ter element determines the rate at which transcription Is initiated and may act in either 
orientation. Of particular use as mammalian pronx^ters are the promoters firom mammalian viral 
genes, since the viral genes are often highly expressed and have a broad host range. Examples 
include the SV40 early promoter, mouse mammary tumor virus LTR promoter, adenovirus major late 
promoter, herpes simplex virus promoter, and the CMV promoter. 

Typically, transcription tenmination and polyadenylation sequences recognized by mammalian cells 
are regulatory regions located 3' to the translation stop codon and thus, together with the promoter 
elements, flank the coding sequence. The 3' terminus of the mature mRNA is formed by site-specific 
post-translational cleavage and polyadenylation. Examples of transcription terminator and 
polyadenylation signals include those derived from SV40. 

The methods of Introducing exogenous nucleic acid into mammalian hosts, as well as other hosts, Is 
well known in the art, and will vary with the host cell used. Techniques include dextran- mediated 
transfection, calcium phosphate precipitation, polybrene mediated transfectlon, protoplast fusion, 
electroporation, viral infection, encapsulatbn of the polynuc!eotide(s) in liposomes, and direct 
microinjection of the DNA into nuclei. 

Bacterial expression systems 

In a preferred embodiment, library proteins are expressed in bacterial systems. Bacterial expression 
systems are well known In the art. 

A suitable bacterial promoter is any nucleic acid sequence capable of binding bacterial RNA 
polymerase and initiating the downstream (3') transcription of the coding sequence of library protein 
into mRNA A bacterial promoter has a transcription Initiation region which is usually placed proximal 
to the 5' end of the coding sequence. This transcription initiation region typically Includes an RNA 
polymerase binding site and a transcription initiation site. Sequences encoding metabolic pathway 
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enzymes provide particulariy useful promoter sequences. Examples include promoter sequences 
derived from sugar metabolizing enzymes, such as galactose, lactose and maltose, and sequences 
derived from biosynthetic enzymes such as tryptophan. Promoters from bacteriophage may also be 
used and are known in the art In addition, synthetic promoters and hybrid promoters are also useful; 
for example, the tac promoter is a hybrid of the trp and lac promoter sequences. Furthemiore. a 
bacterial promoter may include naturally occurring promoters of non-bacterial origin that have the 
ability to bind bacterial RN A polymerase and initiate transcrv)tion. 

In addition to a functioning promoter sequence, an effident ribosome-binding site is desirable. In £ 
coff. the ribosome-binding site is called the Shine-Dalgamo (SD) sequence and includes an initiation 
codon and a sequence 3-9 nucleotides In length located 3-11 nucleotides upstream of the initiation 
codon. 



Baculovirus expression system 

In one embodiment library proteins are produced in insect cells. Expression vectors for the 
transfonnation of insect cells, and in particular, bacutovirus-based expression vectors, are well knovm 
in the art and are described e.g., in O'Reilly et al., Bacutovinis Expression Vectors: A Laboratory 
Afawa/ (New York: Oxford University Press, 1994). 

Yeast expression systems 

In a prefen-ed embodiment, library protein is produced in yeast cells. Yeast expression systems are 
well known in the art, and Include expression vectors for Saccharomyces cerevisiae. Candida 
albicans and C. maltosa, Hansenula polymorpha, Kluyveromycesfragilisand K. lactis, Pichia 
guHlerimondiiand P. pastoris, Schizosaccharomyces pombe, and Yarrowia lipolytica. Preferred 
promoter sequences for expression In yeast include the inducible GAL1 ,10 promoter, the promoters 
fi^om alcohol dehydrogenase, enolase. glucokinase, glucose-6-phosphate isomerase, glyceraWehyde- 
3-phosphate-dehydrogenase. hexokinase, phosphofiructoklnase, 3-phosphogiycerate mufase. 
pyruvate kinase, and the acid phosphatase gene. Yeast selectable markers Include ADE2, HIS4, 
LEU2, TRP1 . and ALG7, which confers resistance to tunicamycin; the neomycin phosphotransferase 
gene, which confers resistance to G418: and the CUP1 gene, which allows yeast to grow in the 
presence of copper ions. 



In Vitro Expression systems 

In one embodiment, the library proteins are expressed in vitro using cell-free translation systems. 
Several commerdal sources are available for this system including but not limited to Roche Rapid 
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Translation System, Promega TnT system, Novagen's EcoPro system, Ambion's ProteinScipt-Pro 
system. In vitro translation systems derived from both prokaryotic (e.g. E. coli) and eulcaryot'c (e.g. 
Wheat germ. Rabbit reticuiocytes) cells are available and may be chosen based on the expression 
levels and functional properties of the protein of interest. Both linear (as derived from a PGR 
amplification) and circular (as in plasmid) DNA molecules are suitable for such expression as long as 
they contain tlie gene encoding the protein operably linked to an appropriate promoter. Other 
features of the molecule that are important for optimal expression in either the bacterial or eukaryotic 
cells (including the ribosome binding site etc) are also included in these constructs. The proteins may 
again be expressed individually or in suitable size pools consisting of multiple library members. The 
main advantage offered by these in vitro systems is their speed and ability to produce soluble 
proteins. In addition the protein being synthesized may be selectively labeled if needed for 
subsequent functional analysis. 

Protein purification 

In a preferred embodiment, the library protein Is purified or isolated after expression. Library proteins 
may be isolated or purified in a variety of ways known to those skilled In the art depending on what 
other components are present in the sample. Standard purification methods include electrophoretic, 
molecular, Immunological and chromatographic techniques, including bn exchange, hydrophobic, 
affinity, and reverse-phase HPLC chromatography, and chromalofocuslng. For example, the library 
protein may be purified using a standard anti-library antibody column. Ultrafiltration and diafiltration 
techniques, in conjunction with protein concentration, are also useful. For genera! guidance in 
suitable purification techniques, see Scopes, R., Protein Purification, Springer- Verlag, NY (1982). 
The degree of purification necessary will vary depending on the use of the libraiy protein. In some 
instances no purification will be necessary. 

Screening of library members 

Library members may be screened using a variety of assays, including but not limited to in vitro 
assays, and in vivo assays such as cell-based, tissue-based, and whole-organism assays. 
Automation and high-throughput screening technologies may be utilized in the screening procedures. 

Ceil'based assays - eukaryotic and prokaryotic 

In a preferred embodiment, the library is screened using cell-based assay systems. 
in vivo selection of library variants 

Host cells transformed with a library representing variants of an enzyme or resistance factor of 
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interest are grown in the presence of the corresponding substrate or antibiotic. Only dones with a 
functional variant of the enzyme or resistance factor will survive. 



Screening based on cell survival, ceti death or expression of reporter genes in cells 

Cells are exposed to individual variants or pools of variants belonging to a library to be assayed. The 
cells are transformed or transfected either transienUy or stably witii the corresponding receptor 
responsive to the figand represented by the library. The receptor is coupled to a signaling pathway 
that either causes cell death, cell sunflval. or triggers expression of a reporter gene. These readout 
modalities nnay be measured using dyes or Immuno-cytochemical reagents ttiat indicate cell death, 
cell vitality (e.g. Caspase staining assay for apoptosis, Alamar blue for cell vitality), or In case of the 
reporter constructs enzymes that convert dyes and cause them to be luminescent (e.g. lucHerase) or 
shift their absorbance or fluorescent properties to wavelengths different from their properties before 
conversion. 



Screening based on cell survival of indhr iduaf clones or clone pools 

Host cells are transformed or transfected with library DNA representing variants of a ligand or 
receptor of interest. The cells are also transformed or transfected either ti^nsientiy or stably witti the 
conresponding receptor responsive to tiie ligand represented by the library or in case of a receptor 
library with ligand signaling through the receptor represented by the library. The receptor is coupled 
to a signaling pathway tiiat causes cell survival. If the sequence of the variant causing cell survival is 
not pre-identified, surviving cell clones may be used to identify the sequence Identity of the 
corresponding variant. 



Screenin g based morphological changes of cells 

All of the above described assay readouts rely on changes that may be measured using absorbance, 
fluorescence or luminescence readers. The assays described may also be read measuring 
morphological changes of the cells as a response to the presence of a library variant These 
morphological changes may be registered using microscopic Image analysis systems (e.g. Cellomics 
ArrayScan technology) now available commercially. 



Screening based on candidate bioactive agents 

Candidate agents are obtained from a wide variety of sources, as will be appreciated by those In tiie 
art, including libraries of syntfietic or natural compounds. As will be appreciated by tiiose in the art. 



96 



wo 03/014325 



PCT/US02/25588 



the present invention provides a rapid and easy method for screening any library of candidate agents, 
Including the wide variety of known combinatorial chemistry -type libraries. 

In a preferred embodiment, candidate agents are synthetic compounds. Any number of techniques 
are available for the random and directed synthesis of a wide variety of organic compounds and 
biomolecules, including expression of randomized oligonucleotides. See for example WO 94/24314, 
hereby expressly incorporated by reference, which discusses methods for generating new 
compounds, including random chemistry methods as well as enzymatic methods. As described in 
WO 94/24314, one of the advantages of the present method is that it is not necessary to characterize 
the candidate bioactive agents prior to the assay; only candidate agents that bind to the target need 
be identified. In addition, as is known in the art coding tags using split synthesis reactions may be 
done, to essentially Identify the chemical moieties on the beads. 

Alternatively, a preferred embodiment utilizes libraries of natural compounds in the form of bacterial, 
fungai. plant and aninnal extracts that are available or readily produced, and can be attached to beads 
as is generally Icnown in the art 

Additionally, natural or synthetically produced libraries and compounds are readily modified through 
conventional chemical, physical and biochemical means. Known pharmacological agents may be 
subjected to directed or random chemical modlficafions, including enzymatic modifications, to produce 
structural analogs. 

In a preferred embodiment, candidate bioactive agents include proteins, nucleic acids, and chemical 
moieties. 

In a preferred embodiment, the candidate bioactive agents are proteins. In a preferred embodiment, 
the candidate bioactive agents are naturally occurring proteins or fragments of naturally occurring 
proteins. Thus, for example, cellular extracts containing proteins, or random or directed digests of 
proteinaceous cellular extracts, may be attached to beads as is more fully described below. In this 
way libraries of procaryotic and eucaryotic proteins may be made for screening against any number of 
targets. Particularly preferred in this embodiment are libraries of bacterial, fungal, viral, and 
mammalian proteins, with the latter being preferred, and human proteins being especially pretended. 

In a preferred embodiment, the candidate bioactive agents are peptides of fi^om about 2 to about 50 
amino acids, with from about 5 to about 30 amino adds being preferred, and from about 8 to about 20 
being particularly preferred. The peptides may be digests of naturally occurring proteins as is outlined 
above, random peptides, or "biased" random peptides. By"randomlzed" or grammatical equivalents 
herein is meant that each nucleic acid and peptide consists of essentially random nucleotides and 
amino acids, respectively. Since generally these random peptides (or nucleic acids, discussed below) 
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are chemically synthesized, they may incorporate any nucleotide or amino acid at any position. The 
synthetic process can be designed to generate randomized proteins or nucleic adds, to allow the 
tbrmation of ail or most of the possible combinations over the length of the sequence, thus fbnning a 
library of randomized candidate bioactive proteinaceous agents. In addition, the candidate agents 
may themselves be the product of the invention; that is. a library of proteinaceous candidate agents 
may be made using the methods of the invention. 

Hlqh4hro uqhput screening technology 

Fully robotic or microfluidic systems include automated liquid- , particle-, cell- and organism-handling 
including high throughput pIpetBng to perfomn aD steps of gene targeting and recombination 
applications. This includes liquid, particle, cell, and organism manipulations such as aspiration, 
dispensing, mixing, diluting, washing, accurate volumetric transfers; retrieving, and discarding of 
pipette tips; and repetitive pipetting of Identical volumes tor mulbple deliveries from a single sample 
aspiration. These manipulations are cross-contamination-firee liquid, particle, cell, and organism 
transfers. This instrument performs automated replication of micrxjplate samples to filters, 
membranes, and/or daughter plates, high-density transfers, full-plate serial dilutions, and high 
capacity operation. 



In addition, as will also be appreciated by those in the art biochips may be part of the HTS system 
utifizing any number of components such as biosensor chips with protein arrays to measure pix)tein- 
protein interactions or DislA^ensor chips to measure protein-DNA interactions. Microfluidic chip 
anays (e.g.. technology developed by Caliper) may also be utiBzed in the context of automated HTS 
screening. 



The automated HTS system used may include a computer workstation comprising a microprocessor 
programmed to manipulate a device selected from the group consisting of a thermocycler. a 
multichannel pipetter, a sample handler, a plate handler, a gel loading system, an automated 
transformation system, a gene sequencer, a colony picker, a bead picker, a cell sorter, an incubator, a 
light microscope, a fluorescence microscope, a spectrofluorimeter. a spectrophotometer, a 
luminometer, a CCD camera and combinations thereof. 



In vitro assay s 

In a prefened embodiment, different physical and functional properties of the library members are 
screened in an in vitro assay. In vitro assays aHow a broader dynamic range lor screening protein 
properties of interest that are not limited by cellular viability of the cells expressing ttie Hbrary 
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members or library members acting upon other cells to exert Its effects. Properties of library members 
that may be screened include, but are not limited to, various aspects of stabili^ (including pH, 
thermal, oxidative/reductive and solvent stability), solubility, affinity, activity and specificity. Multiple 
properties may be saeened simultaneously (e.g. substrate specificity in organic solvents, receptor- 
Itgand binding at low pH) or individually. 

Protein properties may be assayed and detected in a wide variety of ways. Modality of detection 
could include, but are not limited to, chromogenic, fluorescent, luminescent, or isotopic substrates for 
protein library members. Any of these detection modalities are utilized in several assay methods 
including, but not limited to, FRET (fluorescence resonance energy transfer) and BRET 
(bioluminescence resonance energy transfer) based assays, AlphaScreen (Amplified Luminescent 
Proximity Homogeneous Assay), SPA (scintillation proximity assay), ELISA (enzyme-linked 
immunosorbent assays), or enzymatic assays. 

Additional characterization 

In a preferred embodiment, a library member or members isolated from a cell positively selected for 
any number of protein properties by in-vivo or in-vitro screening methods well known to those in the 
art, are further characterized for said properties by aforementioned screens or other methods 
including physical, structural, kinetic, and thermodynamic analysis. Thus, for example, a selected 
library variant may be subjected to physical characterization through gel electrophoresis, reverse- 
phase HPLC, MS, LC-MS. RP-HPLC. SEC-HPLC. LC-MS peptide mapping, CD, analytical ultra- 
centrifugation, and proteolysis. Structural analysis employing X-ray crystallographic techniques, 
NMR, and cross-linking are also useful. In addition, thermodynamic and kinetic characterization of 
proteinaceous moieties are well known in the art. 

EXAMPLES 

The following examples serve to more fully describe the manner of using the above-described 
invention, as well as to set forth the best modes contemplated for carrying out various aspects of the 
invention. It is understood that these examples in no way serve to limit the true scope of this 
invention, but rather are presented for illustrative purposes. All references cited herein are 
incorporated by reference. 

Example 1 

Computational Prescrecning on P-lactamase TEM-1 
Experiments were performed on the p-lactamase gene TEM-1 . Brookhaven Protein Data Bank entry 
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1 BTL was used as the starting structure. All water molecules and the S04=^ group were removed and 
expDcit hydrogens were generated on the structure. The structure was then minimized for 50 steps 
without electrostatics using the conjugate gradient method and the Drekling II force field. These steps 
were performed using the BIOGRAF program commercially available from Molecular Simulations, 
Inc.. San Diego, CA This minimized structure served as the template for all the protein design 
cabulattons. 

Computational Screening 

Connputational screening of sequences was performed using PDA™ technology. A 4 A sphere was 
drawn around the heavy side chain atoms of the four catalytic residues (S70, K73, 8130, and E166) 
and aU amino adds having heavy side chain atoms within this distance cutoff were selected. This 
yielded the following 7 positions: F72. Y105, N132, N136, L169, N170, and K234. Two of these 
residues. N132 and K234, are highly conserved across several different ^-lactamases and were 
therefore not included in the design, leaving five variable residue positions (Fll. Y105, N136, L169. 
N170). These designed positions were allowed to change their identity to any of the 20 naturally 
occurring amino acids except proline, cysteine, and glycine (a total of 17 amino acids). Proline is 
usually not allowed since it is difficult to define appropriate rotamers for proline, cysteine is excluded 
to prevent formation of disulfide bonds, and glycine is excluded because of conformational flexibility. 

Additionally, a second set of residues within 5 A of the five residues selected for PDA™ technology 
design were floated (their amino acid identity was retained as wild type, but their conformation was 
allowed to change). The heavy side chain atoms were again used to determine which residues were 
within the cutoff. This yielded the following 28 positions: M68, M69, S70, T71, K73, V74, L76, VI 03. 
E104, S106, P107. 1127. M129, S130. A135. L139. L148. L162. R164. W165. E166. P167, D179, 
IVI21 1 , D214. V216. S235, 1247. The two prolines, P107 and P167. were excluded from the floated 
residues, as were positions M69. R164. and W165, since their crystal structures exhibit highly 
strained rotamers, leaving 23 floated residues from the second set Also. A248 was included instead 
of 1247. The conserved residues N1 32 and K234 from the first sphere (4 A) were also floated, 
resulting in a total of 25 floated residues. 

The potential functions and parameters used in the PDA™ technology calculations were as follows, 
■me van der Vteals scale factor was set to 0.9. and the etecfrostetic potential was calculated using a 
distance dependent dielectric of e=40 R. The well depth for the hydrogen bond potential was set to 8 
kcal/mol with a local and remote backbone scale factor of 0.25 and 1 .0 respectively. The solvation 
potential was only calculated for designed positions classified as core (F72, L169, M68, T71, V74, 
L76. 1127, A135, L139. L148. L162, M211 and A248). Type 2 solvation was used (Street and Mayo, 
1998). The non-polar exposure multiplication factor was set to 1.6, the non-polar burial energy vras 
set to 0.048 kcal/moVA^ and the polar hydrogen burial energy was set to 2.0 kcal/mol. 
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The Dead End Elimination (DEE) optimization nriethod (see reference) was used to find the lowest 
energy, ground state sequence. DEE cutoffe of 50 and 100 kcal/mol were used for singles and 
doubles energy calculations, respectively. 

Starting from the DEE ground state sequence, a iVIonte Carlo (MC) calculation was performed that 
generated a list of the 1 000 lowest energy sequences. The MC parameters were 1 DO annealing 
cycles with 1 ,000,000 steps per cycle. The non-productive cycle limit was set to 50. In the annealing 
schedule, the high and low temperatures were set to 5000 and 100 K respectively. 
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The following probabifity distribution was then calcuialed from the top 1000 sequences in the MC list 
(see Table 31 below). It shows the number of occunences of each of the amino acids selected for 
each position (the 5 variable residue positions and the 25 floated positions). 

Table 1: Monte Carlo analysis (amino acids and their number of occurrences at the designed 
positions resulting from the MC list of the 1000 lowest energy ranked sequences. 



POSITION / 


69 


M:1000 


70 


S:1000 


71 


T:1000 


72 


Y:591 


73 


KrIOOO 


74 


V:1000 


76 


L:1000 


103 


V:1000 


104 


E:1000 


105 


M:183 


106 


S:1000 


127 


1:1000 


129 


M:1000 


130 


3:1000 


132 


N:1000 


135 


A:1000 


136 


D:530 


139 


L:1000 


148 


L:1000, 


162 


LIOOO 


166 


E:1000 



F:365 V:35 E8 L:1 



V:68. £66 S:38 T:33 A:27 Q:6 
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169 


L:698 


E:156 


M:64 


S:37 


D:23 


A:21 


170 


M:249 


L:118 


E:113 


D:112 


T:90 


Q:87 




N:24 


F:21 


K:15 


Y:9 


H:9 


V:8 


179 


D:1000 












211 


M:1000 












214 


D:1000 












216 


V:1000 












234 


KIOOO 












235 


S:1000 












248 


A:1000 
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This probability distribution was tiien transformed into a rounded probability distribution (see Table 2). 
A 10% cutoff value was used to round at the designed positions and the wild type amino acids were 
forced to occur with a probabiKty of at least 10%. An E was found at position 169 15.6% of the time. 
However, since this position is adjacent to another designed position, 170. its closeness would have 
required a more complicated oligonucleotide library design; E was therefore not included for this 
position when generating the sequence Bbrary (only L was used). 



POSITION 



RESIDUE/PROBABILITY 



72 


105 


136 


169 


170 


Y50% 

F50% 

1 20% 

N 10% 

E10% 

S10% 

Y10% 


IW20% 
Q20% 
N10% 


D70% 
M20% 


L 100% 

E20% 
D20% 
N10% 


M30% 
L20% 



As seen from Table 2, the computational pre-screening resulted in an enomious reduction in the size 
of the problem. Originally, 17 different amino adds were allowed at each of the 5 designed positions, 
giving 17= = 1.419.857 possible sequences. This was pared down to just 210 possible sequences — 
a reduction of neariy four orders of magnftude. 



Generation of Sequence Library 

Overlapping oligonucleotides corresponding to tiie full lengttj TEM-1 gene for ^-lactamase and all 
desired mutations were synthesized and used in a PGR reaction as described previbusly (Figure 1). 
resulting in a sequence library containing ttie 210 sequences described above. 
Synthesis of mutant TEM-1 genes 

To allow ttie mutation of the TEM-1 gene. pCR2.1 (commercially available from Invitrogen) was 
digested with Xbal and EcoRI. blunt ended witti T4 DNA polymerase, and religated. This removes the 
Hindlll and Xhol sites within the polylinker. A new Xhol site was then introduced into the TEM-1 gene 
at position 2269 (numbering as of the original pCR2.1) using a Quiclcchange Site-Directed 
Mutagenesis Kit as described by the manufacturer (commercially available from Stratagene). 
Similariy. a new Hindlll site was introduced at position 2674 to give pCR-Xenl. 
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To construct the mutated TEM-I genes, overlapping 40-mer oligonucleotides were synthesized 
conresponding to the sequence between the newly introduced Xhol and Hindlll sites, designed to 
allow a 20-nucleotide overlap with adjacent oligonucleotides. At each of the designed positions (72, 
105. 136 and 170) multiple oligonucleotides were synthesized, each containing a different mutation so 
that all the possible combinations of mutant sequences (210) could be made in the desired 
proportions as shown in Table 3. For example, at position 72, two sets of oligonucleotides were 
synthesized, one containing an F at position 72. the other containing a Y. Each oligonucleotide was 
resuspended at a concentration of 1|jg4Jl, and equal molar concentrations of the oligonucleotides 
were pooled. 

At the redundant positions, each oligonucleotide was added at a concentration that reflected the 
probabilities in Table 3. For example, at position 72 equal amounts of flie two oligonucleotides were 
added to the pool, while at position 136. twice as much M-encoding oligonucleotide was added 
compared to the N-containing oligonucleotide, and seven times as much D-containing oligonucleotide 
was added compared to the N-containing oligonucleotide. 

DNA library assemblv 

For the first round of PGR, 2 pi of pooled oligonucleotides at the desired probabilities (Table 3) were 
added to a 100 pi reaction that contained 2 pMO mM dNTPs. 10 pi 1 0x Taq buffer (commercially 
available from Qiagen), 1 pi of Taq DNA polymerase (5 units/pl: commercially available from Qiagen) 
and 2 pi Pfu DNA polymerase (2.5 units/pl: commercially available from Promega). The reaction 
mixture was assembled on ice and subjected to 94^0 for 5 minutes. 15 cycles of 94''C for 30 seconds, 
52'C for 30 seconds and 72°C for 30 seconds, and a final extension step of 72"C for 10 minutes. 

Isolation of full-lenoth oligonucleotides 

For the second round of PGR. 2.5 pi of the first round reaction was added to a 100 pi reaction 
containing 2 pi 10 mM dNTPs, 10 pi of 10x Pfu DNA polymerase buffer commercially available from 
Promega, 2 pi Pfu DNA polymerase (2.5 units/pl: commercially available from Promega), and 1 pg of 
oligonucleotides corresponding to the 5' and 3' ends of the synthesized gene. The reaction mixture 
was assembled on ice and subjected to 94°G for 5 minutes. 20 cycles of 94DG for 30 seconds, 52X 
for 30 seconds and 72°C for 30 seconds, and a final extension step of 72*'G for 10 minutes to Isolate 
the full length oligonucleotides. 

Purification of DNA library 

The PGR products were purified using a QIAquick PGR Purification Kit commercially available from 
Qiagen, digested with Xhol and Hindlll, electrophoresed through a 1.2 % agarose gel and re-purified 
using a QIAquick Gel Extraction Kit commercially available from Qiagen. 

Verification of Sequence Library Identitv 
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The PGR products containing the library of mutant TEM-1 p-lactamase genes were then cloned 
between a promoter and terminator in a kanamycin resistant plasmid and transformed into E colL An 
equal number of bacteria were then spread onto media containing either kanamycin or amplcilBn. All 
transformed colonies will be resistant to kanamycin. but only those with active mutated P-lactamase 
genes will grow on ampiciilin. After overnight Incubation, several colonies were observed on both 
plates, indicating that at least one of the above sequences encodes an active P-Jactamase. The 
number of colonies on the kanamycin plate far outnumbered those on the amprdllin plate (roughly a 
5:1 ratio) suggesting that either some of the sequences destroy activity, or that the PGR introduces 
en-ors that yield an inactive or truncated enzyme. 

To distinguish between these possibilities. 60 colonies were picked from the kanamycin plate and 
their plasmid DNA was sequenced. This gave the distribution shown in Table 3. 

Table 3: Percentages predicted by PDA™ technology vs. those observed from experiment for the 
designed positions. 



Wild Type 


PDA™ Technology Residues (Predicted Percentage/Observed Percentage) 


72F 


Y50/50 


F50/60 




105Y 


M20/27 


(320/18 


120/21 N10/7 E10/7 S10/10 Y10/10 


136N 


D70/72 


M20/17 


N10/11 


170N 


M 30/34 


L 20/21 


E 20/21 D 20/17 N 10/7 



Note that the observed percentages of each amino add at all four positions ctosely match the 
predicted percentages. Sequencing also revealed that only one of the 60 colonies contained a PGR 
error, a G to C transition. 

This small test demonstrates that multiple PGR with pooled ofigonucleotides may be used to construct 
a sequence library that reflects the desired proportions of amino acid changes. 

Experimental Screening of Sequence Library 

The purified PGR product containing the library of mutated sequences was then ligated into pCR- 
Xenl that had prevfously been digested with Xhol and Hindlll and purified. The ligation reaction was 
transformed into competent TOP10 £ coli cells (Invitrogen). After allowing the cells to recover for 1 
hour at 37X. the cells were spread onto LB plates containing the antibiotic cefotaxime at 
concentrations ranging from 0.1 pg/ml to 50 pg/ml and selected for increasing resistance. 

A triple mutant was found that improved enzyme function by 35-fold in only a single round of 
screening (see Figure 4). Thfe mutant (Y105Q. N136D. N170L) sun/ived at 50 pg/ml cefotaxime. 
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Example 2 

Secondary Library generation of a Xylanase 

PDA™ technology Screening Leads to Enormous Reducfon in Number of Possible Sequences 

To demonstrate that computational screening is feasible and will lead to a significant reduction in the 
number of sequences that have to be experimentally screened, calculations for the B, drculans 
xylanase with and without the substrate were performed. The PDB structures 1XNB of free B. 
drculans xylanase and 1 BCX for the enzyme substrate complex were used. 27 residues inside the 
binding site were visually identified as belonging to the active site. 8 of these residues were regarded 
as absolutely essential for the enzymatic activity. These positions were floated (see Figure 2). This 
means that they couki change their side chain conformation but not their amino acid identity. 

Three of the 20 naturally occurring amino acids were not considered (cysteine, proline, and glycine). 
Therefore, 17 different amino acids were still possible at the remaining 19 positions; the problem 
yields 17^® = 2.4 xlO^^ different amino acid sequences. This number Is 10 orders of magnitude larger 
than what may be handled by state of the art directed evolution methods. Clearly these approaches 
cannot be used to screen the complete dimensionality of the problem and consider all sequences with 
multiple substitutions. Therefore PDA™ technology calculations were performed to reduce the 
sequence space. Starting from the PDA™ technology ground state a list of 10,000 low energy 
sequences was created by Monte Carlo and the probability for each amino acid at each position was 
determined (see Table 4). 



Table 4: Probability of amino acids at the designed positions resulting from the PDA technology 
calculation of the wild type (\A/T) enzyme structure. Only amino acids with a probability greater than 1 
% are shown. 



WT 


PDA™ Technology Probability Distribution 


5Y 


W7.2% 


F5.8% 


Y2.9% 


H4.0% 








7Q 


E9.1% 


LO.2% 












11 D 


11.2% 


DO.7% 


VO.1% 


M7.9% 


L6.4% 


E5.3% 


T4.2% 






Q3.8% 


Y2.6% 


F2.1% 


N1.9% 


S1.9% 


A1.1% 


37 V 


D9.9% 


M9.4% 


VI .4% 


S2.8% 


14.1% 


E1.0% 




39 G 


A9.8% 














63 N 


W1.2% 


Q6.7% 


A1.4% 










65 Y 


El. 7% 


L4.9% 


M3.4% 










67 T 


E1.0% 


D2.3% 


L3.9% 


A1.7% 








71 W 


V7.8% 


F5.5% 


W8.5% 


M6.0% 


D5.8% 


E4.3% 


11.0% 


BOY 


M2.4% 


LI. 5% 


F9.0% 


15.9% 


Y5.7% 


E3.7% 




82 V 


V8.6% 


D1.0% 
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88 Y 


N1.1% 


K6.6% 


W1.3% 


HOT 


09.9% 




115A 


A5.6% 


Y7.8% 


T4.4% 


118E 


E2.2% 


D2.6% 


. 12.0% 


125 F 


F9.4% 


Y1.8% 


M7.3% 


129 W 


E1.3% 


S8.6% 




168 V 


D8.1% 


A1.0% 




170 A 


A8.7% 


S7.6% 


D3.7% 



D10.2% 

A1.7% 

L1.5% 



S9.2% 



F2.6% 



If we consider aU the amino adds obtained from the PDA™ technology calculation, including those 
with probabilities less than 1%. we obtain 4.1 x 10« different amino add sequences. This is a 
reduction by 7 orders of magnitude. If one only considers ttiose amino adds that have at least a 
probability of more than 1% as shown in Table 1 (1% criterion), the problem is decreased to 6.6 x 10» 
sequences. If one neglects all amino acids with a probability of less than 5% (5% criterion) there are 
only 4.0 x 10^ sequences left. This is a number that may be easiV handled by screening and gene 
shuffling techniques. Increasing the list of low energy sequences to 100.000 does not change these 
nunr*ers significantly and the effect on the amino acids obtained at each position is negfigible. 
Changes occur only among the amino acids with a probability of less than 1%. Induding the 
substrate in the PDA™ fechnotogy calculabon further reduced the number of amino adds found at 
each position. If we consider those amino acids with a probabflity higher than 5%. we obtain 2.4 x 10^ 
sequences (see Table 5). 



JSilJ P""»f'>"*y °f amino acids at the designed positions resuHing from the PDA™ 
technology calculation of the enzyme substrate complex. Only those amino acids with a 
probability greaterthan1% are shown. -mmu acias wnn a 



Wr I PDA™ TECHI^OLOGY PROBABILITY DISTRIBUTION 



5Y 
70 

11D 

37V 

39G 

63N 

65Y 

67T 

71W 

SOY 

82V 

88Y 



Y69.2% W17.0% H7.3% F 6.0% 



078.1% 

D97.1% 
V50.9% 
S80.6% 
W92.2% 
E91.1% 
Ed2.8% 
W62.6% 
M66.4% 
V86.0% 

m5.i% 



E18.0% L 3.9% 

S5.4% A 1.2% 



□33.9% 
A19.4% 
D 3.9% 
L8.7% 
L5.2% 
E13.3% 
F13.6% 
D12.8% 
Y15.9% 



02.9% 



Ml 1.0% 
E10.7% 



S 6.9% 
1 6.0% 



L1.0% 



D4.0% 
LI .3% 



N11.4% F9.5% K1.9% Ql.4% D1.4% M1.4% 
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designed sequences. Wang Y. Zhang. H. Scott, RA. A new computational model for protein folding 
based on atomic solvation. Protein Sci. 1 995 Jul;4(7):1402- 1 1 .; Kuhlman B. Baicer, D. Native protein 
sequences are close to optimal for their structures. Proc Natl Acad Sci USA. 2000 Sep 12; 97(19): 
10383-8; Street AG, Datta, D, Gordon, DB, Mayo, SL. Designing protein beta-sheet surfaces by Z- 
score optimization. Phys Rev Lett. 2000 May 22; 84(21):5010-3.; Gordon, DB, Marshall, SA, Mayo. 
SL. Energy functions for protein design. Curr Opin Struct Biol. 1999 Aug; 9(4): 509-13. Review, 

In a preferred embodiment, one or more scoring functions are optimized or "trained" during the 
computational analysis, and then the analysis re-run using the optimized system. For example, the 
results of PDA™ technology calculations, described below, performed on decoy structures, may be 
used to obtain optimal sets of scoring functbn weights. First, the various components of a force field 
are factorized within the computer algorithm. A starting set of parameters, or weights, is defined 
based on best guess or previous knowledge of the parameter space. The current parameter set Is 
used in conjunction with a protein design algorithm to design one or more protein sequences and 
structures. These generated structures are then treated as decoy structures. The optimat set of 
parameters is considered to be that which predicts that decoys with properties very- different from the 
reference stmcture (native structure or prototype structure) are high in energy. In a preferred 
embodiment of the invention, a set of equations, relating the calculated energies of each decoy and 
comparison of each of its energy components to the reference, is used to iteratively optimize the 
parameters. 

In a prefenred embodiment of the Invention, the parameterization simulation begins with the creation 
of a number of decoy structures (e.g. 1 00-200) using random scoring function weights (within a 
predefined range), and a computatbnal protein design algorithm. 

In a preferred embodiment, parameters are modified at each iteration cycle according to the following 
equation: 

mod,,=A. ^\^f"^ P, . 

which details the modification of weight i based on evaluation of decoy d. E|.d and Ei,n represent the 
values of the ith scoring function component for the decoy and reference structures. Pd represents 
the normalized Boltzmann probability of the decoy structure according to its total energy using the 
current weights. The equation may be interpreted as follows. If a decoy structure's ith energy 
component is higher in value than that of the native (Ej ^ vs. the weight of the ith parameter is 
increased to an extent related to the difference in energy components. The extent of increase is 
further related to the current probability (P^) of the decoy structure - only high probability (low energy) 
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decoy structures contribute. The change in parameter will thus by definition lead to an increase in the 
calculated energy of the decoy relative to the reference (higher energy or score = less favorable). 
The value of X determines the rate at which the parameters are varied. The equation is applied to 
each decoy in the set Because the probabflities of the decoys are dynamically related to the change 
in parameters, multiple iterations over the current decoy set (see below) are applied: 



In a preferred embodiment, once the parameterization cycle begins, new decoys are generated with 
the modified parameters. This ensures broad enough coverage of parameter space by creating a 
broad range of decoy structures, and leads to the creation of ever-increasing competition between 
decoys and reference. However, because the actively created decoys win. at some point in the cycle, 
become high quality structures, an additional term may be added to the parameter modification 
scheme: 



In a prefen-ed embodiment, the parameterization is performed independently on a number of protein 
target structures. Parameterization using a number of small protein staictures has revealed, 
importantly, that optimal parameters derived fi-om one protein correlate strongly with those derived 
from different proteins. This result indicates that the invention yields parameter sets that are 
applicable to a wide variety of proteins. In an additional aspect the parameter optimization method is 
applied separately to sets of proteins that exhibit a common desired property (e.g. high solubility, 
thermostabinty). In this manner, force field parameters may be spedfically trained to design proteins 
with desired properties, such as thermostability, solubility, and the like. 

A key discovery associated with the method is that natural proteins appear to have highly conserved 
relative amounte of various energy components (polar group contecte. hydrogen bonding energy, 
ete.). Thus, although the invention is conveniently applied using the native stnjcture as a reference 
stete. analysis of multiple natural proteins has indicated that a prototypfeal protein may readily be 
defined and utilized as a reference state. 



Rnally. a diversity of related scoring function weights may be applied in separate applications of a 
protein design cyde. such that a diversity of sequence solutions are derived. 

In a preferred embodiment, the scoring funcBons outlined above may be biased or weighted in a 
variety of ways that does not invoh/e "training-. For example, a bias towards or away from a reference 
sequence or family of sequences may be incorporated; for exampte. a bias towards wild-type or 
homologue residues may be used. Similarly, the entire protein or a fragment thereof may be biased; 
for example, the active site may be biased towards wiW-type residues. A bias towards or against 
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increased energy may be generated. Furthemiore, biases may be used to design in selectivity. For 
example, a bias against sequences that bind to one or more unwanted substrates or receptors may 
be used. Additional scoring function biases include, but are not limited to applying electrostatic 
potential gradients or hydrophobicify gradients, and biasing towards a desired charge, isoelectric 
point, or hydrophobicity. In addition, experimental data, which may include values for any protein 
property or properties, may be used to generate biases or weights. 

Once the scoring functions to be used are identified for each variable position, the preferred first step 
in the computational analysis is the determinatbn of the interaction of each possible rotamer with all 
or part of the remainder of the protein. That is, the energy of interaction, as measured by one or more 
of the scoring functions, of each possible rotamer at each variable position (or each variable and 
floated position) with the backbone and/or other rotamers, is calculated. In a preferred embodiment, 
the interaction energy of each rotamer with the entire remainder of the protein, i.e. both the entire 
template and all other rotamers, is calculated. However, as outlined above, It is also possible to 
model only a portion of a protein, for example, a domain, motif, or site in a larger protein. 

In a preferred embodiment, two sets of interaction energies are calculated for each side chain rotamer 
at every position: the interaction energy between the rotamer and the template or backbone (the 
"singles" energy), and the interaction energy between the rotamer and all other possible rotamers at 
every other position (the "doubles" energy), whether that position is varied or floated. It should be 
understood that the template in this case includes both the atoms of the protein structure backbone, 
as well as the atoms of any fixed residues, as well as non-protein atoms in the scaffold. In an 
alternate embodiment, singles and doubles energies are calculated for fixed positions as well as for 
variable and floated positbns. 

Some energy terms, such as the secondary structure propensity scoring function, may be a 
component of the singles energy only. As will be appreciated by those in the art, many of the doubles 
energy terms will be close to zero, as many of the energy terms depend on the physical distance 
between the first rotamer and the second rotamer. That is, the farther apart the two moieties, the 
lower the energy typically will be. Furthermore, energy terms are not typically calculated for atoms 
that are separated by less than three, or alternatively less than four, covalent bonds. 

Once the singles and doubles energies are calculated and stored, the next step of the computational 
processing may occur: the identification of one or more sequences that have a low energy or 
favorable score. Alternatively, energies may be calculated as needed during the optimization steps, 
although this is often less computationally efficient. 

COflflBINATORIAL OPTIMIZATION ALGORITHMS 
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The discrete nature of rotamer sets allows a simple calculation of the number of possible rotameric 
sequences for a given design problem. A backbone of length n with m possible rotanners per position 
will have m" possible rotamer sequences, a number which grows exponentially with sequence length. 
For very simple design calculations, it is possible to examine each possible sequence in order to 
Identify the optimal sequence and/or one or more favorable sequences. However, for a typical design 
problem, the number of possible sequences (up to 10®** or more) is sufficiently terge that examinaUon 
of each possible sequence is intractable. A variety of combinatorial optimization algorithms may then 
be used to identify the optimum sequence and/or one or more favorable sequences. 

Combinatorial optimization algorithms may be divided into two classes: (1) those that are guaranteed 
to return the global minimum energy configuration if they converge, and (2) those that are not 
guaranteed to return the global minimum energy configuration, but which will always return a solution. 
Examples of the first class of algorithms include, but are not limited to, Dead- End Elimination (DEE) 
and Branch & Bound (B&B) (including Branch and Temiinate) (see Gordon and Mayo. Stmcture Fold. 
Des. 7:1089-98, 1999), and examples of the second dass of algorithms include, but are not limited to. 
Monte Cario (MC). self-consistent mean field (SCMF), Boltzmann sampling, simulated annealing, 
genetic algorithm (GA) and Fast and Accurate Side-Chain Topology and Energy Refinement 
(FASTER). 

Combinatorial optimization algorithms may be used alone or In conjunction with each other. 
Strategies for applying combinatorial optimization algorithms to protein design problems include, but 
are not limited to, (1) Find the global minimum energy configuration. (2) Find one or more low-energy 
or favorable sequences, and, most preferred. (3) Find the gtobal minimum energy configuration and 
then find one or more low-energy or favorable sequences. For example, as outlined in U.S.S.N. 
09/127,926 and POT US98/07254, preferred embodiments utilize a Dead End Elimination (DEE) step, 
and preferably a Monte Cario step. 

In a preferred enr*odlment when scoring is used, the primary library comprises the optimum 
sequence. That Is. computational processing is run until the simulation program converges on a 
singte sequence which is the global optimum. In a prefen-ed embodiment, the primary library 
comprises at least two optimized protein sequences. Thus for example, the computational processing 
step may eliminate a number of disfavored sequences but be stopped prior to convergence, providing 
a library of sequences of which the global optimum is one. In addition, further computational analysis, 
for example using a different method, may be run on the library, to further eliminate sequences or 
rank them differently. Alternatively, as is more fully described in U.S. Patent Nos. 6,269,312. 
6,188.965, and 6,403,312. which are herein expressly incorporated by reference, the global optimum 
may be reached, and then further computational processing may occur, which generates additional 
optimked sequences. 



46 



wo 03/014325 



PCT/US02/25588 



It should be noted that the preferred methods of the invention result In a rank-ordered list or filtered 
set of sequences; that is, the sequences are ranked or filtered on the basis of some objective criteria. 
However, as outlined herein, it is possible to create a set of non-ordered sequences, for example by 
generating a probability table directly that lists sequences without ranking them. 

In a preferred embodiment, an algorithm that is guaranteed to retum the global minimum free energy 
configuration (GMEC) Is used. However, such algorithms are not guaranteed to converge to a 
solution in a tractable amount of time. That is. the algorithm may get stuck. As a result, alternate 
strategies may be required for some design problems. 

The DEE calculation is based on the assumption that if the worst total interaction of a first rotamer is 
still better than the best total interaction of a second rotamer, then the second rotamer cannot be part 
of the global minimum energy configuration. An additional aspect of DEE states that if the energy of a 
rotamer sequence can always be lowered by changing from a first rotamer to a second rotamer, the 
first rotamer cannot be part of the global minimum. Since the energies of all rotamers have already 
been calculated, the DEE approach only requires sums over the sequence length to test and 
eliminate rotamers, which speeds up the calculations considerably. DEE may also include steps in 
which pairs of rotamers, or combinations of rotamers, are compared in order to identify sets of 
rotamers that are not compatible with the global minimum free energy configuration. In order to use 
DEE, the energy or scoring function must be pairwise-decomposable. That is, the energies or scores 
must be a function of the conformation and/or identity of at most two rotamers. 

In the B&B or A* algorithm, a tree is built, where a rotamer is first picked for one position, then a 
second position, and so on until one complete rotameric sequence is generated. The energy for that 
rotameric sequence is then calculated or obtained from the results of an earlier energy calculation. 
The process is then repeated, adding additional branches to the tree. However, in subsequent steps, 
if at any point the energy of the partially constructed rotameric sequence is worse than the energy of a 
previously identified complete rotameric sequence, all sequences that contain that partial rotameric 
sequence may be eliminated. The process may be completed until the GMEC Is identified. 
Additionally. B&B may be used to generate a list of all sequences that are within sonne energy or 
score of the GMEC (Gordon and Mayo, Stmcture Fold. Des. 7:1089-98, 1999) (Leach and Lemon, 
Proteins 33(2): 227-239, 1998). As for all the techniques listed herein, these algorithms can be used 
to generate a primary library or a secondary computational library. 

In an alternate embodiment, combinatorial search algorithms that are not guaranteed to return the 
GMEC may be used, either alone or following identification of the GMEC. These algorithms may also 
be referred to as sampling techniques. Algorithms that do not retum the GMEC are typically 
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computationally efficient and converge to a solution or solutions in a tractable, predictable amount of 
time. However, the quality of the solutions returned using these algorithms is variable, and may 
sometimes be insufficient These sampling methods may include the use of amino add substitutions, 
insertions or deletions, or recombinations of one or more sequences. 

Sampling techniques use a variety of approaches to jump between different points in sequence space 
(that is, between different possible variant sequences). For an sampling techniques, the kinds of 
altowable jumps may be altered (for example, jumps to random residues, jumps biased away from the 
vwld type sequence, jumps biased towards similar residues, jumps wher© multipie residue positions 
are simultaneously changed, etc). After jumping to a new sequence, the algorithm wfll choose 
whether to accept or reject the jump. The acceptance criteria for each sampling jump may be altered, 
by modifying the temperature factor. As vwll be appreciated by those stalled in the art high 
temperature fectors allow searches across a broad area of sequence space, and low temperature 
fectors allow searches over a narrow region of sequence space. See Metropolis et al.. J. Chem Phys 
v21, pp 1087, 1953, hereby e}^ressly incorporated by reference. 

A preferred embodiment utilizes a Monte Carto search, which is a series of biased, systematic, or 
random jumps in sequence space. Monte Carto searching may be used to explore sequence space 
around the global minimum, to find new local minima distant in sequence space . or to find one or 
more low energy sequences. A Monte Carto search may be performed to generate a ranl<-ordered list 
or filtered set of sequences in the neighborhood of the GMEC. Starting at the GMEC, random 
positions are changed to other resklues or rotamers (that is. the conformation and/or identity is 
changed), and ttie energy of tiie new sequence is calculated. If the new sequence meets the criteria 
for acceptance, it is used as a starting point for anotiier jump. After a predetermined number of 
jumps, a rank-ordered list or filtered set of sequences is generated. Monte Carto searches may also 
be started at sequences that are not the GMEC, including randomly selected sequences. Such 
searches may be used to generate a list of favorable sequences when ttie GMEC is not known. 

In another embodiment, self-consistent mean field CSCMF") metiiods (see Delame et al. Pac. Symp. 
Biocomput 109-21 (1997). Koehl et al., J. Mol. Biol. 239249 (1994); Koehl et al., Nat Stmc. Biol. 
2:163 (1995); Koehl et al.. Curr. Opin. Struct Biol. 6:222 (1996); Koehl etal., J. Mot. Bio. 293:1183 
(1999); Koehl et al., J. Mol. Biol 293:1161 (1999); Lee J. Mol. Biol. 236:918 (1994); and Vasquez 
Biopolymers 36:53-70 (1 995); all of which are expressly incorporated by reference.) are used. SCMF 
works by determining the optimal set of probabilities for all rotemer and residue states in the 
simulation, using a self-consistency criterion that relates ttie mean-field energies of the states to ttieir 
probabillHes. and vice versa. The final probabilities may be used to define a list of a favorable 
sequence combinations that define a combinatorial library of protein sequences. As for al the 
techniques listed herein, SCMF can be used to generate a primary library or a secondary 
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computational library. 

In a preferred embodiment, the sampling technique utilizes genetic algorithms, e.g., such as those 
described by Holland (Adaptation in Natural and Artificial Systems, 1975, Ann Arbor, U. Michigan 
Press). Genetic algorithm analysis generally takes generated sequences and recombines them 
computationally, similar to a nucleic acid recombination event, in a manner similar to gene shuffling . 
Thus the "Jumps" of genetic algorithm analysis generally are multiple position jumps. In addition, as 
outlined below, conrelated multiple jumps may also be done. Such jumps may occur with different 
crossover positions and more than one recombination at a time, and may involve recombination of 
two or more sequences. Furthermore, deletions or insertions (random or biased) may be done. In 
addition, as outlined below, genetic algorithm analysis may also be used after the secondary library 
has been generated. 

In a preferred embodiment. Boltzmann sampling is done. As will be appreciated by those in the art, 
the temperature factor criteria for Boltzmann sampling may be altered to allow broad searches at high 
temperature factors and narrow searches close to local optima at bw temperature tectors (see e.g., 
Metropolis et at, J. Chem. Phys. 21:1087. 1953). 

In a preferred embodiment, the sampling technique utilizes simulated annealing, e.g., such as 
described by Kirkpatrick et al. (Science, 220:671-680, 1983). Simulated annealing alters the cutoff for 
accepting good or bad jumps by altering the temperature factor in a systematk: manner. That is, 
slowly decreasing the temperature factor will slowly increase the stringency of the cutoff. This allows 
broad searches at high temperature factors to new areas of sequence space and narrow searches at 
low temperature factors to explore regions In detail. 

In a preferred embodiment, the FASTER method is used for determination of global optimization of 
the side chain conformations of proteins. The FASTER method focuses on resolving the combinatorial 
side chain packing problem, by converging on the near-optimal minima, (see Desmet. et al.. Proteins, 
48:31-43. 2002). 

TABOO 

In a prefenred embodiment, a diverse set of low-energy sequences is obtained using a class of 
algorithms referred to as tabu search algorithms. Traditionally, tabu search algorithms have been 
used to search for alternative local minima. The present invention presents a novel use of tabu 
search algorithms by using these algorithms to map amino acid sequence subspaces (see Modem 
Heuristic Search Methods, edited by V.J. Rayward-Smith, et al., 1996. John Wiley & Sons Ltd., 
hereby expressly incorporated by reference in its entirety). When used to map sequence space, the 
tabu search algorithms are referred to herein as "Taboo" search algorithms. A Taboo search 
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assumes that alternative optimization methods, such as protein design algorithms incorporating Dead 
End Elimination, genetic algorithms, Monte Carlo searches, have been used to provide the location of 
the global nrunlmum or a local minimum. Thus, a Taboo search is used for finding other regions or 
subspaces of the search space that contain local minima; preferably those that are reasonably low in 
energy compared Id the global minimum. 

Taboo searches are capable of Identifying alternative low energy basins because the search 
incorporates local optima avoidance by recording previously seen solutions by making a Bst of moves 
which have been made in the recent past of ttie search and which are tabu or forbidden for a certain 
number of iterations. That is, if a move in the search space has been made recentiy, that move is 
discouraged for some duration of time during the sampling procedure. The moves may be forbidden 
for some period of time or search (which can be varied), or weighted against but not forbidden. Such 
a mechanism helps to avoid qrcling and senses to promote the identification of alternative low energy 
basins. This concept Is illustrated in Figure 2. For example, by making the low energy basin 
identified by PDA ™ technology taboo, tfie search is forced to discover a different low energy basin. 
This cycle may be repeated until most or all of ttie alternative tow energy basins are identified. These 
alternative low energy basins or subspaces represent regions of the sequence space of a protein that 
should be explored experimentally by creation of secondary libraries (see Figure 3). 

Once a starting sequence is chosen, a taboo search is done. Preferably, tiie taboo search is done by 
applying one or more pseudo energies (pE) and serves to temporarily change the perceived energy 
landscape of ttie sequence space (see Figure 4). For example, if a single protein design simulation 
converges at iteration k to a variable protein sequence and structure ttiat contains amino add aa in 
rotamer state r at position i, tiien ttie matrix of side chain-template energies will be modified at 
iteratbn k+1 as follows: 

Where tiie pseudo energy at tiie very first iteration Is equivalent to the cateulated energy: 

PE°aa/J = Eaa,r,i 

and SEtaboo is defined by the simulation parameters. In some embodiments, ttie SEtaboo magnitude is 
dynamfc (e.g., random and/or slowly decreasing), again as defined by simulation parameters. Ega/j 
represents \he energy calculated by the force field or scoring function (i.e.. Ecaic in ttie Figures). 

Application of ttie pseudo energy increase discourages repeated conveigence to solutions containing 
amino acid aa in rotamer state r at position i in subsequent design simulations. However, if ttie 
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current value of pE is Insufficient to discourage incorporation of tlie rotamer, convergence to the same 
rotamer state in a subsequent sinaulation will result in additional increase of the pseudo energy. 

In a preferred embodiment, calculated and pseudo energies are stored in separate memory locations 
so that the calculated energy of any solutfon may be reported directly. This aspect is important for 
separating the effects of the taboo search from an accurate assessment of protein sequence 
energies. In a prefenred embodiment, the pseudo energy Increase Is applied to only one rotamer state 
of a converged aa at position i. In a preferred embodiment, the pseudo energy increase is applied to 
all rotamer states of a converged aa at position I In a preferred embodiment, the pseudo energy 
increase Is applied to a plurality of amino acid positions, and or a plurality of rotamer states. Thus, a 
taboo search results in the identification of altemate amino acids/rotamer states for at least one and 
preferably more than one amino acid position. Alternate amino acids/rotamer states may be reused in 
a protein design cycle to generate altemate variable protein sequences. 

In a preferred embodiment, the taboo search is done by applying a probability paranneter to at least 
one amino acid position. Such biased sampling is generally applied within so- called heuristic or 
stochastic sampling methods such as Monte Carlo or genetic algorithms. For example, many 
heuristic methods use the concept of Boltzmann sampling, which is based on the energy difference 
between two states. In a preferred embodiment, a probability parameter results in a modification of 
the Boltzmann probability (Pb) such that the sampling probability (Ps) is reduced: 

PB = e'<^^>/RT 

In a prefenred embodiment any or all of the methods described herein may utilize a recency 
parameter. In other words, application of a recency parameter ensures that the most recent moves in 
sequence space are prohibited for a certain number of iterations. Moves that are considered to be 
prohibited are derived from a mnning list, which is an ordered list of all moves performed throughout 
the search. If the length of the running list is limited, recency may be viewed as the equivalent of 
short term memory. As will be appreciated by those of skill in the art, one consequence of limiting the 
length of the running list is that the prohibited moves may be encouraged at a later point in the 
simulation to allow for the exploration of a sequence space that has not been visited for some defined 
duration. Thus, recency may be a fixed paranneter or allowed to vary dynamically during the search. 

In a preferred embodiment, recency is applied to the nr^odified energy matrix by continual application 
of a damping term to all pseudo energies as follows: 

PE'^aaAl = A*pEVr.» 
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This continual dannping has the simple effect of an exponential decay of the pseudo eneigy over 
multiple simulation cycles. 

In a preferred embodiment, the damping is applied at every simulation cyde before or after 
application of additional pseudo energy Increases. As will be appredated by those of skill in the art, 
this approach mathematically enforces a ceiling or upper limit on the magnitude of the pseudo energy, 
defined by the combination of A and 5E. 

In a prefen^ed embodiment, the application of the damping and pseudo energy temis are reversed. 
In an alternative embodiment, recency is applied to the modified energy maWx by continual 
application of a damping term to ail pseudo energies as follows: 

p^avj = A-y 

In a preferred embodiment, the frequency parameter is appRed such that the strengtti of the taboo 
energy increase is dependent on the number of times a given amino acid has occurred at a particular 
position. For example, ttie pseudo energy equation may be modified to include a frequency bias as 
follows: 

In other words, applied in ttiis manner the strength of tiie taboo energy increase depends on the 
frequency of occurrence (faa^j) of that amino acid or rotamer in previous solutions. In a preferred 
embodiment, the frequency parameter is biased against ttie most frequent amino acid residue at a 
particular position. Any or all of these metiiods involving recenqr and frequency parameters may be 
used rerteratively or combined in any order. 

Thus, taboo analysis can be done to generate sequences that are not the GMEC but are local rrnnima 
(low energy) as well. As for all the computational metiiods ouUined herein, this may be done at any 
point during ttie analysis. Thus, for example, a taboo analysis may be done to identify one or more 
starting scaffolds, e.g. even before a primary library is generated. Altematively. taboo analysis can be 
used as ttie computational analysis for primary and/or library. Altematively. taboo analysis can be 
applied in combination witti other computational techniques as either part of ttie primary or secondary 
Bbraiy generation. For example, taboo constraints may be added to a Monte Carlo search. 

SELECTIiyiG SEQUEWCgS F QR THE PRIMARY UBRARY 

In general, some subset of all possible sequences Is used as ttie primary library. However, in some 
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These calculafons show that PDA™ technology may significantly reduce the dimensionality of the 
problem and may bring it into the scope of gene shuffling and screening techniques (see Figure 3). 

Example 3 

Protocol for TNFa Llbraiy Expression and Purification 

Overnight culture preparation: 

Competent Tuner(DE3)pLysS cells in 96 well-PCR plates were transformed with 1 ul of TNFa library 
DIsiAs and spread on LB agar plates with 34 g/ml chloramphenicol and 100 jig/ml ampicfllin. After ai 
overnight growth at 37-C. a cotony was picked from each plate in 1.5 ml of CG media with 34 ^g/ml 
chloramphenicol and 100 ng/ml ampidllin kept in 96 deep well block. The btock was shaken at 250 
rpm at 37*Covemrght 

Expression: 

Colonies were picked from the plate into 5 ml CG media (34 ^ig/ml chloramphenicol and 100 ^g/m! 
ampiciilln) In 24-well block and grown at 37°C at 250 rpm until OD600 0.6 were reached, at which 
time IPTG was added to each well to l^iM concentration. The culture was grown 4 extra hours 
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Lysis: 

The 24-well block was centrifuged at 3000 rpm for 10 minutes. The pellets were resuspended in 700 
ul of lysis buffer (50 mM NaH2P04, 300 mM NaCI, 10 mM imidazole). After freezing at-80«»C for 20 
minutes and thawing at 37"C twice, MgCb was added to 10 mM, and DNase I to 75 ng/ml. The 
mixture was incubated at 37X for 30 minutes. 

Ni^* NTA column purification: 

Purification was earned out following Qiagen Ni NTA spin column purification protocol for native 
conditton. The purified protein was dialyzed against 1 X PBS for 1 hour at 4**C four times. Dialyzed 
protein was filter sterilized, using Millipore muitiscreenGV filter plate to allow the addition of protein to 
the sterile mammalian cell culture assay later on. 

Quantification: 

Purified protein was quantified by SDS PAGE, followed by Coomassie stain, and by Kodak digital 
image densitometry. 

TNFa Activitv Assay: 

The activity of variant TNFa protein samples were tested using Vybrant Assay Kit and Caspase Assay 
kit. Sytox Green nucleic acid stain is used to detect TNF-induced cell permeability in Actinomycin-D 
sensitized cell line. Upon binding to cellular nucleic acids, the stain exhibits a large fluorescence 
enhancement, which is then measured. This stain is excluded from live cells but penetrates cells with 
compromised membranes. 

Caspase assay is a fluorimetric assay, which may differentiate between apoptosis and necrosis in the 
cells. This kit measures the caspase activity, triggered during apoptosis of the cells. 

WEHI cells (Var-13 Cell Line from ATCC) were plated at 2.5 x 10® cells/mL. 24 hrs prior to ttie assay 
(100 ^Uwell for the Sytox assay and 50 ^L/well for the Caspase assay). 



Table 6. Activity Assay Results for Sytox vs. Caspase 
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Table 7. Activity Assay Results for Sytox 



i/VEHI TNF-a Activity Assay(%); 
(Sytox); 1,18.01 



Oligo Name 


Conc.(ng/ui) 


% Activity 






10 


100 1000 10000 


IM30Df 


51.40 


108 


100 99 


70 


K65Ef 


0.00 


85 


60 


31 


45 


G66Qf 


50.99 


90 


88 


57 


58 


QSTWf 


0.00 


94 


50 


27 


68 


K112Df 


14.06 


17 


15 


15 


42 


D143Ef 


35.86 


25 


15 


13 


21 


D143Nf 


38.62 


14 


13 


14 


26 


D143Qf 


19.28 


14 


13 


14 


18 


D143Sf 


45.91 


19 


14 


12 


75 


A145Rf 


41.18 


14 


11 


11 


15 


A145Kf 


56.47 


13 


12 


12 


81 


A145Ef 


53,29 


20 


13 


12 


64 


E146Kf 


33.28 


12 


12 


11 


17 


E146Rf 


32.22 


12 


11 


12 


21 


K65E/D143K 


B.39 


34 


18 


15 


19 


K65E/D143R 


9.97 


16 


15 


14 


81 


WT 53.92 


130 


117 112 


19 


A84V I 


27.84 


15 


15 


14 


33 



Table 8. Activity Assay Results Sytox vs. Caspase 



WEHi TNF^ Actr 
(Sytox); 02.01.01 


^ity Assa^ 




WEHI TNF-a AcUvity Assay(%); 
(Caspase); 02.01.01 


Oligo 
Name 


Cone, 
(ng/ul) 
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Neat" 


1:10 


1:100 


1:1000 


Neat** 


1:10 


1:100 


1:1000 . 


G66Q 


29.18 


122 


61 


30 


18 


90 


48 


35 


30 


D143N 


18.92 


13 


11 


13 


12 


37 


25 


25 


24 


D143Q 


0.00 


12 


13 


12 


12 


27 


23 


24 


27 


D143S 


13.39 


12 


12 


12 


12 


39 


26 


25 


25 


A145R 


39.62 


14 


13 


13 


12 


38 


23 


23 


23 


A145K 


23.67 


13 


11 


12 


13 


34 


24 


23 


22 


E146R 


34.11 


14 


14 


14 


13 


49 


31 


28 


28 


K65eD143K 


0.00 


12 


12 


13 


12 


31 


24 


24 


22 


K65E/D143R 


0.00 


12 


13 


14 


13 


32 


27 


26 


23 


WT 


0.00 


149 


116 


68 


32 


96 


80 


39 


27 



Neat**: Normalized to 500 ng/mL 



Figure 



Sytox Vs. Caspase; % Activitiy Based on Highest Std. Pt. 




♦ Sytox 012501 
■ Caspase 012501 
. Sytox 011801 
X Sytox 020101 
o Caspase 010201 



Data Analysis: Fluorescence vs. TNFa standard concentration was plotted to make a standard curve. 
Compare the fluorescence obtained from the highest point on the standard curve (5 ng/mL) to the 
fluorescence obtained from the unknown samples, to determine the % activity of the samples. The 
data may be analyzed using a four-parameter fit program to determine the 50% effective 
concentration for TNF (EC50). % Activity of unknown samples = (Fluor Of unknown samples/fluor. of 
5 ng/mL std. Point) x 100. 
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Example 4 

PDA™ Technology Calculations for soluble TNF-R (p55) 

Using pubUcly available protein three-dimensional structures for the p55 TNFR (Protein Data Bank 
codes 1ext. Incf, 1nr) both alone and complexed with its ilgand, PDA^'^ technology was used to 
design optimized soluble p55 receptors as TNF-a antagonists. For the library shown below, the 
sequences shown were generated using PDA™ technology relative to the Protein Data Bank 1 ext 
numbering scheme. Amino add residues known from the structure of the rcceptor-TNFa complex to 
be critical for p55 binding to TNFa were designed around. The results shown in Table 1 are an 
example of a library in which 15 positions from the wild -type p55 receptor were used for PDA™ 
technology design. Four of the positions chosen were nonpolar. 7 of the positions were charged, and 
4 were polar. The library shown in Table 1 was pooled from five independent designs, and a 1 5% 
cutoff was applied for each position in the library. The size of the library for a single mutation is 78 
and the entire library is 1.5 x 10^° sequences. The wild-type (\A^ sequence is shown in the first line 
of the table. The mutation pattem for soluble p55 receptors at given position is shown in the 
remainder of the table. 



Table 9. 
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EXAMPLE 5 

Computational Stabilization of Human Growtii Hormone (hGH) 

Human Growth Hormone (hGH) was computationally redesigned to improve its thermostability. The 
computational design was performed using a previously developed combinatorial optimizatbn 
algorithm based on the dead-end elimination theorem. The algorithm uses an empirical free energy 
function forni scoring designed sequences. This function was augmented with a tenm that accounte 
for the loss of backbone and side chain conformational entropy. The weighting factors for this term, 
the electrostatic interaction term, and the polar hydrogen burial term were optimized by minimizing the 
number of mutations designed by the algorithm relative to wild- type. Forty-five residues in the core of 
the protein were selected for optimization with the nfK)dified potential function. The proteins designed 
using the developed scoring function contained six to ten mutations, showed enhancement in the 
melting temperature of up to IB^C, and were biologically active in cell proliferation studies. (See 
Rlikov, et at. Computational stabilizatbn of human growth hormone. Protein Science (2002), 11: 
1452-1461, Cold Spring Harbor Laboratory Press, hereby expressly Incorporated by reference in its 
entirety.) 

EXAMPLE 6 

DEVELOPMENT OF A CYTOKINE ANALOG WITH ENHANCED STABILITY USING 
COMPUTATIONAL ULTRAHIGH THROUGHPUT SCREENING 

An ultra high throughput, computational screening method was used to Improve the physkx>-chemical 
characteristics of Granulocyte-Colony Stimulating Factor (G-CSF). Residues in the buried core were 
selected for optimization to minimize changes to the surface, thereby maintaining the active site and 
limiting the designed protein's potential for antigenicity. Using a structure that was homology modeled 
from bovine G-CSF, core designs of 25-34 residues were completed, con-espondlng to 10^^ - 10^® 
sequences screened. The optimal sequence fi^om each design was selected for biophysical 
characterization and experimental testing; each having 10-14 mutations. The designed proteins 
showed enhanced thermal stabilities of up to13 ^'C, displayed 5- to 10-fold improvements in shelf life, 
and were biologically active in cell proliferation assays and in a neutropenic mouse model. 
(See Luo, et at, Development of a cytokine analog with enhanced stability using ultrahigh 
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computational tliroughput screening, Protein Science (2002), 11: 1218-1226. Cold Spring Harbor 
Laljoratory Press, hereby expressly incorporated by reference in its entirety.) 
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WE CLAIM: 

1 . A method executed by a computer under the control of a program, said computer including a 
memory for storing said program, said method comprising the steps of: 

a) receiving a scaffold protein structure with 

residue positions; 

b) selecting a collection of variable residue 

positions from said residue positions; 

c) establishing a group of potential rotamers for each 

of said variable residue positions, and wherein a 
first group for a first variable residue position 
has a first set of rotamers from at least two 
diffierent amino acid side chains, and wherein a 
second group for a second variable residue 
position has a second set of rotamers from at 
least two different amino acid side chains; and, 

d) analyzing the interaction of each of said rotamers 

in each group with all or part of the remainder 
of said protein to generate a set of optimized 
protein sequences. 

2. A method according to claim 1 wherein said first and second sets of rotamers are different 

3. A method according to claim 1 wherein said first and second sets of rotamers are the sanrte. 

4. A method executed by a computer under the control of a program, said computer Including a 
memory for storing said program, said method comprising the steps of. 

a) receiving a scaffold protein with residue positions; 

b) selecting a collection of variable residue positions from said residue positions; 

c) establishing a group of potential amino acids for each of said variable residue positions, 
wherein a first group for a first variable residue position has a first set of at least two amino 
acid side chains, and wherein a second group for a second variable residue position has a 
second set of at least two different amino acid side chains; and. 

d) analyzing the Interaction of each of said amino acids with ait or part of the remainder of 
said protein to generate a set of optimized protein sequences. 

5. A method according to claims 1-4 wherein after step d) a library of said optimized protein 
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sequences is generated. 



6. A method according to claim 5 further comprising physicaBy generatang at least one member of 
said set of optimized protein sequences and experimentally testing said sequences for a desired 



function. 



method for generating a secondary library of scaffold protein variants comprising: 

a) providing a primary library comprising a filtered set of scaffold protein primary variant 
sequences; 

b) generating a list of primary variant positions in said primary library; 

c) combining a pluraBty of said primary variant positions to generate a secondary library of 
secondary sequences. 



8. A method for generaBng a secondary library of scaffold protein variants comprising: 

a) providing a primary library comprising a filtered set of scaffold protein primary variant 
sequences; 

b) generating a probability distribution of amino add residues In a plurality of variant positions; 

c) combining a plurality of said amino acid residues to generate a secondary Obrary of 
secondary sequences. 

9. A method according to claim 7 further comprising synthesizing a plurality of said secondary 
sequences. 



10. A method according to daim 8 wherein said synthesizing is done by muWpte PGR with pooled 
oligonucleotides. 

11. - A method according to daim 10 wherein said pooled ofigonucleotides are added In equimolar 
amounte. 



12. A method according to daim 10 wherein said pooled oligonucleotides are added in amounts 
that correspond to the firequency of the mutation. 

13. A composition comprising a plurality of secondary variant proteins comprising a subset of said 
secondary library. 

14. A composition comprising a plurality of nudeic acids encoding a plurality of secondary variant 
proteins comprising a subset of said secondary library. 
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15. A method for generating a secondary library of scaffold protein variants comprising: 

a) providing a first library rank-ordered list of scaffold protein primary variants; 

b) generating a probability distribution of amino acid residues in a plurality of variant positions; 

c) synthesizing a plurality of scaffold protein secondary variants comprising a plurality of said 
amino acid residues to fomia secondary library; 

wherein at least one of said secondary variants is different from said primary variants. 

16. A computational method comprising: 

a) receiving a scaffold protein with residue 

positions; 

b) selecting a collection of variable residue 

positions from said residue positions; 

a) providing a sequence alignment of a plurality 
of related proteins; 

b) generating frequencies of occurrence for 

individual amino acids in at least a plurality of positions with said alignments; 
e) creating a pseudo-energy scoring function 
using said frequencies; 

f) using said pseudo-energy scoring function and at feast one additional 
scoring function to generate a set of optimized protein sequences. 

17. A method according to claim 16 wherein said frequencies are weighted. 

18. A method according to claim 17 wherein said frequencies are weighted using a diversity 
weighting function. 

19. A method according to claim 17 wherein said frequencies are weighted using a sequence 
homology weighting function. 

20. A method according to claim 17 wherein said frequencies are weighted using a structural 
homology weighting function. 

21. A method according to claim 17 wherein said frequencies are weighted using a weighting function 
based on physical properties. 

22. A method according to claim 17 wherein said frequencies are weighted using a functional-based 
weighing function. 
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23. A method according to claims 19 or 20 wherein if said homology is high, said weighting is high. 

24. A method according to claim 19 or 20 wherein if said homotogy Is high, said weight is low. 

25. A method according to claim 17 wherein said multiple sequence alignment comprises proteins 
with related three-dimensional structures. 

26. A method according to claim 16 wherein pseudo-energy is based on bgarithms of said 
frequencies. 

27. A method according to claim 26 wherein said pseudo energy scoring function is based on log- 
odds ratios. 

28. A method according to claims 16-27 wherein after step f) a library of said optimized protein 
sequences is generated. 

29. A method according to claim 28 further comprising physically generating at least one member of 
said set of optimized protein sequences and experimentally testing said sequences for a desired 
function. 

30. A computational method comprising: 

a) receiving a scaffold protein with residue 
positions; 

b) selecting a collection of variable residue 
positions from said residue positions; 

c) providing a sequence alignment of a plurality 
of related proteins; 

d) generating a frequency of occurrence for 

individual amino acids in at least a plurality of positions with said proteins; 

e) selecting a group of potential amino acids for 

each of said variable residue positions, wherein a first group for a first variable residue 
position has a first set of at least two amino acid side chains, and wherein a second 
group for a second variable residue position has a second set of at least two different 
amino acid side chains according to their frequency of occurrence; and, 
0 analyzing the interaction of each of said 

amino acids at each variable residue position with all or part of the remainder of 
said protein using at least one scoring function to generate a set of optimized 
protein sequences. 
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31. A computational method according to claim 28, wherein amino acids with a frequency of 
occun^ence of at least 1 % are selected. 

32. A computational method according to claim 28, wherein amino acids with a frequency of 
occurrence of at least 5% are selected. 

33. A computational method according to claim 28, wherein amino acids with a frequency of 
occurrence of at least 10% are selected. 

34. A computational method according to claim 28, wherein amino acids with a frequency of 
occurrence of at least 20% are selected. 

35. A method according to claim 28 wherein said frequency is weighted. 

36. A method according to claim 28 wherein said frequencies are weighted using a diversity 
weighting function. 

37. A method according to claim 28 wherein said frequencies are weighted using a sequence 
homobgy weighting functbn. 

38. A method according to claim 28 wherein said frequencies are weighted using a structural 
homology weighting function. 

39. A method according to claim 28 wherein said frequencies are weighted using a weighting functbn 
based on physical properties. 

40. A method according to claim 28 wherein said frequencies are weighted using a functional-based 
weighing function. 

41. A method according to claims 37 or 38 wherein if said homology is high, said weighting is high. 

42. A method according to claims 37 or 38 wherein if said homology is high, said weight is low. 

43. A method according to claim 28 wherein said multiple sequence alignment comprises proteins 
with related three-dimensional stmctures. 

44. A computational method according claims 16-43 wherein said analyzing step further comprises at 
least two scoring functions. 
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45. A method according to claim 28 wherein said scoring function is selected from the group 
consisting of van der Waals potential scoring function, a hydrogen bond potential scoring function, an 
atomic solvation scoring function, an electrostatic scoring function, a secondary structure propensity 
scoring function and a pseudo-energy scoring function. 

46. A method according to claims 28^5 wherein after step f) a libraiy of said optimized protein 
sequences is generated. 

47. A method according to claim 28 further comprising physically generating at least one member of 
said set of optimized protein sequences and experimentally testing said sequences for a desired 
function. 

48. A computational method comprising: 

a) receiving a scaffold protein with residue 
positions; 

b) selecting a collection of variable residue 
positions from said residue positions; 

c) providing an amino acid substitution matrix; 

d) creating a pseudo-energy scoring function 
using said matrix; 

e) using said pseudo-eneigy scoring function and 

at (east one additional scoring function to generate a set of optimized protein 
sequences. 

49. A computational method according to claim 48 wherein said substitution matrix is selected from 
the group consisting of PAM, BLOSUM. and DAYHOFF. 

50. A method according to claims 46^9 wherein after step e) a library of said optimized protei n 
sequences is generated. 

51 . A method according to claim 50 further comprising physically generating at least one member of 
said set of optimized protein sequences and experimentally testing said sequences for a desired 
function. 

52. A method executed by a computer under the control of a program, said computer including a 
memory for storing said program, said method comprising the steps of: 

a) receiving a scaffold protein with residue 
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positions; 

b) selecting a collection of at least one 

variable residue position from said residue positions; 

c) importing a set of coordinates for a scaffold protein, said scaffold protein 
comprising amino acid positions; 

d) analyzing the interaction of each of said amino acids with all or part of the 
remainder of said protein; 

e) utilizing a plurality of scoring functions, at 

least a first a scoring function having a first weight and a second scoring function 
having a second weight, to generate at feast one variable decoy sequence; and, 

f) comparing the scores from said scoring functions of 

said variable decoy sequence to the scores of a reference state to generate modified 
weights, wherein each weight is increased if the corresponding score of the decoy is 
higher than the conresponding score of the reference state and each weight is 
decreased if the corresponding score of the decoy is lower than the corresponding score 
of the reference state and, wherein the extent of increase or decrease is based on the 
relative individual and total scores of the decoy and reference states. 

53. A method according to claim 52 comprising repeating steps a) and e) at least one or 
more times to generate a final modified weight for each scoring function. 

54. A method according to claim 52 wherein the collection of variable residue positions is modified 
repeating steps a through f). 

55. A method according to claim 52 wherein said final modified weight for each scoring function is 
used b generate a set of optimized protein sequences. 

56. A method according to daim 52, wherein the reference state is based on the native sequence 
and structure. 

57. A method according to claim 52. wherein the reference state is a prototypical protein, 

58. A method according to claim 52, wherein said prototypical protein Is derived from a set of proteins 
with similar physical or functional properties. 

59. A method according to claim 52. wherein said weights are optimized on a set of said scaffold 
proteins. 
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60. A method according to claim 52, wherein the extent of increase or decrease of said weights Is 
based on the total Boitzmann probabBities of said reference and decoy states. 

61. A method according to claim 52, wherein the extent of increase or decrease of said weights is 
based on the difference between individual scores of said decoy and reference states. 

62. A method according to daim 52 comprising replacing at least a single amino add in said scaffold 
protein to create a variable sequence and analyzing said variable sequence using said scoring 
fiincitons. 

63. A method according to claim 52 further comprising replacing a subset of amino adds, said subset 
selected from the group comprising core, boundary, and surface amino adds. 

64. A method according to dalm 52 further comprising protein design automation. 

65. A method according to daim 52 further comprising sequence prediction algorithm. 

66. A method according to clainre 52-65 wherein after step d) a libiary of said optimized protein 
sequences is generated. 

67. A method according to dalm 66 further comprising physically generating at le ast one member of 
said set of optimized protein sequences and experimentally testing said sequences for a desired 
function. 

68. A method executed by a computer under the control of a program, said computer induding a 
memory for storing said program, said method comprising the steps oft 

a) receiving a scaffold protein with residue positions; 

b) selecting a collection of variable residue positions from said residue positions; 

c) importing a set of coordinates for a scaffold protein, said scaffold protein comprising amino 
add positions; 

d) generating a variable protein sequence comprising a defined energy state for each amino 

acid position; 

e) applying an energy increase to at least one of said defined energy states tor a least one of 
said amino acid positions; and, 

0 generating at least one alternate variable protein sequence. 

69. A method executed by a computer under the control of a program, said computer including a 
memory for storing said program, said method comprising the steps of: 
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a) receiving a scaffold protein with residue positions; 

b) selecting a collection of variable residue positions from said residue positions; 

c) importing a set of coordinates for a scaffold protein, said scaffold protein comprising amino 
acid positions; 

d) generating a variable protein sequence comprising a defined energy state for each amino 
acid position; 

e) applying a probability parameter to at least one of said amino acid positions; and 

f) generating at least one alternate variable protein sequence. 

70. A method according to claim 68 or 69 wherein said energy increase is applied to a plurality of 
amino acid positions. 

71. A method according to claim 68 or 69 generating a plurality of alternate optimized variable protein 
sequences. 

72. A nnethod according to claim 68 further comprising applying a recency parameter with said 
energy increase. 

73. A metfiod according to claims 68-72 further comprises comparing said alternate optimized 
variable protein sequences. 

74. A method according to claim 68 further comprising applying a frequency parameter with said 
energy increase. 

75. A method according to claim 74 wherein said method comprises biasing said frequency 
parameter against the most frequent amino acid residue at a particular position. 

76. A method according to claim 68 wherein said energy increase includes the energy increase of a 
set of rotamers for at least one amino acid position. 

77. A method according to claim 68 wherein said energy increase includes the energy increase of a 
set of rotamers for a plurality of amino acid positions. 

78. A method according to claim 68 wherein said protein design cyde comprises applying protein 
design automation technology. 

79. A method according to claim 68 wherein said protein design cycle comprises applying the 
sequence prediction algorithm. 
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80. A method according to daim 68 wherein said protein design cyde comprises applying a force field 
calculabon. 

81. A method according to claims 68-80 wherein after step 0 a library of said optimized protein 
sequences is generated. 

82. A method according to claim 81 further comprising physically generating at least one member of 
said set of optimized protein sequences and experimentally testing said sequences for a desired 
function. 

83. A method executed by a computer under the control of a program, said computer including a 
memory for storing said program, said method comprising the steps of. 

a) receiving a scaffold protein with residue positions; 

b) selecting a collection of variable residue positions 
from said residue positions; 

c) importing a set of coordinates for a scaffold protein, 
said scaffold protein comprising amino add positions; 

d) generating a set of optimized variant protein sequences 
comprising one or more variant amino acids; and, 

e) applying a clustering algorithm to duster said set into 
a plurality of subsets. 

84. A method according to daim 83 further comprising applying a taboo search. 

85. A method according to claim 83 wherein said dustering algorithm comprises a single- linkage 
clustering algorithm. 

86. A method according to claim 83 wherein said dustering algorithm comprises a complete linkage 
dustering algorithm. 

87. A method according to claim 83 wherein said clustering algorithm comprises an average 
linkage clustering algorithm. 

88. A method according to daim 83 wherein saW subsets are dustered according to sequence 
similarity. 

89. A method according to daim 83 wherein said subsets are dustered according to energetic 
similarity. 
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90. A method according to claim 83 wherein DNA shufRing is applied with said subsets to generate a 
library of optimized protein sequences. 

91. A method according to daim 83 wherein said protein design cycle comprises protein design 
automation technology. 

92. A method according to clainn 83 wherein said protein design cycle comprises the sequence 
prediction algorithm. 

93. A method according to claim 83 wherein said protein design cycle comprises a force field 
calculation. 

94. A method according to claims 83 or 90 wherein said subsets are used to generate secondary 
libraries comprising related sequences. 

95. A method according to claims 83-90 wherein after step e) a library of said optimized protein 
sequences is generated. 

96. A method according to claims 94 or 95 further comprising physically generating at least one 
member of said set of optimized protein sequences and experimentally testing said sequences for a 
desired function. 

97. A method for identifying proteins that have a similar conformation to a target protein, said method 
comprising: 

a) receiving at least one scaffold protein 

structure with variable residue positions of a target protein; 

b) computationally generating a set of primary variant amino acid sequences that 
adopt a conformation similar to the conformation of said target protein; and, 

c) identifying at least one protein sequence that 

is similar to at least one member of said set of primary variants, but is dissimilar 
to said target protein amino acid sequence. 

98. A method according to claim 97, further comprising the step of confirming that said protein will 
adopt said conformation of said target protein. 

99. A method according to claim 97 wherein an amino acid sequence with less than 30% sequence 
identity is dissimilar. 
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1 00. A method according to claim 97 wherein an amino acid sequence with less than 20% sequence 
identity is dissimilar. 

101 . A method according to daim 97 wherein a similar conformation is a protein comprising a 
position for a given fold. 



.102. A method according to claim 97 wherein said computationally generating is applying a protein 
design algorithm. 

103. A method according to claim 102 wherein said computationally generating is applying protein 
design automation. 

104. A method according to claims 97 or 102 wherein said computationally generating step 
comprises a taboo search. 

105. A method according to daim 102 wherein said computationally generating step comprises 
applying a sequence prediction algorithm. 

106. A method according to claims 97 or 102 wherein said computationally generated sequences are 
used to create a Position Spedfic Scoring Matrix. 

107. A method according to claim 102 wherein said computationally generating indudes the use of at 
least two scoring functions. 

108. A method according to daim 107 wherein said scoring functions are seleded from the group 
consisting of a van der Waals potential scoring function, a hydrogen bond potential scoring function, 
an atomic solvaton scoring function, an electrostatic scoring function and a secondary stmdure 
propensity scoring fundion. 

109. A method according to daim 97 wherein the method for identifying said protein comprises 
searching pubfic databases. 

1 1 0. A method according to daim 97 wherein the method for identifying said protein comprises using 
a dynamic programming algorithm. 

11 1. A method according to daim 97 wherein said confimiing is seleded from the group consisting of 
x-ray crystallography, NMR spectroscopy, and combinations thereof. 
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112. A method for generating variant protein sequence libraries comprising: 

a) providing populations of at least two double stranded 

donor fragments corresponding to a nucleic acid template; 
• b) adding polymerase primers capable of hybridizing to end 
regions of each of said population of donor fragments; 

f) generating a population of hybrid double stranded 
molecules wherein one strand comprises a 5 -purification 
tag and the other strand comprises a 5'-phosphorylated 
overhang; 

g) enriching for variant strands by removing strands 
comprising a 5 -biotin moiety; 

h) annealing said variant strands to fornn at least two 
double stranded ligation substrates; and, 

1) ligating said ligation substrates to form a double 

stranded ligation product wherein said ligation product 
encodes a variant protein. 

113. A method according to claim 1 12 wherein one of said polymerase primers generates a variant 
nucleic acid strand. 

1 14. A method according to claim 1 12 wherein said template generates a variant nucleic acid strand. 

1 1 5. A method according to claim 1 1 2 wherein step e) precedes step d). 

1 16. A method according to claim 1 12 wherein steps a) through f) are repeated to generate a variant 
protein. 
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oecfUencG 


14 


0 


SGDNGWTDKESADSRYKFYFNDQTNESTSSMPTN 


18 


0 


SGPNGWMMMMLDSSimmrYDLLTNLSTSSEPAS 


54 


0 


SLPAGWTDMELDSSNYKYYFNMLYQQSTETMPTS 


66 


0 


SGDNGWTLMMLDSSNVMYYYNKETQLSTMSQPTS 


138 


0 


SLPNGWMLMELNSARYMYYYNMLYNLSTNSMPTN 


156 


0 


SGPNGWMLMEMDSSSYMFYYNMETQESTWSRPTS 








112 


1 


SGPAGVm^EMLQNADFWYTFNKIYEQMLEQMPTS 


119 


1 


NGPSGWTWKMMDSSNVWYTFNDSTNTSTMSMPTN 


127 


1 


sgpsgwmwmi^dnsqfyytfnnetntsleqmp™ 


129 


1 


SGPDGWTWKM^^ttroiU^YWYTFNl^ 


137 


1 


ngpsgwtwkmmdsadfwytfndatntmleempts 


141 


1 


sgppgwtwkmmdsadfwytynmatnemtesmpnd 


157 


1 


sgpsgwawkmldnaqywytfndlyqqstssmpnd 


169 


1 


sgpsgwmwmmaqsanfwftynlltntsleqmpts 


176 


1 


sgpdgwtwmmldsanfwytfnketqlaleqmptn 


182 


1 


sgpdgwmwkmmdsadfwft^ 


184 


1 


ngpdgwawmmlqnanywytynnetnestqempts 


186 


1 


SGPDGWTWKMMDNAQFWYTFNl^TNTAQESMPiro 


189 


1 


sgpdgwlwkmlqnanywytfndltqlatesmptn 


191 


1 


NGPSGWNTAMMDNAQFl/^nfTFNTIYQQ]^ 


193 


1 


sgppgwtwmmmndanvwftfnklteeatmsmptn 


196 


1 


SGPNGWTWMl^DSANYWYTYNMMYNLATASMPTS 


199 


1 


AGPSGWTWMMIJDSAlSrWFTYNMATNTATETMPTS 


203 


1 


SGPDGWTWMtMQSANYWYMFNMATQETTASMPTS 


208 


1 


SGPDGWAWMMLDNAQFWFMFNNETNTAQESMPTD 


211 


1 


NGPDGWMWMKLDNADYWYTFNLLTNTATETMPTN 


213 


1 


SGPSGWTWMMMQSANFWYTFNMATNTMLEQMPTS 


215 


1 


NGPPGWTV\ny[MLQSADYWTFNmSOT 


218 


1 


SGPSGWNTAny[MLQSAlSnrWFTFNTLYQOAQQEMPND 


220 


1 


NGPSGWTWKMMDSAQFWYTFNDATDTMQEQMPTS 


222 


1 


SGPDGWTWMKLDSADYWFTFNKLTNTAQASMPND 


224 


1 


NGPDGWlS^flJMM^D3NANFWYTFNLATNTAQEEM 


228 


1 


SGPSGWAWMMLDNANFWFTFNNETNTSQSEMPND 


235 


1 


SGPNGWMWMMLQNANFWYTFNNMYQQAQEQMPTN 


238 


1 


SGPDGWTWMMLDNAQFWFTFNMETNTAQESMPND 




1 


SGPDGWMWKMimNANYWYTFNNETNTMQAEMPTS 


242 


1 


SGPDGWMmMAQSADFWYTFNLATNTSQEQMPTS 


245 


1 


NGPrX5WTl^^^^lLDNANFWFTFNMATNTSQESMP^ 


248 


1 


SGPSGWTWMMMQSANFWYTFNMATQESTTSMPTS 


253 


1 


SGPDGWTWMMADNADYWFTFNKLTNTAQAEMPTS 


257 


1 


SGPPGWMWiynyiMQSANFWYTFNMATNESTEQMPTS 


259 


1 


NGPSGWTWMMLQSANYWYTFNMLYQQMQAEMPTN 


262 


1 


SGPI>3WTWMMLDNADFWFTFNKLTNTAQESMPND 


269 


1 


SGPPGWTWKMMDSANYWYTFNDLTNTAQASMPND 
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271 


1 


NGPSGWMWMMLDNANFWFTFNLATNTAQEQMPTS 


277 


1 


SGPDGWTWMMMQSANFWYTFNMATNTMQEQMPTN 


279 


1 


SGPSGWTWMMLQSANFWYTFNMATNTSQESMPND 


280 


1 


NGPDGVinyrWKMMDNADFOTTFNNETNTAQ 


281 


1 


SGPDGWTWMMLQSANYWYTFNMLTNTAQAEMPTS 


283 


1 


SGPDGWMWKMMDNANFWYTFlSrDATNTAQESM 


285 


1 


SGPDGVMVMMLQNANYWYTFNLATNTAQAEMPTS 


290 


1 


SGPDGWMWMMLQNANFWYTFNLLTNTAQESMPND 


299 


1 


SGPDGWMWlynyLLQNANFlftnfTFNLLTNTAQESMPT 


300 


1 


SGPDGWTWMMLQNANFWYTFNMLTNTAQESMPTS 








4 


2 


SGPDGWTWMMANSSQLQYTFNKETNTSTMENPTN 


8 


2 


NGPDGWMWMKLDSSSEEYTYNMATNESTKEEPTS 


10 


2 


AGPNGWTWLEANSSQYKFTYDMLYQQSTSSEPTS 


31 


2 


SLPSGWMWMMASDAQLEYTYNM^^ 


34 


2 


SGPNGWTWMKNNSSEYWYTYDNETQQSTESNPTN 


44 


2 


NGPDGWTWEMDQSSNRYYTFDKATEQATDTEPTS 


51 


2 


AGPPGWTWLMLDNSDVLYTYNLETQTATMEEPTS 


56 


2 


SGPDGWTWAKADSSELLYTYNNETQEATMSEPAS 


77 


2 


SGPSGWTWMESNDSQTLYMFNNETQQTTMTMPTS 


79 


2 


NLPSGWSWLKMDSASYEYTYNMSSNTATDSEPTS 


85 


2 


SGPSGWTWEMTNNSQDMYMFNKSTMTSTKSMPTS 


88 


2 


SGPAGWTMAMLQNADLRFAFNDATNTTTESDPAN 








1 , 


3 


NGPDGWTWl#lWQNASYEYTFlSnyiLTNTAQWSMPTS 


121 


3 


SGPPGWTWEiyMQNADYWYTFNKMTEFWLWEMPTS 


135 


3 


sgpsgwawmmwqsanywytynne™ 


139 


3 


SGPDGWTWEMWQNAQYYYTFNKMTNTMTWTMPTS 


143 


3 


SGPSGWTWMl^QSANYWYTFNKIYEQMTWTMPTS 


155 


3 


NG P S GWTWMMWQ S ANYWYTFNMATNTAQWEMPTS 


159 


3 


AGPSGWNVmiyn/tfQSANYYYVFNTIYQW 


161 


3 


SGPDGWTWMMMDSAQYWYTFl^LTDTiy^ 


163 


3 


NGPAGWLWM^QNANLWYTFNLATNTWTWSEPTS 


167 


3 


NGPDGWTWMKWNSADYYYTFNKATNTMTWSMPTS 


180 


3 


SGPPGWTWMMWDNANLWYTFNKATNTMTWTMPTS 


197 


3 


NGPSGwMWKMwENADYYYTFNDLTNTML 


201 


3 


SGPSGWTWMMWDSADYWYTFNKMYEQMQWSMPND 


206 


3 


SGPDGVraWKMVroSAimjYTFNDATNTWTWTMPTS 


230 


3 


SGPNGWMWMMWDNAQYWYTFNMMTNTMQWSQPTS 


232 


3 


SGPSGlAflyiWMMWDNANYWYTFNLETOT 


246 


3 


SGPSGWISFWMMWQSASYWYMFNTLYQQMTWTMPTS 


250 


3 


NGPNGWMWJyMWDSANYWYTra 


255 


3 


SGPNGWTWMMWQSANYWYTFimDT^ 


263 


3 


NGPSGWTWMMWDSAlvlYOTTFN^ 


265 


3 


SGPDGWTWEMWQSANYWYTFNKATNTMTWTMPTS 
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3 
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O Q o 

298 


3 


NGPIXSWMWMMWDNANYWYTFNLLTNTMQWSMPTS 








93 


4 


AGPNGWMWKMWDSASYEYTFNDQTDTATWSEPAS 


151 


4 


ngpixswtwmmmqnasveytfnmatnta™ 


187 


4 


NGPPGWTWEMWQSASYEYTFNKMTEQMLWEEPTS 








15 


5 


NGPTOWTMEMMQNSQFMYAFNKMYEQTQSSDPTS 


21 


5 


NGPPGWTMLMMQSAQVMYAYDMTTQQSTMSMPTN 


32 


5 


SGPIX5iAFLMMELDNSREMYAFNLLTQQSQWSMPTS 


47 


5 


SGPSGWTMAimDNADLWYSFNNETMLAQMENPTS 


52 


5 


SGPDGWLMMRMQSAETDYAFNDQTNETTDSMPTN 


c o 
68 


5 


SGPDGWTMAELSNSQFKYAFNNMSQQSQSAMPTS 


O 

82 


5 


SGPNGWTMAELDNSRVMYAFNNETQQSTMTMPTS 


91 


5 


NGPDGWNMMKMNSAEFMYAFNTl^^ 


100 


5 


NGPNGWLiynyiMADSANYMYAYNDMYlSnJAQSS^ 








2 


6 


SGPNGWmEMMDSSNVKYYFNKDTQQTTMTMPTS 


5 


6 


SLPSGWTMLMKQNADEWYYYNMMSNTSQWSEPAS 


7 


6 


NGPAGWADKEMADARVKFYFNDQTQLATMTMPTN 


23 


6 


SGPSGWADMELSSANYKYYFNNETNTTTDQMPTS 


28 


6 


SGPDGWTDLMLDNSEYRYYYNMLYNWSQSSMPTS 


30 


6 


NGPNGWTKEMMQSSDYWYYFNKASEQSTASEPAS 


33 


6 


AGPSGWTDAMLDSANFKFYFNDSTMTTTESEPTS 


39 


6 


SGPSGWTDLESENANYKYYYNLETQQSQSSMPAS 


43 


6 


DGPSGWNLMEKSDAREMYYFNNETDTSQWSKPND 


48 


6 


SGPDGWADMMAQSAI^RFYYNDMYQQSTSSMPTS 


58 


6 


SGPNGWTDMELDNSRYKYYFNMDTQQTTSSRPTS 


60 


6 


NGPDGWAMMMLQSLNYWYYFNDATNESTWSEPAS 


64 


6 


SGPDGWMNEEMQNSRFKYYFNKDTQQTQESMPTN 


74 


6 


NGPNGWTDLELDNARYKFYYNKATNLSQSSMPTS 


76 


6 


AGPDGWMNEMLQSAEVRYYYNALTNESTMSQPTS 


84 


6 


SGPDGWMDMMLOSANFKYYFNMLTEOTQSAMPND 


86 


6 


AGPNGWTKKKLDSAEYEYYYNMMSNESTWSEPTS 


89 


\j 


ovjir r^muiJriiSijjJooW VKi i YIvJjLioNTSTMSMPTN 


94 
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SGPEGWMDEMAQS SlSPf KYYYNKLYNL 


96 


6 


NGPPGWAMI^LNS SOFWF YFNNMTNESTESMPTN 


101 


6 


SLPDGWMDMEWNSSQYKYYFNKDTQTMTWTMPTN 


105 


6 


NLPPGWTDMEWPDARLKYYFNMMTNTAKWERPT55 


107 


6 


SGPSGWMIjynyiLDSANVMYYYl^TQESTMTMPTS 


109 


6 


DLPSGWLDMEKSNASEKYYFNLDTQQTTWSRPTN 


111 


6 


NGPDGWTLKQMDSSQYMYYFNDLTNTMQWAEPTS 


113 
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SLPSGWNDMELNSASYKYYFNTAYQLATSSMPTS 


118 


6 


SLPDGWMDMELDNASFKFYYNLLYNWAQESMPTS 
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134 
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bvjFWijW iL>JXlililjJJoAI:>YJ\x irJ^rLLxycjol ooKlr i o 


136 
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oVjFDUWMDMMMJjJNI o c* V Kr x r JM isxri y y 1 y WciN ir 1 o 


140 


6 


ACjjPAGWMNMRADbbNxDr YrNNbiyyAybbyFi b 


1 ii >i 
144 


6 


SGPNGWMIjKKMDbbEVMr x rNNEuNIjAiMbbFl b 


148 
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bli FDO W 1 UWJbAl^N b K X K X X r W isA A i*^ 


t C A 

154 


O 


SG PDGWLiDMEAbN Ab Y K. Y Y FN Liij^ y Jj lib brlF 1 b 


158 


6 


JMvaFlJoW 1 JsJ?iJ\JjNbAJJVlJx x r I>JJ\A1JM i 1 iiXll iiF i b 


160 


b 


bVjFJNIGWi>lUJ\rir(UJJNIbKr Jvr x xiNJJNIrtlNyAifibJyLFliNi 


162 


b 


DGPPGVWyiljMKijWbADYMY Yr 1 y bbJvFJMIJ 


168 


6 


DGFFGWADMbLiUJNIbKY K Y X r NJNM X yjjAy b bi\FJN 


179 


6 


NCjFDljWriUJXLbbbJJAKYrS.Y Yr J>JrsJJiyi:»Ax AbKF^I b 


lot 

181 


6 


S G P b WN KMiMLMy b Ab V r. Y 


185 


6 


AGPSGWMLiMKADSADYTlYYFNLL YyyTl bby F 1 ^ 


188 


6 


SGPNGWMDMEMDS S S VKF YFNNETNES 


200 


6 


SGPDGWADMESSNAyYKYYFNLLTQEbTbbMPTb 


209 


6 


o^T^o/^T»TrnT^Airn» A ATO A OTT v^T'vrTr'ATA/i'A rn/^mcfn A OAif Turner 

SGPSGWTDMEANSASYKYYFNMATyTSTAbMPTS 


223 


6 


SG P SGWTD AES SDAQYK YY FNNETQ E ST S Sy PTS 


225 


6 


AGPPGlfMDMELQ SAS YKYi FNMMYEyTTS SRPTS 


229 


6 


NG PDGWTDEEMy SARYKYY FNKljTEy bTb bMFTb 








3 


7 


SGPSGWLNMEIjDS AS FMFMYWLiATNE b TE SMFTb 


6 


7 


SGPPGWNWDlKKDbbEYbYMrNlLYEy i lAbMFib 


12 


7 


NGPDGWMKJMMMDN bb VEYMYNJjAlXs|L(b i M IMF 1 b 


19 


7 


(-* T^T'x^r.TKTATmm A ATT~> A T^\n^\/TUf 'CATATATm/*\'C* CTI A OAAT*lATT*^ 

SGPDGwNDJKE<AIvrUADYDYMFIm i AbMFJMU 


20 


7 


SGPSGWTWKEMDNSRFLiYTFNDSSNTAyESNPTN 


24 


7 


SGPDGWTWMKLDNSEFQ YroYTJKLYlslJjbySbJyiF i b 


25 


7 


AT/^T%C/^«TKifA.113*A>f A C* A T^X/TiTl/TTIK/rMrn A rn'C(|MillHL,»OC*T> A O 

NGPbGWMNEMAybADYWYTrlNTATE ill IsbEFAb 


35 


7 


SGPPGWTNERMySADVDYWFNKASNEATMTMFWL) 


36 


7 


NG PDGWNwMMADNS S YEFTFNTLYy LiTyAERFTD 


37 


7 


SG PSGwMDMML PDAQ FMYvYtJMMTNTSTS^ 




/ 


iNwIrl.AjVvAX\XTilTUT1^00 O V d X ITIX^ IN^A X \^1JA X 1*1 J. I'lf X O 


46 


7 


SGPNGWTWMELNSSSFLYTFNLLTQESTTSNPAN 


50 


7 


SGPEGWTNEESNDAQFKYMFNKLYQWSQSSMPND 


71 


7 


SGPSGWTOKKADNADYWYTYNDLYNWAQAENPAN 


75 


7 


SGPPGWLWKMADSSNFWYTFNDMYKFALEQMPND 


80 


7 


SGPDGWMKiyQyiLQSSNYWFMyNLLTNESTASMPAS 


83 


7 


NGPSGWAWMKANSADYLFTYNDAYNLATESEPTS 


92 


7 


SGPSGWTNMELDNSSYKFMYNKETNLSTASMPTS 


98 


7 


SGPDGWTDEEMQSASYKYMFNKLTEQTTSSMPTS 


99 


7 


NGPSGWTWMKLDNSEFYFTFNMATNQSTSEEPTS 
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-L W _) 


7 


bCjFSGWNNMRLQSAJOT^DYMYNLLYNTST^ 


JL \J \J 


7 


bu FNoWAVmjyujQNS QFEYTFI^ 


X X o 


7 


bCyir'JNIijWAKMKl^SADYDYMFNNMTQQSTASMPND 




7 


SGPSGWTWbOKMDSSDVWFTYNLIYNLATMSMPTS 


J. 


7 


SGPEGWTNEMLiNDSQYWYMFNKIYEQTQASMPND 


J. Oft 


7 


SGPDGWTDMEMDS S S VKFMYNMLTNEATMTMPTN 


J. / J 


7 


SGPDGWTDDffiliNSSSYKFMYNMLTNESTESMPTN 


1 7ft 


7 


SGPNGWMWMMADSSQYYFTYNLMyNQSQSEMPND 








Q 


Q 


SGPSGWTMEMLQNANFMYAFNKETDTTQESMPND 




Q 
O 


SGPPGWNwMMLQSANFWYTFNTLTQTSQEQNPAN 




O 

o 


bGFbGwTWMKLESAEYLYTYNMESQFMLWEMPTS 


1 7 
X / 


p 


bUFbGWAWMMMDSASVEFTFNNATNTATMTQPTS 


97 


o 


AGPAGWTWKMMDSANVYYTFNDQTNEATMSQPTS 


^ 27 


p 
o 


bOFbbWNV\ny[ESNSASFLFTFNNETDTSLEQMPTN 


f4 X 


p 

o 


AGPPGWMWKMLDSADFWFTFNNLSNTSLEQQPTN 




p 


SGPNGWTDMKWDNSEYRYiyrrrN^MTNESTWSMPTS 


f± -7 


p 

O 


NGPSGWTWMKKDSSEEYYTFNMASNTSTWSEPTS 




p 


NGPSGWNNMMMENAQFWYMFNTATNTTTESMPTS 




p 


SGPPGWMWMKNNSASFEFTYNMLTNFALEQMPTS 


D X 


p 
O 


AGPSGWMWMEMDNASVLYTFNLESMTSQMEQPTS 




p 

o 


bGFSGWNWMMLNDAQEYFTFNTLTNTATKENPTN 




p 


NGPAGWMWKKWDSADYWYTFNDATDTSLWQMPTS 


Q -7 


p 
o 


SGPPGWNVMIKQSANELYTFNLATNTSTWTQPTS 


79 


p 


btjFiJbWAMMMLQSASYMYAFNNETQQSTWSMPTS 


ft7 


p 

25 


NGPDGWNWMMMDSSSEWYTFDNETQQAQKEMPTS 




p 
O 


SGPSGWTWMMKSDAQELYTFNMDTQEATWTMPTS 


-7 J 


O 

o 


bCjFJJLiW i WMKMDNADVDYTFDMASQTTTMTMPTS 


-7 / 


p 


SGPSGWTWLMLNDASDEYTYNLETQTAQKSMPND 


1 09 

X 


p 


SGPAGWTWKMMENADFWYTFNNETNEMTEQEPAS 


1 Oil 


p 
o 


AGPDGWTWMMMDSS S YWFTFNKATEQSQWSNPND 


1 OP 
X\J o 


p 
o 


SG PDG WTWKKMDS S DF WYTFNMATNTAQ SEMPTN 


Tin 

xxu 


p 

o 


SGPNGWMWMMWNSANYYFTYNNMTNQWTWSMPND 


117 
XX / 


p 

o 


SGPSGWTWEMWQSLQYYYTFNKASNTMTWTQPTS 


1 9'^ 
X^ J 


o 

o 


SGPAGWMWMKMNSASYEYTFNMLSNTSTETMPTS 


131 


8 


AG P SG WLWMMLOS ADYW7TYNT .'nTKFrQ'Pi?'T»ni>nic 


132 


8 


NGPDGWTNiyOCWESAEYMYMFNKATEQSTWSMPTN 


133 


8 


SGPPGVmVi«CMMDNSQFYYTFNDMTNTMQEQNPTN 


145 


8 


SGPDGWTWMMMPNASFEYTFNMATQTSLEQMPTS 


147 


8 


NGPPGWNWMMLOSAim^TYNMMTNEATASMPTS 


149 


8 


SGPSGWLWMKNNSADFWFTYNLETNTSQEEMPTS 


153 


8 


NGPPGWMWKKMDSADFYFTFNNMTNTSTEQEPTS 


165 


8 


NGPSGWMWKMKDSAQYYYTFNDSTNTEQAEEPTS 


171 


8 


NGPSGWTWiyttmDNASVEYTFNKETNTTTMTMPTS 


174 


8 


NGPAGWNWMMMENADYWYTFNTATEQSTOEEPTS 
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194 


8 


NG PNGWAWMKAQNADF Y YVFNNETQ QTTESMPT S 








11 


9 


SGPSGWLDMRWDSADYl^YFNNMYNTATWARPTD 


78 


9 


SGPNGWNDMMWPNADYKYYFNTATETSQWENPTS 


115 


9 
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These two views of representative structures from cluster 1, 3, and 9 illustrate 
the relationship between clusters and couplings of amino acids. In the upper 
set, it is clear that clusters 1 and 3 share a coupling pattern that is distinct from 
that of cluster 9. However, the bottom set reveals that on the opposite side of 
the protein, cluster 1 and 3 have a different coupling pattern. 
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BOX n. OBSERVATIONS WBDEKE UNITY OF INVEiSTION IS LACKING 

This application contains the following inventions or groups of inventions which are not so linked as to fonn a single general 
inventive concept under PCX Rule 13.1. In order for all inventions to be examined, the appropriate additional examination fees must 

be paid. 

Group I, claim(s) 1-3, 5-6, drawn to a method of determining an optimum product enoploying a conqjuter program wherein a group of 
potential rotamers are determmed or calculated. 

Group li, claim(s) 4-6, drawn to a method for generating a set of optimized sequences executed by a con^uter program vdierein a 
potential amino acids in the variable positions of the sequence is determined. 

Group in, claim(s) 7,9, drawn to a method for generating a secondary library by generating a list of prunary variant positions. 

Group IV, claini(s) 8, 10-12, drawn to a method for generating a secondary library by generatmg a probability distribution of amino 
acids. 

Group V, claim(s) 13, drawn to a con?)osition corcprising a plurality of secondary variant proteins. 
Group VI, claim(s) 14, drawn to a composition comprising of a plurality of nucleic acid. 

Group Vn, claim(s) 15, dravm to a method for generating a secondary library of scaffold protem by using a first library rank-ordered 
list of scaffdd protein primary variants. 

Group VIII, claim(s) 16-29, 31-47, 50-51, drawn to a method of generatmg an optimized protein sequence executed by a con^ 
program using sequence alignment of the proteins and a pseudo-energy ftinction calculation or steps. 

Group IX. claim(s) 30, 44, 46, 50-5 1 , drawn to a method of generathig an optimum protein by a computational method without using 
a computer program employing numerous process steps such as sequence alignment of related proteins, frequency of occurrence and 
other steps as recited. 

Group X, ciaim(s) 48-51, drawn to a metliod of generating optimized protein sequences using (manual) computational method ushag 
substitution matrix. 

Group XI, claun(s) 52-67, drawn to a method executed by a conq)uter program including the step of using coordinates for the scaffold 
protein and scoring fiinctions to generate a decoy sequence and comparative scoring. 

Group XII, claim(s) 68, 70-82, drawn to a method of generating an alternate variable protein sequence using a con5)Uter program 
including the step of applying a defined energy state to a variable amino acid position of the protein scaffold. 

Group XIII, claim(s) 69-71 , 73, 81-82, drawn to a method of generating a variable protein sequence executed by computer program 

including the step of applying a probability parameter. 

Group XIV, claim(s) 83-96, drawn to a method executed by a computer program wherein a set of coordinates are imported into a 
scaffold protein and employing a clustering algorithm. 

Group XV, claim(s) 97- 1 1 1, drawn to a method of identifying a protem with a simflar confonnation as the target protein. 

Group XVI, claim(s) 1 13-1 16, drawn to a method of generating a variant protein sequence libraries by polymerize reaction. 

The inventions listed as Groups l-II, VIII, XI, XII, Xni, XIV do not relate to a single general inventive concept under PCT 
Rule 13.1 because, under PCT Rule 13.2, they lack the same or corresponding special technical features for the following reasons: 
each of these methods would eniploy different programs since different parameters are required to achieve an optimum protein 
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prodact For example, the method of Group I would require a program to calculate the rotamers that would lead to an optimam 
product The method of Group H requires a program for the presence or absence of an amino acid residue in a scaffold protein 
position. Thus, these methods vwjuld lack the same q)ecial technical features of a resultant prodict and/or method steps. 



The inventions listed as Groups in, IV and VII do not relate to a single general inventive concq>t under PCX Rule 13, 1 
because, under PCT Rule 13.2, they lack the same or corresponding special technical features for tiie following reasons: the method 
of GroTip in generates a secondary library en5)loying different method steps from those of Groups IV and VH. Thus, there are no 
corresponding special technical features for flie different method stg}s resulting in differaxt secondary library compositions. 

The inventions listed as Groups V and VI do not relate to a single general inventive concept under PCT Rule 13.1 because, 
under PCT Rule 13.2, they lack the same or corresponding special technical features for the following reasons: Group V is drawn to a 
con^josition conprising of a phirality of variant proteins that lack the technical features of Group VI composition relating to a 
structurally different phiral nucleic acid conspounds Because of the degeneracy of the nucleic acid triplet differcait protein products 
may be encoded by the nucleic adds. . 

The inventions listed as Groups IX and X do not relate to a single general inventive concept under PCT Rule 13. 1 because, 
under PCT Rule 13.2, fhey lack the same or correspondmg special technical features for the foUowmg reasons: Gronps IX en5)loys * 
different process steps such as alignment and frequency of occurrence \*ich does not conespond to the process step of Group Xnsmg 
substitution of residues in the scaffold protein region. 



The inventions listed as Gronps XV and XVI do not relate to a single general inventive concept under PCT Rule 13. 1 because, under 
PCT Rule 13.2, they lack the same or corresponding special technical features for the following reasons: Groups XV and XVI are 
drawn to two differait methods resulting in two different products. Group XVI uses polymerase reaction while Gronp XV employs 
the conformation of a target that would generate two different libraries. 

The inventions listed as Groups (V. VI) and (MV and VII-XVI) do not relate to a single general inventive conceptunder PCT Rule 
13.1 because, under PCT Rule 13.2, they lack the same or corresponding special technical features for the foUowing reasons: Groups 
V and VI are drawn to CQnq)ositions vAach lack (he technical features of the processes of Groups I-IV and VII-XVI. 



The inventions listed as Groups I-IV and VII-XVI do not relate to a smgle general inventive concept under PCT Rule 13. i because, 
under PCT Rule 13.2, they lack the same or correspondmg special technical features for the foUowmg reasons: these groi5)s are 
drawn to different methods resulting in different optimized protein or Ubraries. For example. Group I en5)loys con?)uter to generate 
proteins. Group III, for example, relate to a non-con5)uter method of geieratmg secondary libraries. (See tiie above groupings of the 
different method Groups). 



This application contams dauns directed to more than one ^cies of the generic mvention. These species are deemed to lack unity of 
invention because they are not so linked as to form a single general mventive concept under PCT Rule 13, 1 . 

In order for more than one species to be examined, flie apprqpriate additional exammation fees must be paid. The species are as 
follows: 

Species: 

A) . Frequency Calculation (claims 17-22). 

B) . Scoring function (chum 45). 

C) . Subset (claun 63). 

D) . Identifying means (clauns 109 or 1 10). 

The species listed above do not relate to a single general inventive concept under PCT Rule 13. 1 because, under PCT Rule 13.2, the 
species lack the same or corresponding special technical features for tiie following reasons: each of the species recited in each 
subgroups A-D do not have the same or corresponding technical features as each of tiie species in e.g.. Group A involves different 
steps of calculating tiie frequency occurrence of an amino acid residue in tiie scaffold position. For example, ttie metiiod by which a 

homologous amino acid re sidue in one or more position of tiie scaffold protem would not detemune or relate to the mefliod of 
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determimng the structure (i.e. , conformatioBal) or function of the variable protein. Likewise, the species recited in subgroup D 
wherein a protein is identified by searching the pubUc databases would lack the special technical features of identifying said protein 
xisiog dj^namic conipiiter programming algorittuiL 
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